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The  growing  use  of  single  mode  fibers  in  local-area 
networks  and  customer  premises  networks  will  increase  the  need 
for  passive  optical  components,  such  as  branching  devices, 
mixers  etc.  Integrated  optical  devices  are  potentially  ideal 
for  these  applications,  provided  that  they  can  be  made 
compatible  with  single  mode  fibers.  The  use  of  Ge02  as  the 
core  dopant  and  Si02  as  the  substrate  ensures  that  these 
waveguides  will  have  virtually  identical  characteristics  to 
single  mode  fibers.  Additionally,  glasses  in  the  form  of 
waveguides  have  recently  been  used  to  study  various  nonlinear 
optical  phenomena,  which  provide  great  potential  applications 
such  as  data  storage  and  information  processing. 

The  present  study  has  for  the  first  time  demonstrated  the 
feasibility  of  employing  both  sol-gel  multiple  dip-coating  and 
low  pressure  chemical  vapor  deposition (LPCVD)  in  the  produc- 
tion of  Ge02-Si02  waveguiding  films  with  various  germania 
contents.  The  thin  film  characteristics  were  studied  by 
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various  analytical  techniques (e.g.  ellipsometry , waveguiding 
Raman  spectroscopy,  FTIR,  XPS,  SEM/TEM) . The  composition 
dependence  of  the  linear  refractive  index  of  Ge02-Si02  films 
follows  that  predicted  by  the  Lorentz-Lorenz  model.  Vibra- 
tional spectroscopy  revealed  the  existence  of  Si-O-Ge  linkages 
in  Ge02-Si02  glass  network.  The  addition  of  Ge02  in  Si02  caused 
a decrease  in  the  size  of  both  the  D1  and  D2  defect  bands  in 
the  Si02  Raman  spectra.  The  structure  of  the  LPCVD  film 
appears  to  be  dominated  by  D1  and  D2  defect  bands. 

Using  a three-prism  loss  measurement  technique,  the 
propagation  losses  were  found  to  be  3.31  dB/cm  and  2.59dB/cm 
for  sol-gel  and  LPCVD  films,  respectively.  These  losses  are 
attributed  to  various  scattering  processes  in  the  films.  The 
mode  indices  of  the  waveguide  were  measured  using  a prism 
coupling  technique.  The  measured  mode  indices  were  found  to 
agree  with  the  calculated  value  based  upon  a step-index 
profile  assumption.  The  theoretical  electromagnetic  field 
distribution  profiles  for  a step-index  planar  waveguide  has 
been  calculated  and  compared  to  the  experimentally  measured 
mode  profiles  using  a near  field  technique. 

The  nonlinear  refractive  indices  of  the  sol-gel  films 
(Ge02-Si02,  and  Ge02-Ti02)  were  measured  using  a third-harmonic- 
generation  (THG)  interferometry  fringe  technique.  The  relation 
between  n2THG  and  n,  was  found  to  follow  that  predicted  by  the 
empirical  BGO  model.  An  additive  model  was  used  to  calculate 
the  linear  refractive  indices,  Abbe  numbers,  and  n,  dispersion 
curves  of  the  films. 
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CHAPTER  I 


INTRODUCTION 


1.1  Importance  of  Ge02-Si02  Waveguides 

1.1.1  Materials  for  Optical  Fibers 

The  transmission  and  processing  of  signals  carried  by 
optical  beams  rather  than  by  electrical  currents  or  radio 
waves  has  been  a topic  of  great  interest  ever  since  the  early 
1960s,  when  the  development  of  the  laser  first  provided  a 
stable  source  of  coherent  light  for  such  applications.  There 
are  many  aspects  required  for  optical  fiber  materials.  The 
major  features  to  be  considered  are 

(1) .  Optical  loss  properties:  intrinsic  loss  properties 

and  potentiality  of  extrinsic  loss  reduction  are  of 
vital  importance. 

(2)  . Refractive  index  control:  precise  control  of 

refractive  index  profile  in  radial  direction  and 
the  reduction  of  index  fluctuation  along  the  axial 
direction  are  necessary. 

(3)  . Shape  control:  cross-sectional  shape  and  size,  the 
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surface  finish  and  the  fluctuation  of  the  cross- 
sectional  dimensions  along  axial  direction  of 
fibers  strongly  affect  the  transmission 
characteristics  of  the  fibers. 

Laser  beams  can  be  transmitted  through  the  air,  but 
atmospheric  variation  causes  undesirable  changes  in  the 
optical  characteristics  of  the  path  from  time  to  time.  From 
the  historical  view  point,  many  kinds  of  light-transmission 
waveguides  were  studied  for  their  practical  applications. 
These  materials  cover  gas,  liquid,  crystal  and  glass.  Many 
attempts  to  produce  gas  waveguides  have  been  made  [Ber64, 
Eag62].  The  fatal  defect  in  the  use  of  gaseous  materials  is 
the  difficulty  of  controlling  the  refractive  index. 

Some  liquids  are  also  low  loss  media  for  optical 
waveguides  [Mar64,  Ston72].  However,  the  refractive  index  of 
a liquid  changes  considerably  so  that  it  is  fairly  hard  to 
maintain  the  precise  waveguide  characteristics.  Furthermore, 
only  step-index  waveguides  can  be  fabricated  for  practical 
use. 

Solid  state  materials  have  stronger  interaction  with 
light  waves  than  gases  and  liquids.  Therefore,  in  general, 
the  transmission  loss  may  be  higher;  however,  the  optical 
properties  are  much  more  stable  compared  to  gases  and  liquids. 

Solid  state  materials  have  strong  absorption  at 
ultraviolet  wavelength  and  infrared  wavelength  regions,  which 
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are  attributed  to  electronic  transitions  and  molecular 
vibrations,  respectively.  These  absorption  peaks  affect  the 
absorption  loss  in  the  visible  and  near-infrared  wavelength 
regions.  In  order  to  find  low-loss  materials,  it  is  better  to 
select  materials  whose  strong  absorption  wavelengths  are  far 
from  each  other, i.e.  having  a wider  transmission  window. 

Solid  state  materials  are  classified  into  two  groups: 
crystalline  state  and  amorphous  state.  The  optical  characte- 
ristics of  the  two  states  are  quite  different,  even  if  the 
atomic  composition  is  the  same.  Crystalline  materials,  such 
as  quartz,  sapphire,  and  alkaline  halides  have  essentially 
very  low  loss  compared  to  amorphous  materials,  because 
Rayleigh  scattering  loss,  which  originates  from  the  density 
fluctuations,  is  less  than  one-tenth  that  of  amorphous 
materials.  Therefore,  it  can  be  said  that  the  best  materials 
may  be  found  in  single  crystals  from  the  viewpoint  of 
intrinsic  loss.  However,  high  speed  crystal  growth  or  mass 
production  of  single  crystal  fibers  is  very  difficult  [TangSO, 
Bri80,  Mim80,  0ka80],  and,  moreover,  the  refractive  index 
control  of  crystalline  materials  is  also  difficult.  High 
doping  of  a refractive  index  modifier  to  the  host  crystal 
sometimes  makes  it  difficult  to  grow  a large  single  crystal. 
Furthermore,  as-grown  crystals  have  specific  crystal  surfaces, 
this  leads  to  roughness  on  the  surface  and  core/clad  interface 
of  the  waveguide  introducing  a large  scattering  loss.  The 
major  loss  increase  in  polycrystalline  materials  is  caused  by 
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light  scattering  at  the  grain  boundary.  The  scattering 
intensity  is  very  strong.  Therefore,  the  transmission  loss  of 
polycrystalline  waveguide  is  larger  than  expected  from  the 
intrinsic  material  property. 

Amorphous  materials  are  ideal  for  optical  waveguides, 
because  some  glasses  are  highly  transparent  in  visible  and 
infrared  wavelength  regions.  These  glasses  permit 
comparatively  easy  fabrication  of  fibers  without  introducing 
extrinsic  scattering  sources  such  as  specific  crystal  surface 
roughness  or  grain  boundary.  Moreover,  the  refractive  indices 
of  these  glasses  can  be  easily  controlled  by  changing  the 
constituent  composition,  which  does  not  severely  affect  the 
glass  formation  condition.  Rayleigh  scattering  loss, 
originating  from  the  density  and  composition  fluctuations  of 
amorphous  materials,  is  larger  than  that  of  crystalline 
materials.  However,  the  Rayleigh  scattering  loss  is  much 
smaller  than  the  loss  caused  by  grain  boundary  scattering  in 
polycrystalline  materials. 

1.1.2  Ge02-Si02  Waveguides 

Oxide  glasses,  especially  silicate  glass,  have  high 
transparency  in  the  visible  and  near-infrared  wavelength 
region.  Chemical  and  mechanical  durabilities  are  also  high. 
The  refractive  index  of  silicate  glass  can  be  changed  easily 
by  doping  with  other  oxides  such  as  Ge02,  A1203,  and  Ti02. 
Furthermore,  the  raw  materials  are  relatively  cheap. 
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Therefore  silica  glass  is  one  of  the  most  suitable  materials 
for  low  loss  optical  waveguides  and  has  been  widely  used  in 
various  optical  technologies  today. 

With  regard  to  the  possible  refractive  index  modifiers 
that  can  be  used  in  fabricating  low-loss  waveguides,  most 
transition  metal  oxides  in  oxide  glass  introduce  strong 
absorption  in  the  visible  and  near-infrared  wavelength  region. 
Therefore,  transition  metals  should  be  removed  from  the 
waveguide  materials,  in  order  to  make  low-loss  waveguides. 
Some  compounds  of  S,  Se,  Te,  and  P tend  to  form  colloids, 
which  may  cause  coloration  of  the  materials;  therefore,  they 
are  not  preferable  as  waveguide  material  either.  Alkali  metal 
oxides  are  sometimes  used  to  decrease  the  glass  softening 
temperature.  Chemical  durability  of  the  glass  containing 
alkali  metal  oxides  is  comparatively  poor,  because  alkali  ions 
in  the  glass  dissolve  in  water.  The  use  of  alkali-metal  and 
alkaline-earth  metal  and  heavy  metal  atoms  increases  the 
difficulty  of  reducing  the  possible  transition  metal 
contaminants  to  low  levels. 

The  kinds  of  oxides  suitable  for  making  low-loss 
waveguides  are  limited  strongly  by  the  reasons  mentioned 
above.  Ge02  glass  has  an  atomic  structure  similar  to  that  of 
silica  glass  and  has  a strong  absorption  band  at  11.4  jum 
[Bor 69]  due  to  the  v 3 vibration  mode.  The  absorption  peak  for 
Ge02  on  the  UV  side  is  around  185nm[ Yue82 ] ; an  absorption  band 
at  245nm  has  been  frequently  observed,  which  is  attributed  to 
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defect  structures  within  the  glass  [Jac85].  Addition  of 
germania  to  silica  doesn't  affect  optical  characteristics  too 
much  except  for  a slight  shift  of  absorption  peaks  to  longer 
wavelength.  Therefore,  germania-doped  silicate  glasses  are 
the  most  important  materials  used  in  optical  fiber  today. 

One  of  the  major  reasons  why  the  high  silica  glass 
systems  are  the  most  preferable  materials  for  optical  fiber  is 
that  high  purity  glass  can  be  easily  prepared  from  chloride 
raw  materials  through  vapor  reaction  and  deposition  processes. 
The  vapor  phase  process  has  an  important  feature:  the 
refractive  index  control  or  dopant  concentration  control  is 
comparatively  easier.  It  is  also  easy  to  make  a precise 
control  in  shape  and  size  of  fiber  preforms  obtained. 

Today,  Ge02  is  used  as  a dopant  to  increase  the 
refractive  index  of  the  silica  core  because  it  seldom  produces 
light  absorbing  ions  in  the  process  of  glass  fabrication.  The 
one  disadvantage  of  using  Ge02  is  that  its  vapor  pressure  at 
high  temperatures  is  too  high  to  control  the  concentration 
precisely.  When  Ge02-Si02  glass  is  heated  up  to  high 
temperatures,  Ge02  evaporates  from  the  surface.  Therefore,  it 
is  very  difficult  to  make  a transparent  glass  directly  from 
the  raw  material  vapors. 

This  difficulty  was  solved  by  the  development  of  a two 
step  process  or  soot  process.  The  soot  process  is  the  method 
of  making  a transparent  article,  which  includes  vaporizing 
hydrolyzable  compounds  of  silicon  and  doping  element  into  a 
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flame  of  combustible  gas  to  decompose  the  vapor  and  to  form 
finely  comminuted  doped  silica  and  then  vitrifying  the  silica 
to  a transparent  body  by  subsequent  heat  treatment.  There  are 
three  major  practical  soot  processes  for  making  preforms: 
Modified  Chemical  Vapor  Deposition (MCVD)  [Mac80] , Outside 
Vapor  Deposition (OVD)  [Sch80],  and  Vapor-phase  Axial 
Deposition (VAD)  [Iza79].  High  silica  fibers  are  made  by 
drawing  down  a preform,  which  has  a structure  and  composition 
similar  to  that  of  the  fiber  to  be  obtained. 

The  difference  between  these  three  processes  lies  in  the 
geometrical  growth  direction  of  the  fine  glass  particles  or 
soot.  Glass  particles  are  deposited  on  a substrate  tube  in 
the  radial  direction  in  case  of  MCVD  and  OVD  processes.  On 
the  other  hand,  glass  particles  are  deposited  on  a substrate 
rod  in  the  axial  direction  in  the  VAD  process.  Another 
difference  is  the  heat  source.  Although  usually  oxy-hydrogen 
gas  flames  are  used  for  the  chemical  reaction,  hydrocarbon  gas 
flames  such  as  CH4-02,  C3H8-02  or  RF  discharge  gas  plasma 
torches  can  be  substituted  for  oxy-hydrogen  gas  flames. 

The  sol-gel  method  has  been  one  of  the  promising  methods 
for  synthesizing  high  purity  silica  glass[Hen90,  Sak80]  which 
should  be  applicable  to  fiber  optics [Han80 ] . In  addition,  the 
method  has  an  advantage  in  saving  energy  over  conventional  CVD 
processes . 

It  has  been  shown  that  the  method  is  also  useful  for  the 
fiber  optics  application,  using  undoped  silica  glass  [Sus82, 
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Mat84 ] . However,  the  processing  of  doped  silica  has  not  been 
successful  because  of  poor  transmission  loss  properties 
[Puy82].  The  8%Ge02  doped  Si02  fiber  prepared  from  the  preform 
made  by  the  sol-gel  method  possesses  a loss  minimum  around  22 
dB/km  at  830nm.  The  high  overall  loss  was  believed  to  be  due 
to  transition  metal  impurities  from  the  germanium  ethoxide 
precursor  or  from  the  residual  water  in  the  fiber.  Both 
preform  rods  with  uniform  refractive  index  as  well  as  gradient 
refractive  index[Shi86]  have  been  reported.  Susa  et 
al. [ Sus90 ] also  discussed  the  foaming  problem  upon  heating 
glasses  at  high  temperatures  after  consolidation  of  gels, 
which  has  to  be  solved  before  fabricating  low-loss  fibers. 

1.2  Literature  Review 

1.2.1  The  Basic  Concepts  of  Waveguiding 

Dielectric  waveguides  are  the  structures  that  are  used  to 
confine  and  guide  light  in  guided-wave  devices  and  circuits  of 
integrated  optics.  Several  examples  of  waveguide  geometries 
are  shown  in  Figure  1.1.  A well-known  dielectric  waveguide  is 
the  optical  fiber  which  usually  has  a circular  cross-section 
as  shown  in  (a)  . In  contrast,  the  guides  of  interest  to 
integrated  optics  are  usually  planar  structures,  such  as 
planar  films  or  strips,  as  shown  in  (b)-(d).  The  simplest 
dielectric  waveguide  is  the  planar  waveguide  shown  in  Figure 
1.2,  where  a planar  film  of  refractive  index  nf  is  sandwiched 
between  a substrate  and  a cover  material  with  lower  refractive 
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Figure  1.1  Examples  of  guided  wave  geometries.  (a)  fiber 
waveguide,  (b)  planar  or  slab  waveguide,  and  (c)  (d)  channel 
waveguides . 
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cladding 


n 


nf  > ns,  nc 

An  «nf  - ns=  0.001  ~0.1 
nc  = 1 for  air 
d « 1//m 


Figure  1.2 


Schematics  of  a planar  waveguide. 
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indices  ns  and  nc  (nf>n5>nc)  . Often  the  cover  material  is  air, 
in  which  case  nc=l.  Typical  differences  between  the  indices 
of  the  film  and  the  substrate  range  from  0.001  to  0.1,  and  a 
typical  film  thickness  is  1 /urn.  The  light  is  confined  by 
total  internal  reflection  at  the  film-substrate  and  film-cover 
interfaces. 

Since  light  is  a form  of  electromagnetic  wave,  its 
properties  must  be  explained  in  terms  of  Maxwell's  equations. 
Therefore  the  rigorous  way  to  calculate  the  intensity  and 
phase  of  a light  wave  is  to  derive  the  wave  equation  from 
Maxwell's  equations  and  solve  directly.  However,  because  the 
wavelength  of  light  is  very  short,  it  is  often  convenient  for 
demonstration  purposes  to  use  the  concept  of  a "light  ray"  and 
geometrical  optics,  or  the  ray  equations,  which  are  obtained 
by  starting  from  Maxwell's  equations  and  then  taking  the 
approximation  X=>0  (where  X is  the  wavelength)  . With  the 
simplification  of  ray  approximation,  the  refraction  of  the 
light  at  an  interface  separating  two  lossless,  isotropic, 
homogeneous  dielectric  media  of  refractive  index  n,  and  n2,  as 
shown  in  Figure  1.3,  can  be  described  by  Snell's  law  and 
Fresnel  formulas,  respectively. 

When  a coherent  light  wave  is  incident  at  an  angle  0, 
between  the  wave  normal  and  the  normal  to  the  interface.  In 
general,  the  wave,  having  a complex  amplitude  A,  is  partially 
reflected  and  refracted.  The  exit  angle  02  of  the  refracted 
wave  C is  given  by  Snell's  law 
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refracted  ray 


incoming  ray  reflected  ray 


Figure  1.3  Ray  picture  of  light  at  the  interface  of  two 

media  of  different  refractive  indices. 
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(1.2.1) 


j^sinfl,  = n2sin02 


The  reflected  wave  has  a complex  amplitude  B at  the  interface 
linearly  related  to  A by  a complex  reflection  coefficient  R 


The  reflection  coefficient  depends  on  the  angle  of  incidence 
and  the  polarization  of  the  light,  and  is  given  by  the  Fresnel 
formulas  [Kni76].  For  TE  polarization(i.e.  electric  fields 
perpendicular  to  the  plane  of  incidence  spanned  by  the  wave 
normal  and  the  normal  to  the  interface) 


(1.2.2) 


B=R  A 


n,c os0j  -n2 cos02 


njcosflj  + n2cos02 


(1.2.3) 


The  corresponding  formula  for  TM  polarization  (with  the 
magnetic  fields  perpendicular  to  the  plane  of  incidence)  is 
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n2cos0,  -nlcosd2 
n2cosdl  +nl  cos02 

n2c os0,  -nx\jnl  -n^sin^, 

n2cos0j  +nxyn%  -n12sin201 


(1.2.4) 


At  the  critical  angle  of  incidence,  0C/  transmission  of  the 
light  wave  into  medium  2 ceases  completely,  and  all  the  energy 
is  totally  reflected  at  the  boundary;  where 


(1.2.5) 


When  the  input  angle  exceeds  the  critical  angle,  no  light  rays 
propagate  into  the  medium  2;  that  is,  there  is  no  refracted 
ray.  Region  2 then  contains  only  an  evanescent  component 
which  is  guided  along  the  boundary  surface  in  the  z direction. 
The  totally  reflected  light  has  value  |r|=1,  R is  complex  and 
a phase  shift,  <p , is  imposed  on  the  reflected  light. 


R = eW 


(1.2.6) 


According  to  the  Fresnel  formulas  the  phase  shifts  0^  and  0TM 
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corresponding  to  the  two  polarization  states  can  be  expressed 
as [Kog79 ] 


This  phase  shift  is  dependent  on  the  angle  of  incidence  and 
the  refractive  index  of  both  media,  n,  and  n2. 

Consider  a three-layer  system  as  in  asymmetric  planar 
waveguide  structure,  as  shown  in  Figure  1.2.  In  general,  it 
requires  nf>ns>nc  and  two  critical  angles  of  interest,  0,  for 
total  reflection  from  the  film-substrate  interface  and  6C<6S 
for  total  reflection  from  the  film-cover  interface.  When  the 
angle  of  incidence  6 is  increased,  three  distinct  cases  will 
happen  as  shown  in  Figure  1.4[Kog79]. 

For  small  angles  0<0S,  0C,  light  incident  from  the 
substrate  side  is  refracted  according  to  Snell's  law  and 
escapes  from  the  guide  through  the  cover,  see  (a) . There  is 
no  confinement  of  light,  and  the  electromagnetic  mode 
corresponding  to  this  picture  is  called  a "radiation  mode". 

When  6 is  increased  such  that  6t>6>6cl  the  situation  (b) 


(1.2.7) 


(1.2.8) 
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Figure  1.4  Zig-zag  wave  pictures  for  the  radiation 

inodes  (a),  the  substrate  inodes  (b),  and  the  guided  inodes  (c)  in 
a slab  waveguide [Kog79 ] . 
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can  be  observed.  The  light  incident  from  the  substrate  is 
refracted  at  the  film-substrate  interface,  totally  reflected 
at  the  film-cover  interface,  refracted  back  into  the  substrate 
through  which  the  light  escapes  from  the  structure.  In  this 
case,  there  is  no  confinement  either,  and  is  called  a 
•'substrate  radiation  mode”. 

Finally,  when  9 is  large  enough  that  0S,  9C<9 , the  case 
(c)  will  occur,  i.e.  total  internal  reflection  at  both 
interfaces.  The  light,  once  it  is  inside,  is  trapped  and 
confined  in  the  film  and  propagates  in  a zig-zag  path.  This 
case  corresponds  to  a ”guided  mode”  of  propagation. 

An  optical  wave  guided  successfully  by  total  internal 
reflections  may  be  represented  by  bundles  of  rays  called 
"modes”.  By  using  the  idea  of  a mode,  the  phenomenon  of  light 
guidance  can  also  be  treated  as  a wave  motion.  The  formation 
of  a mode  in  a dielectric  waveguide  can  be  shown  schema- 
tically as  in  Figure  1.5. 

The  physical  picture  of  guided  light  propagation  is  that 
of  light  traveling  in  zig-zag  fashion  through  the  film;  each 
ray  represents  an  uniform  plane  wave  from  the  wave  point  of 
view.  More  precisely,  it  is  a picture  of  two  superimposed 
uniform  plane  waves  with  wavenormals  which  follow  the  zig-zag 
path  indicated  in  the  figure  and  which  are  totally  reflected 
at  the  film  boundaries.  These  waves  are  monochromatic  and 
coherent  with  angular  frequency  w,  free-space  wavelength  X, 
and  they  travel  with  a wavevector  k in  the  direction  of  the 
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Figure  1.5 
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Ray  diagram  of  a guided  mode  in  the  planar 
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wave  normal,  This  wave  vector  k of  the  wave  within  the  film 
can  be  expressed  as 


k = k„  nf 


(1.2.9) 


where  the  absolute  value  of  kQ  is 


. = 2n  _ o>  (1.2.10) 

0 ~T  c 

where  kc  is  the  propagation  (or  phase)  constant  of  the  plane 
wave  in  a vacuum,  w is  the  angular  velocity  of  light,  and  c 
the  velocity  of  light  in  vacuum. 

This  plane  wave  can  be  decomposed  into  two  component 
plane  waves,  propagating  in  the  axial  (or  z)  and  transverse 
(or  x)  directions,  as  shown  in  Figure  1-5.  the  component  of 
the  propagation  constant  in  the  axial  direction  is  j3,  that  is 


j3  = — =ka  nf  sin0 

\r  J 

P 


(1.2.11) 


On  the  other  hand,  the  transverse  component  of  the  propagation 
constant  y can  be  given  by 
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(1.2.12) 

y = ±k0nf  cosd 

As  the  ray  travels  through  path  A-B-C  in  Figure  1.5,  it 
suffers  phase  shifts  at  both  interfaces  due  to  total  internal 
reflection  as  well  as  a phase  change  while  propagating  inside 
the  film.  If  the  total  accumulated  phase  change  is  an 
integral  multiple  of  2n , the  phase  front  after  total  internal 
reflection  at  A as  it  advances  to  C comes  in  phase  with  the 
phase  front  which  is  reflected  at  C,  resulting  in  a 
constructive  interference.  Under  these  conditions,  the  film 
supports  a ray  corresponding  to  a specific  incident  angle  d, 
a standing  wave  is  established  in  the  transverse  direction 
therefore,  the  field  distribution  in  the  transverse  direction 
does  not  change  as  the  wave  propagates  in  the  axial  direction. 
This  transverse  resonance  condition  can  be  represented  as 
[Tie71 , Kog79] 


(1.2.13) 

2ka  d nf  cosd  - 2 <ps  - 2 <pc  = 2mn 


where  d is  the  film  thickness,  m is  an  integer  and  the  phase 
shifts  <p,  and  <pc  for  TE  mode(i.e  the  electric  field 
polarization  perpendicular  to  the  plane  of  incidence)  can  be 
expressed  in  terms  of  film  parameters 
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(1.2.14) 


(1.2.15) 


while  for  TM  mode(i.e.  the  electric  field  polarization 
parallel  to  the  plane  of  incidence) 


For  a given  set  of  film  parameters,  the  ray  propagates  at 
specific  discrete  angles  (or  sometimes  none) . 

1.2.2  Review  of  the  Progress  of  Integrated  Optics 

The  guiding  of  light  beams  along  dielectric  layers, 
first  realized  in  the  early  sixties  [Mil69],  has  stimulated 
the  growth  of  a new  class  of  passive  and  active  components 


(1.2.16) 


and 


(1.2.17) 
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using  light  guided  by  films.  By  utilizing  various  masking 
and/or  lithography  techniques,  a variety  of  integrated  optical 
components  capable  of  performing  a wide  range  of  operations  on 
light  waves  have  been  realized.  Because  of  their  small 
dimensions,  light  weight  and  low-power  requirements,  it  was 
then  projected  that  such  optical  components  could  replace 
electrical  circuitry  in  integrated  electronic  equipment.  In 
addition,  the  optical  elements  would  provide  the  advantages  of 
greater  bandwidth  and  immunity  from  electromagnetic 
interference.  These  advantages  provided  by  integrated  optics 
can  benefit  areas  such  as  telecommunications,  spectrum 
analyzer,  gyroscopes,  digital  correlators,  optical  switches 
multiplexers,  and  signal  processors  etc.  [Hal86,  Car87,  Hut88, 
Alf86,  Gar86].  However,  the  major  disadvantage  of  the  optical 
integrated  circuits  is  the  high  cost  of  developing  new 
fabrication  technology. 

The  research  can  be  generally  divided  into  two 
categories: 

(1) .  Development  of  waveguide  materials  and  fabrication 

processes . 

(2)  . Development  of  new  devices  or  new  configuration  for 

the  existing  devices. 

A variety  of  waveguide  materials  such  as  polymers,  fused 
silica,  silicate  glasses,  lithium  niobate,  lithium  tantalate, 
zinc  oxide,  III-V  semiconductors  etc.,  have  been  used  to 
fabricate  waveguides.  The  associated  fabrication  processes 
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includes  ion  diffusion,  ion-exchange,  ion-implantation  and 
various  thin  film  deposition  methods,  such  as  plasma 
polymerization,  sputtering,  dipping (or  spinning) , epitaxial 
layer  growth  etc. 

Kaminow  and  Carruthers [Kam73 ] have  used  out-diffusion 
of  Li  ions  to  form  waveguiding  layers  in  LiNbOj  and  LiTa03 
crystals.  In-diffusion  of  various  transition  metal  elements 
has  been  used  by  Schmidt  and  Kaminow [Sch7 4 ] to  prepare  low- 
loss  planar  waveguides  in  LiNb03.  Low-temperature  proton- 
exchange  process  has  been  used  to  fabricate  stable  low-loss 
LiNb03  waveguides  by  Suchoski  et  al.[Suc88].  Various 
researchers  have  fabricated  glass  waveguides  with  different 
ion-exchange  processes  on  silicate  glasses [And73 , Iza72, 
Poy90,  Ram88],  Schneller  et  al.  have  applied  a proton- 
implantation  process  to  create  waveguiding  layers  in  fused 
silica [ Sch68 ] . Several  polymer  materials  have  been  processed 
into  different  waveguide  structures  by  solution-deposited 
method [Ulr72 ] , plasma-polymerization  method [Tie72 , Vas72, 
Tie73],  and  Photo-induced  polymerization  process [Die90 , 
Boo89].  Goell  and  Standley  have  reported  using  an  RF 
sputtering  technique  to  form  a waveguiding  layer  of  Corning 
7059  glass  on  a microscope  slide[Goe69] . 

In  all  these  cases,  the  propagation  loss  of  the  waveguide 
and  ease  of  fabrication  have  been  the  major  criteria  for 
accepting  or  rejecting  a waveguide  material.  To  date,  Ti- 
indiffused  LiNb03  waveguides  and  the  various  silicate  glass 


24 


waveguides  are  two  of  the  most  important  systems  under 
studied. 

1.2.3  Sol-Gel  Processes 
1 . 2 . 3 . 1 Background 

The  history  of  sol-gel  processing  begins  in  the  1800's 
with  the  studies  by  Ebelman  [Ebe46,  Ebe47]  and  Graham 
[Grah64].  These  early  investigations  observed  that  the 
hydrolysis  of  tetraethylorthosilicate,  Si(OC2H5)4,  at  low  pH 
resulted  in  the  formation  of  Si02.  It  was  recognized  by 
Geffcken  [Gef51]  in  the  early  1930's  that  alkoxides  could  be 
used  in  the  preparation  of  oxide  films.  This  process  was  soon 
developed  by  the  Schott  glass  company  in  Germany  as  reviewed 
by  Schroeder  [Sch69],  Around  the  same  time,  mineralogists 
became  interested  in  the  use  of  sols  and  gels  for  the 
preparation  of  homogeneous  powders  in  studies  of  phase 
equilibria  [Ewe35,  Barr50] . This  method  was  popularized  in 
the  ceramic  community  by  Roy  [Roy56,  Roy69,  McC71]  for  the 
preparation  of  homogeneous  powders.  They  used  this  technique 
and  synthesized  a number  of  oxide  compositions  that  cannot  be 
produced  by  the  traditional  powder  methods.  The  concurrent 
work  of  Her  [Ile55]  and  Stober  [Sto68]  revealed  the 
importance  of  catalysts,  alkoxide  types  and  concentration, 
solution  pH  and  reaction  temperature  as  processing  parameters. 

During  the  past  decade,  an  enormous  growth  of  literature 
has  been  associated  with  the  development  of  the  sol  gel 
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process.  This  was  stimulated  largely  by  the  works  of 
Phalippou  and  Zarzycki [Pha82 ] as  well  as  Yoldas  [Yol75].  They 
demonstrated  the  feasibility  of  preparing  large  , dense, 
monolithic  objects  of  silica  and  alumina,  respectively.  Sol- 
gel  processing  is  a method  for  producing  glass  or  ceramics 
from  solutions  of  metal-organic  or  inorganic  compounds  through 
a series  of  chemical  reactions  (e.g.  hydrolysis  and 
condensation)  which  cause  the  formation  of  the  polymeric 
species  (or  colloids)  dispersed  in  solution (or  called  sol) . 
This  sol  then  gels  into  a three-dimensional  network.  The 
final  products  could  be  obtained  after  drying  and  densifying 
the  gel.  At  the  functional  group  level,  three  general 
reactions  can  be  used  to  describe  the  sol-gel  process: 

(1) .  Hydrolysis  versus  esterification 

sM-OR  + H20  <===>  =M-OH  + ROH 

(2)  . Alcohol  condensation  versus  alcoholysis 

sM-OR  + HO -Ms  <===>  =m— 0-Ms  + ROH 

(3) .  Water  condensation  versus  hydrolysis 

sM-OH  + HO -Ms  <===>  =M-0— M=  + H20 

where  R is  an  alkyl  group,  CxH2x+1;  M is  metal  such  as  Si  or  Ge. 

Such  a technology  offers  enormous  opportunity  for 
preparing  crystalline  ceramics  and  glasses  with  high  purity, 
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chemical  homogeneity,  and  novel  composition  and  microstructure 
at  relatively  low  temperatures.  Applications  for  sol-gel 
materials  cover  optical,  electronic,  structural  ceramics,  and 
recently  discovered  ceramic  superconductors,  in  the  form  of 
fibers,  coatings,  powders,  and  shaped  monolithic  parts  [Uhl84, 
Bri88a] . 


1.2. 3. 2 Film  formation 

One  of  the  most  technologically  important  aspects  of  sol- 
gel  processing  is  that,  prior  to  gelation,  the  fluid  sol  or 
solution  is  ideal  for  preparing  thin  films  by  dipping, 
spinning,  or  spraying.  Compared  to  conventional  thin  film 
forming  processes,  such  as  CVD,  evaporation,  or  sputtering, 
sol-gel  film  formation  requires  considerably  less  equipment 
and  is  potentially  less  expensive.  However,  the  most 
important  advantage  of  this  process  is  the  ability  to  control 
precisely  the  composition  and  the  microstructure  of  the 
deposited  film,  i.e.  the  pore  volume,  pore  size,  and  surface 
area. 

Scriven[Bet88]  described  the  dip-coating  process  into 
five  stages:  immersion,  start-up,  deposition,  drainage,  and 
evaporation.  Figure  1.6[Bri91]  is  a schematic  representation 
of  the  steady  state  deposition  stage  of  the  dip-coating 
process.  Dilute,  non-interacting,  polymer  species  that  make 
up  the  coating  bath  are  concentrated  on  the  substrate  surface 
by  gravitational  draining  accompanied  by  vigorous  evaporation 
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Figure  1.6  Schematics  of  the  sol-gel  film  deposition 
during  the  dip-coating  process [Bri91 ] . 
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and  further  condensation  reactions.  Correspondingly,  the 
viscosity  progressively  increases  due  to  the  increasing 
concentration  and  further  condensation  reactions  promoted  by 
the  increasing  concentration.  This  polymer  growth  process  is 
believed  to  proceed  similarly  to  a cluster-cluster  aggregation 
model  proposed  by  Meakin[Mea83 ] . Gelation  occurs  at  the 
moment  when  the  condensing  network  is  sufficiently  stiff  to 
withstand  flow  due  to  gravity  yet  is  still  filled  with 
solvent. Several  factors  distinguish  structural  evolution  in 
films  from  that  of  bulk  gels [Bri88b] . 

(1)  . Overlap  of  the  deposition  and  evaporation  stages 
establishes  a competition  between  evaporation  (which  compacts 
the  structure)  and  continuing  condensation  reactions  (which 
stiffen  the  structure) . In  bulk  systems,  the  gelation  and 
drying  stages  are  normally  separated. 

(2)  . Compared  to  bulk  systems,  aggregation,  gelation,  and 
drying  occur  in  seconds  to  minutes  during  dipping  or  spinning 
rather  than  days  or  weeks. 

(3)  . The  short  duration  of  the  deposition  and  drying 
stages  causes  films  to  experience  considerably  less  aging 
(crosslinking)  than  bulk  gel.  This  generally  results  in  a 
more  compact  dried  structure. 
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(4) . Fluid  flow,  due  to  draining,  evaporation,  or  spin- 
off, combined  with  attachment  of  the  precursor  species  to  the 
substrate,  impose  a shear  stress  within  the  film  during 
deposition.  After  gelation,  continued  shrinkage  due  to  drying 
and  further  condensation  reactions  creates  a tensile  stress 
within  the  film.  Bulk  gels  are  not  constrained  in  any 
dimension. 

1.2. 3. 3 Attempts  at  making  GeO?-SiO-,  gel  glass 

As  mentioned  earlier,  a lot  of  research  has  been 
conducted  on  making  Ge02-Si02  glass  preforms  for  fiber  drawing 
in  the  past  few  years.  Puyane  et  al.[Puy82]  first  prepared 
silica-germania  gels,  containing  up  to  30wt%  Ge02.  Sols  were 
made  up  of  10%  TEOS,  43%  ethanol,  4%  methanol,  33%  H20,  and  a 
variable  percentage  of  germanium  ethoxide  depending  on  the 
final  oxide  composition  required. 

Mukherjee[Muk84 ] prepared  Ge02-Si02  gels  up  to  80wt%  Ge02 
with  two  gel  preparation  approaches: 

(1)  . The  starting  compounds  were  TEOS  (or  TMOS)  and 
germanium  ethoxide.  The  general  procedure  consisted  of  the 
following  steps:  (a),  mix  alkoxysilane  with  three  times  its 
volume  of  anhydrous  ethanol  at  4 0°C.  (b)  . partial  hydrolysis 
at  40°C  with  1 mole  water  per  mole  of  silicon  alkoxide, 
acidified  with  HC1  (0,003  mole  per  mole  Si(OR)4),  stir  for  two 
hours  at  40°C,  pH=2 . (c)  . cool  to  room  temperature  and  add 
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dilute  ethanolic  solution  of  germanium  ethoxide,  stir  for  0.5 
hour  to  1 hour  at  room  temperature.  (d)  . add  4 moles  of 
H20(diluted  with  ethanol)  per  mole  of  M(OR)4.  (e)  . stir  at 
room  temperature  for  0.5  to  1 hour  then  cast  for  gel 
formation. 

(2)  . The  starting  compounds  were  TMOS  and  germanium 
ethoxide.  The  general  procedure  consisted  of  the  following 
steps:  (a) . mix  TMOS  with  methanol  at  room  temperature.  (b) . 
add  aqueous  NH4OH  solution  (0.0002  mole  NH40H  per  mole  TMOS  and 
1 mole  H20  per  mole  TMOS).  (c) . stir  for  0.75  hour  at  room 
temperature,  pH=8 . (d)  . add  ethanolic  solution  of  Ge(OC2H5)4 
and  stir  for  0.25  to  0.5  hour.  (e) . add  water  diluted  with 
ethanol,  4 moles  per  mole  M(OR)4,  mixed  with  0,0001  mole  NH4OH 
per  mole  alkoxide.  (f) . stir  for  0.25  to  0.5  hour  at  pH=7. 
(g) . cast  for  gelation. 

Consolidated  samples  could  not  be  produced  with  these 
methods.  Crystallization  of  a-cristobalite  phase  in  high 
silica  gel  and  Ge02  hexagonal  phase  in  gels  containing  Ge02 
larger  than  56wt%  were  observed. 

Preparation  of  graded-index (GRIN)  Ge02-Si02  glasses  by 
sol-gel  method  was  reported  by  Shingyouchi  et  al.[Shin86]. 
The  wet  gel  was  prepared  by  using  TMOS,  germanium  ethoxide, 
ethanol,  water,  and  use  HC1  as  catalyst.  The  gel  was  then 
leached  in  a water  bath,  followed  by  washing,  drying,  and  the 
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subsequent  heat-treatment. 

Shibata  et  al.[Shi86]  fabricated  Ge02-Si02  core  optical 
fibers  by  using  TEOS  and  germanium  butoxide  as  starting 
compounds.  The  starting  compounds  were  diluted  with  ethanol 
of  commercial  grade,  and  ammonium  hydroxide  was  added  as 
catalyst.  Pure  silica  cladding  was  formed  on  the  Ge02  doped 
Si02  core  glass  by  a conventional  VAD  process.  The  minimum 
optical  loss  of  the  fiber  was  measured  as  9 dB/km  at  1.07/im, 
which  was  attributed  to  the  remaining  OH  impurity  (>200ppm) . 

1.2.4  Low  Pressure  Chemical  Vapor  Deposition  Process (LPCVD) 
1 . 2 . 4 . 1 Background 

The  growth  of  thin  silicon  dioxide  as  well  as  doped 
silica  layers  by  chemical  vapor  deposition  has  been  one  of  the 
most  important  methods  for  semiconductor  microelectronic 
device  processing  and  has  been  investigated  intensively 
[Ben87,  Pav86,  Tob80,  wah75,  Ker77] . These  oxide  layers  are 
used  mainly  for  interlayer  dielectrics,  protective,  and 
passivation  layers. 

The  reasons  for  the  widespread  acceptance  of  LPCVD  for 
producing  thin  dielectric  film  are  reduced  processing  cost  and 
superior  product  quality  that  can  be  achieved,  i.e.  good 
uniformity,  low  particulate  inclusions  and  low  pin  hole 
density  [Mat80]  . 

Two  basic  parameters  govern  the  deposition  rate  and 
uniformity  of  films  grown  by  all  CVD  processes: 
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(1)  rate  of  mass  transfer  of  reactant  gases  to  the 
substrate  surface. 

(2)  rate  of  surface  reaction  of  the  reactant  gases  at  the 
substrate  surface. 

The  mass  transfer  of  gases  involves  their  diffusion 
across  the  slowly  moving  boundary  adjacent  to  the  substrate 
surface.  The  thinner  this  boundary  layer  and  the  higher  the 
diffusion  rate,  the  greater  the  mass  transport  that  results 
[Ros77 ] . Surface  reaction  rates,  on  the  other  hand,  depend 
mainly  upon  reactant  concentration  and  deposition  temperature. 
By  lowering  the  pressure  in  the  reactor,  the  diffusivity  for 
LPCVD  is  several  orders  of  magnitude  («1000)  higher  than  that 
of  an  atmospheric  reactor.  The  LPCVD  boundary  layer  is 
thicker  by  a factor  of  3 to  10  which,  however,  does  not  offset 
the  huge  diffusivity  difference.  Thus  low-pressure  deposition 
conditions  enhance  mass  transfer  greatly,  making  close-spaced 
substrate  positioning  a practical  and  highly  economical 
technique [Ker79 ] . Moreover,  low-pressure  deposition  minimizes 
the  formation  of  homogeneously  gas  nucleated  particles,  i.e. 
soot  formation  is  depressed,  and  surface  reactions (or 
heterogeneous  reactions)  become  dominant  which  make  high 
quality  thin  film  deposition  possible.  Vertically-positioned 
substrates  may  eliminate  codeposition  of  fine  particles  in  the 
reactor  into  the  film. 

The  key  limitation  of  this  process  is  its  inherently  low 
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deposition  rate  due  to  the  fact  that  the  operating  total 
pressure  is  low  («0.001  atm).  However,  higher  deposition 
rates  are  attainable  by  using  more  concentrated  reactive 
gases. 

Several  different  reactant  gases  have  been  used  to  form 
silica  in  various  CVD  processes.  The  deposition  temperature 
in  the  reactor  determines  the  extent  of  oxide  conversion 
efficiency,  and  sometimes  the  kinds  of  the  reaction  products 
formed.  Some  reactions  are  described  as  follows: 

(1) .  Halide  process(or  soot  process) 

SiCl4  + 02  <===>  Si02  + 2 Cl2  T«1 000-13 00°C 

(2) .  TEOS  process 

Si  ( OC^Hj),*  + 02  <===>  Si02  + 8C02  + 10H2O  T«700°C 

(3)  . Hydride  process 

SiH4  + 02  <===>  Si02  + 2H2  T«300-500°C 

SiH4  + 202  <===>  Si02  + 2H20  T~600-1000°C 

Among  the  reactions  described  above,  the  hydride  process 
is  a low  temperature  process,  and  the  high  concentration  of 
bonded  hydrogen  can  be  eliminated  in  the  oxides  if  deposition 
of  the  material  is  conducted  in  the  appropriate  temperature 
range.  H-bonds  in  the  Si02  films,  existing  in  the  form  of  Si- 
OH  and  Si-H,  can  produce  recombination  centers  and  cause  deep 
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trapping [Kim88]  for  electronic  applications,  and  cause 
undesirable  optical  absorption  in  the  near-infrared  range. 

1.2. 4. 2 Deposition  of  GeCX-SiO?  Film  bv  LPCVD  Method 

Conventionally,  most  of  the  Ge02-Si02  binary  systems  were 
prepared  by  either  a melting  process  or  various  soot  processes 
for  fiber  preform  applications. 

There  are  very  few  studies  of  the  preparation  of  Ge02-Si02 
glass  by  the  LPCVD  process.  Huffman  et  al.[Huf86]  attempted 
to  prepare  Ge02-Si02  binary  amorphous  films  by  LPCVD  method, 
using  SiH4,  GeH4,  and  02  as  reactants.  The  02/ (Si02+Ge02)  ratio 
was  1:1,  and  a reactor  pressure  of  «400  mtorr  and  deposition 
temperature  of  430°C  was  used.  The  film  was  made  to  a 
thickness  of  about  3/m  on  stainless  steel  substrates.  The 
film  contains  3.6  mole%  Ge02.  Owing  to  the  low  deposition 
rates,  films  containing  higher  Ge02  content  could  not  be  made. 
The  film  quality  was  poor  according  to  the  reported  SEM 
micrographs. 

1.2.5  Photosensitive  Response  in  Ge02-Si02  Waveguides 

It  was  first  discovered  by  Hill  et.al. [Hil78]  that 
certain  germania-doped  glass  optical  fibers  can  show  a 
photosensitive  response  to  a single  longitudinal  mode  Argon 
laser  beam,  leading  to  the  formation  of  light-induced  Bragg- 
matched  phase  grating.  This  phenomenon  has  been  observed  in 
both  germanium  silicate  fibers [Hil78 ] and  pure  germania  planar 
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waveguides [Yin83 ] . Despite  a recent  resurgence  of  interest  in 
this  phenomenon,  there  are  still  many  questions  surrounding 
the  growth  of  these  gratings.  Two  types  of  models  have  been 
proposed  to  explain  this  phenomenon. 

An  early  experiment  provided  evidence  that  the  process 
begins  with  two  photon  absorption  (TPA)  of  the  Argon  laser  beam 
[Lam81] . Later  it  was  suggested  by  Meltz  et.al. [Mel89]  that 
this  TPA  results  in  the  bleaching  of  an  ultraviolet  defect 
absorption  band  (-240  nm)  of  the  glass,  resulting  in  a 
modification  of  the  index  of  refraction  in  the  visible  region 
of  the  spectrum.  It  was  experimentally  demonstrated  by  Simmons 
et  al.  [Sim90]  that  a similar  wavelength  dependent  response 
for  both  the  generation  of  Ge  E'  centers  and  the  formation  of 
photoresponsive  gratings  in  germania-doped  silica  fibers. 
This  supports  the  model  [Mel89]  for  the  grating  formation 
whereby  TPA  of  laser  radiation  into  the  UV  bleachable 
absorption  band  (-240  nm)  of  an  oxygen-deficient  germania  site 
creates  hole  traps,  modeled  as  Ge  E'  center,  which  can  be 
detected  by  Electron  Spin  Resonance (ESR) . Recent  work  by  Hand 
and  Russell  [Han90]  suggested  that  new  absorption 
bands (Ge (1) :281nm,  and  Ge(2) :213nm)  are  formed  which 
contribute  to  the  index  change,  perhaps  even  dominating  over 
the  simple  effect  of  bleaching  the  240nm  band.  These  models 
are  referred  to  as  local  bleaching  models,  in  which  the 
modification  of  the  index  of  refraction  of  the  glass  is  a 
function  only  of  the  local  intensity  of  the  light.  This  type 
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of  model  was  demonstrated  by  Mizrahi  et . al . [Miz90]  to  suit 
experimental  results  better. 

The  other  type  of  model  is  a photorefractive-like  model, 
as  was  suggested  recently  by  Payne [ Pay89 ] , which  involves 
carrier  transport  over  macroscopic  distances.  As  Mizrahi 
et. al. [Miz90]  pointed  out  such  a model  predicts  a modulation 
depth  of  the  grating (in  steady  state)  independent  of  the 
intensity  of  the  laser  that  formed  it.  This  is  in  sharp 
contrast  with  experiment [Lam81] . 

In  the  local  bleaching  model [Miz90b] , a change  in  the 
real  part  of  the  dielectric  constant,  Ae(z,t),  in  the  visible 
region  of  the  spectrum  with  respect  to  the  local  intensity  of 
the  light  in  the  waveguide,  I(z,t),  can  be  written  as 


dAe (z, t) 
Jt 


= AI2{z,  t) 


(1.2.18) 


where  A is  an  material-related  factor,  which  is  related  to  the 
two  photon  absorption (TPA)  cross-section  of  the  245nm  band,  as 
well  as  other  material  parameters  which  are  at  present 
unknown . 

The  absorption  band  at  about  245nm  is  often  attributed  to 
an  oxygen  vacancy  defect  of  the  Ge02[Gar58b,  Yue82].  However, 
the  precise  microscopic  origin  of  this  band  is  still  the 
subject  of  discussion [ Sku84 , Rau87,  Tsa87].  Moreover,  silica 
glass  itself  may  also  contribute  two  absorption  bands  near 
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248nm(~5 . OeV)  , namely  247nm(5. 02eV)  for  B2a  band  and 
240.7nm(5. 15eV)  for  B2)3  band.  Both  of  them  were  observed  in 
oxygen-deficient  samples  and  can  be  distinguished  with 
photoluminescence  measurement  as  well  as  optical  absorption 
spectra  obtained  from  samples  heat-treated  under  different 
atmospheres [Toh89 ] . This  has  complicated  the  analysis  of  the 
245nm  absorption  band  for  Ge02-Si02  glass  system. 

1.2.6  Frequency  Doubling  in  Glass  Waveguides 

The  anomalous  photoinduced  second  harmonic  generation  in 
germania-doped  silica  fiber  was  first  observed  by  Osterberg 
and  Margulis [Ost86] . A similar  phenomenon  was  also  later  seen 
in  pure  Ge02  waveguides  [Kas89 ] . This  phenomenon  is  still 
largely  unexplained  today. 

Second  harmonic  generation (SHG)  is  one  of  the  best  known 
of  the  nonlinear  optical  effects.  First  observed  in  1961  from 
a quartz  crystal [Fra61] , it  is  the  basis  of  the  most  important 
nonlinear  optical  devices.  But  the  observation  of  efficient 
second  harmonic  generation  from  "glass"  was  totally 
unexpected.  Very  general  symmetry  arguments  prove  that  a 
material  with  a center  of  inversion  symmetry  cannot  exhibit 
dipole-allowed  second  harmonic  generation.  A number  of  models 
have  been  proposed  to  explain  this  effect.  Some  models 
suggest  that  dipoles  already  present  in  the  glass  are  somehow 
reoriented  by  interaction  between  infrared  radiation  and  the 
glass  to  break  its  average  inversion  symmetry [Sto8 7 , Les90]. 
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Recent  experimental  evidence[Miz90a,  Hib90a,  Hib90b]  of 
temperature  dependence  study,  polarization  study  of  SHG,  and 
infrared  erasure  of  self-organized  xa)  gratings  in  germania- 
doped  silica  fiber  argues  against  dipole  orientation  models, 
and  favors  carrier  transport  models.  Several  carrier 
transport  models  have  been  proposed  by  various  research 
groups [And89 , Que88,  Dia89,  Luc90] . These  authors  suggest 
that  efficient  SHG  from  an  initially  isotropic  glass  fiber 
implies  a laser-induced  modification  in  the  glass  structure, 
which  involve  some  form  of  carrier  transport,  perhaps  by  a 
hopping  process  or  by  diffusion,  leading  to  an  internal  d.c. 
electric  field  which  acts  on  the  x0)  of  the  glass  and  results 
in  an  effective  xa>*  The  carriers  might  come  from  charges 
released  from  certain  defect  sites  in  the  glass,  and  trapped 
at  other  sites  due  to  the  interaction  between  laser  and 
waveguide.  Tsai  et.al. [Tsa89]  have  observed  an  apparent 
increase  of  the  concentration  of  the  Ge  E'  defect  center  in 
fibers  which  are  conditioned  for  SHG. 

The  "defects"  in  the  glass  have  been  thought  to  play  an 
important  role  in  the  process.  Defect  (as  distinct  from 
impurity)  densities  depend  on  such  parameters  as  raw 
materials,  processing  method,  processing  temperatures,  fiber 
drawing  conditions,  and  atmosphere,  which  can  vary  between 
nominally  similar  fibers  made  by  different  manufacturers. 
This  in  part  explains  the  problem  of  the  conflicting  or 
confusing  experimental  results  found  in  the  literature  to  be 
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due  to  the  wide  variety  of  optical  fibers  actually  being 
studied  in  different  laboratories  around  the  world.  Another 
source  of  the  complication  is  that  it  would  be  desirable  to 
conduct  this  study  in  waveguides  that  are  always  single  mode 
at  both  the  fundamental  and  the  second-harmonic  wavelengths, 
but  standard  telecommunication  fibers  are  never  single  mode  in 
the  visible  region. 


CHAPTER  II 


RESEARCH  OBJECTIVES 
2 . 1 Statement  of  Problem 

As  mentioned  earlier,  the  fabrication  of  high  quality 
thin  films  and  devices  has  received  much  attention  in  recent 
years.  The  increasing  penetration  of  single  mode  fibers  into 
local-area  networks (LAN)  and  customer  premises  networks (CPN) 
will  inevitably  increase  the  requirement  for  low  cost  passive 
components,  such  as  branching  devices,  mixers,  etc.. 
Integrated  optical  devices  are  potentially  ideal  for  these 
applications,  provided  that  they  can  be  made  compatible  with 
single-mode  fibers.  One  of  the  factors  which  inhibits  the 
widespread  use  of  integrated  optical (10)  devices  is  the 
incompatibility  of  the  materials  used  (e.g.  LiNb03,  GaAs)  with 
silica-based  optical  fibers.  This  leads  to  problems  of 
Fresnel  reflection  and  mode  mismatch  at  the  interface,  and 
generally  makes  it  very  difficult  to  achieve  low-loss 
permanent  coupling  between  10  devices  and  optical  fibers.  The 
use  of  germania  as  the  core  dopant  and  silica  as  the  substrate 
ensures  that  these  waveguides  will  have  virtually  identical 
characteristics  (optical,  mechanical,  thermal,  etc.)  to 
single-mode  fibers.  Moreover,  glasses  in  the  form  of 
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waveguides,  both  fiber  and  planar,  have  recently  been  used  to 
study  various  nonlinear  optical  phenomena,  such  as  all-optical 
switching[Blo89]  and  spatial  soliton  propagation[Ait90] . 
Interesting  phenomena  such  as  photosensitive  grating  effect 
and  anomalous  photoinduced  second  harmonic  generation  (SHG) , 
have  been  observed  in  various  germanium  silicate  waveguides. 
The  further  understanding  of  these  nonlinear  optical  phenomena 
will  be  of  great  scientific  interest  and  may  lead  to  an 
important  technological  impact  in  fields  such  as  data  storage 
and  information  processing. 

2.2  Objectives  of  Present  Study 

The  objectives  of  this  work  are  to 
(1)  . explore  new  and  reliable  methods  to  fabricate  Ge02-Si02 
thin  films  which  can  be  used  as  planar  waveguides:  The  thin 

film  characteristics  were  studied  by  various  analytical 
techniques (e.g.  ellipsometry , wave-guiding  Raman  spectroscopy, 
XPS , electron  microscopy,  etc.).  Two  thin  film  deposition 
methods  were  successfully  applied  to  achieve  this  goal;  the 
details  of  thin  film  deposition  set-ups  and  procedures  are 
described  in  Chapter  III. 

(a)  . Sol-Gel  method:  This  involves  synthesis  of  the 

germanium  alkoxide  precursors,  sol  preparation,  film 
formation,  and  densif ication  of  the  films.  The  details  of 
this  study  and  the  various  characteristics  of  the  films  are 
presented  in  Chapter  IV.  With  this  method,  planar  waveguides 
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of  high  Ge02  content  can  be  easily  controlled  and  fabricated. 

(b) . Low  pressure  chemical  vapor  deposition (LPCVD) 
method:  This  involves  investigation  of  appropriate  film 
formation  conditions (e.g.  temperature,  flow  rate,  pressure, 
etc.).  By  using  this  method,  Ge02-doped  waveguides,  with  low 
Ge02  contents,  can  be  prepared  at  very  low  temperature 
compared  to  conventional  methods.  The  details  of  this  part  of 
the  study  are  described  in  Chapter  V. 

(2)  . characterize  waveguide  properties:  The  planar  waveguides 
were  characterized  by  a prism  coupling  technique.  the 
propagation  characteristics  (number  of  modes  supported  and 
their  propagation  constant,  propagation  loss,  etc.)  were 
measured.  Near-field  mode  profiles  of  the  waveguides  were 
measured  and  compared  to  the  calculated  results  using  a model 
calculations  for  planar  waveguide  with  a step-index  profile. 
The  details  of  experimental  set-up  are  given  in  Chapter  III 
and  the  results  are  discussed  in  Chapter  VI. 

(3) .  measure  nonlinear  optical  properties  of  the  films:  Two 
methods  were  used  to  characterize  the  nonlinear  optical 
property  of  the  films. 

(a) . Empirical  approach:  nonlinear  optical  coefficient 
n2  can  be  related  to  Abbe  number  and  linear  refractive  index 
through  an  empirical  equation  derived  by  Boling[Bol78] . The 
Abbe  number  of  the  films  was  calculated  by  using  an  additive 
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model  developed  by  Huggins  and  Sun(Hug43) . 

(b)  . Third  harmonic  generation (THG)  measurement:  This  is 
an  ideal  technique  for  the  measurement  of  relatively  small 
values  of  x®/  such  as  for  thin  films  deposited  on  a much 
thicker  substrate,  as  was  the  case  in  this  study.  By  measuring 
the  changes  in  the  THG  Interferometry  Fringes,  both  the 
magnitude  and  the  phase  of  x0)  can  be  obtained.  The  details 
of  the  experimental  set-up  and  results  will  be  described  in 
Chapter  III  and  VII,  respectively. 


CHAPTER  III 


EXPERIMENTAL 

3.1  Thin  Film  Preparation 

3.1.1  Sol-Gel  Method 

3. 1.1.1  Synthesis  of  germanium  alkoxides 

Because  of  the  extremely  moisture-sensitive  character- 
istics of  the  germanium  alkoxides,  all  reactions  and 
manipulations  were  carried  out  under  dry  nitrogen  using 
standard  Schlenk  techniques [ Shr86 ] . Ethanol,  2-propanol  and 
n-butanol  were  distilled  before  use.  2-methoxyethanol  and  n- 
hexane  were  used  as  received. 

Germanium  ethoxide,  Ge(OC2H5)4,  used  in  this  study  is 
prepared  using  an  ammonia  method  described  by  Bradley  et.al 
[Bra56 ] , or  is  commercially  available  from  Alfa  Products.  In 
a typical  reaction,  25  g of  germanium  tetrachloride,  GeCl4 
(Alfa)  , is  initially  mixed  with  160  ml  of  hexane,  250  ml  of 
ethanol  is  then  added  to  the  mixture  with  stirring  to  obtain 
a clear  solution.  Since  GeCl4  does  not  react  with  alcohol 
automatically,  the  reaction  needs  to  be  triggered  by  passing 
anhydrous  ammonia,  NH3 , through  the  reaction  vessel  to  act  as 
a proton  acceptor.  Germanium  ethoxide  forms  according  to  the 
reaction  expressed  as  follows: 
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GeCl4  + 4 C2H5OH  + NH3  ===>  Ge(OC2H5)4  + NH^l^ 

This  reaction  is  exothermic,  so  the  bubbling  of  ammonia  is 
continued  until  the  reaction  mixture  cools  down  to  room 
temperature.  Precipitates  of  NH4C1  are  removed  by  filtration. 
The  filtrate  is  then  refluxed  to  expel  excess  ammonia  and  then 
is  concentrated  by  removing  the  excess  solvent  and  alcohol 
under  ambient  pressure.  Vacuum  distillation  of  the  solution 
gives  a colorless  liquid  of  germanium  ethoxide  (boiling  point: 
-80  °C  at  3 torr)  . 

Various  germanium  alkoxides  with  higher  alkyl  groups, 
such  as  propyl,  butyl,  and  methoxyethyl  groups  etc.,  were 
synthesized  by  an  alcohol-interchange  method [Bra56 ] , using 
germanium  ethoxide  as  a starting  compound.  In  a typical 
reaction,  the  mixture  of  germanium  ethoxide,  higher  alcohol, 
and  hexane,  in  a molar  ratio  of  1:25:60,  was  refluxed  in  a 
column  packed  with  glass  beads  as  shown  in  Figure  3.1.  An 
azeotropic  distillation  was  followed  to  expel  the  azeotrope 
which  contains  21  wt%  of  ethanol  and  79  wt%  of  n-hexane  and 
boils  at  58 . 68°C  at  ambient  pressure.  The  formation  of  a low 
boiling  point  azeotrope  facilitates  the  reaction  in  the 
alcohol-interchange,  which  can  be  expressed  as  follows: 

Ge(OR')4  + 4 ROH  <===>  Ge  (OR) 4 + 4 R'OH 
where  R'is  -C2H5,  R is  a higher  alkyl  group,  such  as  -C3H7, 
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Figure  3 . 1 
reaction. 


Experimental  setup  for  alcohol-interchange 
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-C^,  or  -CH3OCH2CH2  etc..  After  removal  of  solvent,  the 
germanium  ester  was  vacuum  distilled  to  produce  high  purity 
alkoxide  for  the  preparation  of  Ge02-Si02  sols.  The 
replacement  of  the  ethyl  groups  in  germanium  ethoxide  with 
higher  alkyl  groups  is  expected  to  increase  the  stability  of 
the  germanium  alkoxide  against  hydrolysis/condensation 
reactions  due  to  steric  effects.  These  modified  germanium 
alkoxides  are  expected  to  increase  the  controllability  of 
hydrolysis  reactions  during  sol  preparation  and  thereby 
enhance  the  possibility  of  obtaining  stable  and  clear  sols  for 
the  coating  process. 

3. 1.1. 2 Sol  preparation 

The  Ge02-Si02  sols  were  prepared  by  mixing  a prehydrolyzed 
tetraethylorthosilicate,  TEOS (Fisher) , and  germania  sol 
together.  The  procedure  for  preparing  the  coating  solution  is 
shown  in  Figure  3.2.  In  a typical  procedure,  TEOS  is 
dissolved  in  ethanol  and  then  prehydrolyzed  using  water  in  the 
concentration  range  of  0.2-4  mole  of  water /mole  of  alkoxide. 
A small  amount  of  hydrochloric  acid,  HC1,  was  added  as 
catalyst.  Reactions  were  allowed  to  continue  for 
approximately  one  day  after  the  water /HC1  addition.  Germania 
sols  were  prepared  simply  by  dissolving  the  desired  amount  of 
germanium  alkoxide,  for  the  multicomponent  sol  of  pre- 
determined composition,  in  ethanol.  The  multicomponent  sol 
was  obtained  by  adding  germania  sol  to  silica  sol  under 
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Figure  3.2  Flow  diagram  for  the  preparation  of  binary 

Ge02-Si02  thin  films  using  sol-gel  method. 
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vigorous  stirring.  The  coating  solution  is  then  aged  for  8-12 
hours  before  the  coating  process.  Typical  formulations  for 
the  preparation  of  various  Ge02-Si02  sols,  using  germanium 
methoxy-ethoxide  as  precursor,  in  the  present  study  are  given 
in  Table  3.1.  Clear  and  stable  coating  solutions  may  be 
obtained  for  an  alkoxide  content  lower  than  0.5  M,  and  water 
to  alkoxide  molar  ratio (R)  smaller  than  3. 

3. 1.1. 3 Film  formation 

Thin  films  were  prepared  by  a dip  coating  method,  in 
which  a precleaned  fused  silica  substrate  was  immersed  into 
the  sol  followed  by  an  immediate  withdrawal  in  the  vertical 
direction.  Figure  3.3  shows  a schematic  of  the  dip-coater 
used  in  our  laboratory.  The  withdrawal  speed  used  in  the 
current  study  was  typically  2 to  7 cm/min. 

After  each  coating,  the  film  was  dried  in  air  for  a short 
time  before  applying  the  next  coating  to  allow  for  evaporation 
of  the  solvent  from  the  film.  After  several  dipping  cycles, 
the  film  was  then  heat-treated  in  the  temperature  range  of  300 
°C  to  1110°C  to  remove  the  organic  residues  as  well  as  to 
densify  the  film.  The  thickness  of  the  film  can  be  built  up 
by  repeating  the  multiple  dipping  and  firing  cycles. 

Several  organic  additives,  such  as  polyvinylpyrolidone 
(PVP  K-3  0 , GAF  Chemicals  Co.),  cp2(Daxad  cp2 , W.R.  Grace  & 
Co.),  and  hydroxypropylcellulose (HPC,  Aldrich  Chemicals  Co.), 
were  used  in  an  attempt  to  adjust  the  viscosity  of  the  coating 
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Table  3.1  Typical  Compositions  for  the  Preparation  of  Ge02-Si02 
sols. 


SAMPLE 

EtOH(ml) 

TEOS(ml) 

Ge(OR)4(ml) 

WATER  (ml) 

PH 

10GS 

20 

2 

0.4 

0.3 

1.58 

30GS 

20.5 

0.9 

0.7 

0.2 

3.63 

50GS 

21 

0.3 

0.6 

0.1 

1.91 

75GS 

19 

0.5 

2.6 

0.1 

1.07 

100GS 

21 

0 

1.1 

0 

2.71 
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substrate 

film 


controller 


Figure  3 . 3 


Schematic  diagram  of  the  dip-coating  system. 
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solution  to  enhance  thin  film  formation.  The  amount  of 
organic  additive  used  in  the  sol  is  usually  less  than  2 wt%. 

3.1.2  LPCVD  method 

A Low  Pressure  Chemical  Vapor  Deposition  System (LPCVD, 
Thermco)  was  used  to  prepare  samples  for  this  study.  The 
system  is  shown  schematically  in  Figure  3.4.  In  this  study  , 
~1  mm  thick  high  purity  grade  fused  silica  slides  were  used  as 
substrates.  The  substrates  were  cleaned  according  to  the 
following  procedures: 

(1)  rinsed  with  soaped  water,  and  immediately  dried  with  N2. 

(2)  sonicated  in  a deionized  water  bath  and  immediately  dried 
with  dry  N2. 

(3)  rinsed  in  electronic  grade  acetone  and  dried  in  dried  N2 
gas  jet. 

(4)  rinsed  in  absolute  ethanol  and  dried  immediately  in  dry 
N2  jet. 

The  precleaned  substrates  were  positioned  vertically  in  a 
quartz  boat  with  a separation  of  3~5  mm  between  substrates. 
The  boat  was  then  placed  in  the  center  of  the  three-zone 
resistively-heated  furnace,  which  acts  as  a hot  wall  reactor, 
and  the  temperature  was  allowed  to  equilibrate.  Usually  30-40 
substrates  were  loaded  in  the  boat  for  each  run. 

In  the  apparatus  used,  the  input  gas  flow  rates  were 
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Figure  3.4  Schematic  diagram  of  the  low  pressure  chemical 

vapor  deposition (LPCVD)  system. 
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measured  and  controlled  by  mass  flow  controllers (MFC)  . The 
gases  used  for  the  Ge02-Si02  film  deposition  were  0.2%  GeH4/Ar , 
SiH4/  and  02.  The  reagent  gases  were  introduced  into  the 
reactor  independently  through  different  injectors  located 
below  the  quartz  boat.  Nitrogen  gas  was  used  to  purge  the 
reactor. 

In  a typical  run,  the  reactor  was  first  pumped  down  to  a 
pressure  below  50  mtorr,  then  the  entire  reactor  was  back- 
filled with  nitrogen.  Before  starting  deposition,  the  N2 
background  atmosphere  was  replaced  by  02  gas  and  then  the 
hydride  gases  were  passed  into  the  reactor  at  predetermined 
flow  rates.  Films  were  usually  deposited  at  a reactor 
pressure  below  500  mtorr.  Table  3.2  shows  the  deposition 
parameters  used  for  the  Ge02-Si02  film  deposition  in  this 
study. 


3.1.3  Other  methods 

According  to  the  literature,  many  other  thin  film 
deposition  techniques  have  been  used  to  produce  Ge02-Si02 
planar (or  channel)  waveguides.  A brief  description  of  each 
technique  will  be  given  in  the  following  sections. 

3. 1.3.1  VAD  soot  particle  method  for  glazing  method^ 

Mori  and  Shimizu [Mor82 ] have  reported  the  use  of  a VAD 
soot  process (using  halide  precursors)  to  deposit  a doped 
silica  glass  particle  (Si02-Ge02-P205)  layer  as  a core  material 
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Table  3 . 2 LPCVD  Deposition  Parameters  for  the  Ge02-Si02  Film 
Formation. 


Sample 

Deposition 

temp.(°C) 

Pressure 

(mTorr) 

Flow  ratio 
0^(MH4) 

Dep.  rate 
(A/min.) 

Other 

LP1 

450 

245+20 

1.305 

13.75 

light  brown 

LP2 

LP3 

430 

450 

380+10 

330+10 

1.381 

1.494 

95.16 

95.67 

slightly 

hazy 

clear 
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on  a fused  quartz  substrate.  In  this  process,  silicon 
tetrachloride  is  used  as  the  main  raw  material.  Halides,  such 
as  GeCl4,  PC13,  and  BBr3,  are  used  as  dopants  which  adjust  the 
refractive  index,  softening  temperature,  and  thermal  expansion 
coefficient  for  the  deposited  glass. 

The  oxidation  reactions  take  place  in  a three-zone 
furnace,  which  acts  as  an  hot  wall  reactor.  The  gas  flows  and 
the  substrates  are  positioned  in  a horizontal  reactor 
configuration.  The  oxidation  reactions  are  usually  complete 
around  1200°C.  The  glass  particles'  deposition  rate  depends 
on  the  raw  material  supply,  furnace  temperature,  and  total  02 
gas  flow  rate.  Since  the  deposition  of  these  soot  particles 
occurs  via  a thermophoresis  process [Sim79,  Wal80],  the 
deposition  rate  depends  most  on  the  temperature  gradient  along 
the  gas  flow  direction. 

The  deposited  glass  particle  layer  was  made  transparent 
by  heating  above  the  glass  softening  temperature.  In  this 
process  P205  is  usually  added  to  lower  the  softening 
temperature  of  the  core  layer.  For  a glazing  temperature 
above  1450°C,  the  waveguide  loss  around  1 dB/cm  at  633  nm  was 
achieved.  The  typical  deposition  rate  of  this  process  is 
0.3-1  /m/min,  therefore  it  is  suitable  for  multimode  waveguide 
fabrication.  The  fabricated  waveguides  can  be  further 
processed  into  various  desired  devices,  by  going  through 
techniques  such  as  masking,  patterning  with  reactive-ion- 
etching, depositing  cladding  layer(e.g.  Si02-B203)  etc.. 
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By  using  this  process,  the  relative  refractive  index 
difference  between  core  layer  and  fused  silica  substrate  is 
typically  less  than  1%. 

3. 1.3. 2 Flame-hydrolysis  deposition (FHD)  method 

Flame  hydrolysis  deposition  is  one  of  the  promising 
technigues  for  mass-production  of  cost  effective  integrated 
optical  waveguide  devices  for  telecommunication  because  of  its 
compatibility  with  optical  fiber  technology.  Several  high- 
silica  based  waveguides  have  been  fabricated  by  this  method 
[Kaw83 , Kom90,  Max89,  Bar90]. 

Kominato  et . al.  [Kom90]  has  fabricated  Si02-Ge02-P205-B203 
waveguides  on  Si  substrates.  In  this  method,  a gaseous 
mixture  of  SiCl4-GeCl4-BCl3-PCl3  of  predetermined  flow 
proportion  is  fed  into  an  oxy/hydrogen  torch.  Fine  glass 
particles  synthesized  by  flame  hydrolysis  are  deposited  on  the 
substrates  placed  on  a turn  table.  During  the  deposition,  the 
torch  transverses  the  table  repeatedly  in  the  radial 
direction. 

Small  amounts  of  BC13  and  PC13  are  needed  in  the  raw 
material  gas  composition  to  lower  the  softening  temperature  of 
the  synthesized  glass  particles.  The  refractive  index  of  the 
glass  particles  is  controlled  by  the  GeCl4/SiCl4  flow  ratio. 
After  the  deposition,  the  substrates  with  porous  glass  layer 
are  heated  up  to  1250°C  for  one  hour  in  an  electric  furnace 
for  densif ication.  Deposition  parameters  have  to  be  optimized 
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to  obtain  the  proper  glass  composition,  particle  sizes,  and 
packing,  in  order  to  facilitate  the  subsequent  sintering  step 
to  achieve  bubble  free,  transparent,  low-loss  waveguides 
without  warping  the  substrate.  A typical  refractive  index 
difference  between  core  layer  and  the  substrate  is  0.2-1%. 
The  propagation  loss  of  the  waveguides  can  be  as  low  as 
0.01-0.03  dB/cm  at  1.55  jixm.  The  typical  deposition  rate  of 
this  process  is  as  high  as  -1  /xm/ transverse. 

3. 1.3. 3 Plasma-activated  chemical  vapor  deposition (PACVD) 
method 

Nourshargh  et.al.  used  a plasma  activated  technique  to 
produce  Ge02-Si02  and  Si3N4  waveguides [Nou8 6 , Nou89].  Because 
it  is  a plasma  process,  the  introduction  of  a wide  range  of 
dopants  into  the  silica  guide  is  possible.  Like  other  halide 
CVD  processes,  this  method  uses  SiCl4,  GeCl4,  and  02  as 
reactants.  The  reactor  consists  of  a heated  reaction  zone 
maintained  at  a low  pressure («0 . 75  Torr) , into  which  are  fed 
vapors  of  SiCl4,  GeCl4,  and  oxygen  through  independent 
injectors  from  one  end  of  the  reactor.  Ar  gas  is  used  to 
purge  the  reactor  and  to  form  the  plasma,  because  of  it  is 
easily  ionizable.  A microwave  generator («500  W at  f=2.45  GHz) 
is  used  to  excite  a plasma  in  the  reaction  zone.  A fused 
silica  substrate  is  resistively  heated  and  situated 
horizontally  within  the  plasma  region.  Over  the  surface  of 
the  substrate,  a layer  of  glassy  doped  silica  is  deposited. 
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Because  this  is  a low  pressure  process,  formation  of 
particulates  is  avoided  and  homogeneous  waveguiding  layers  are 
grown  on  the  substrate  by  oxidation  of  the  halide  reactants. 
The  detailed  reaction  mechanisms  are  not  well  understood. 

The  refractive  index  of  the  film  can  be  varied  by 
adjusting  the  ratio  of  the  flow  rates  of  GeCl4  and  SiCl4.  The 
thickness  of  the  deposited  layer  is  a function  of  both  the 
deposition  time  and  total  chemical  flow  rates.  The  waveguides 
may  be  fabricated  at  a temperature  ranging  from  1100°C  down  to 
room  temperature.  Typically  a waveguide  prepared  at  1100  C 
with  the  flow  rates  of  02  through  the  SiCl4  and  GeCl4  bubblers 
being  50  and  30  seem  respectively.  The  deposition  rate  is 
around  0.5  jum/min  and  a refractive  index  change  about  0.011 
was  reported.  The  propagation  loss  is  about  0.3  dB/cm  at 
0.633  nm. 

This  technique  is  capable  of  operating  at  low 
temperatures (=s400°C)  , however  the  waveguides  produced  at  low 
temperatures  possess  a higher  loss  rate (-3  dB/cm) . 

Single  mode  pure  Ge02  waveguides  were  successfully 
deposited  at  a deposition  rate  of  0.03  /nm/min  and  with  a loss 
of  approximately  1~2  dB/cm. 

3. 1.3. 4 Reactive  sputtering  method 

Both  Ge02-Si02  and  pure  Ge02  waveguides  have  been  prepared 
by  using  rf  reactive  sputtering  deposition  technique  [Yin82, 
Hua85 ] . In  this  method,  gases  Ar  and  02  are  introduced  into 
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the  vacuum  chamber  through  the  control  of  two  gas  flowmeters. 
Proper  flow  ratio  between  Ar  and  02  is  needed  to  compensate 
for  the  inherent  oxygen-deficient  characteristic  of  the  oxide 
films  prepared  by  this  vacuum  deposition  technique.  With 
insufficient  02  present  in  the  vacuum  chamber  during  film 
deposition,  a highly-reduced  film  with  a dark  yellow 
appearance  will  be  obtained. 

To  start  the  deposition  process,  the  sputtering  system  is 
pumped  down  to  a pressure  in  the  region  of  4x10-6  Torr  and 
then  backfilled  with  gases  of  Ar  and  02(02/Ar=l: 1)  to  a total 
gas  pressure  of  2x10-3  Torr.  Colorless  low-loss  Ge02  thin 
film  optical  waveguides  were  prepared  using  rf  reactive 
sputtering  with  a Ge02  target  and  a low  deposition  rate  (< 
12.7  A/min) . The  deposited  films  were  annealed  at  250°C  for 
6 hours.  Films  with  a thickness  up  to  3800A  have  been 
prepared  with  a waveguide  attenuation  of  <0.7  dB/cm  at 
0.6328 /xm.  The  average  refractive  index  of  the  Ge02  films  was 
measured  to  be  1.6059  at  X=546lA  by  an  ellipsometry  technique. 

3 .1.3. 5 Laser-induced  CVD  method 

Photoinduced  chemical  vapor  deposition  has  been  used  to 
form  various  kinds  of  films [Sol81,  Boy82,  Joh82,  Sol83], 
including  Ge02-Si02  glass  films [Tat86 ] . In  this  method,  a 
pulsed  193nm  ArF  excimer  laser  provided  a rectangular  beam 
focused  by  a cylindrical  lens  telescope.  The  laser  was 
operated  at  a 3 Hz  repetition  rate.  This  optical  setup 
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provided  a sheet  of  photons  parallel  to  and  at  a small 
distance (typically  1mm)  above  the  Si02  substrate  surface. 
Substrates  were  positioned  horizontally  and  were  maintained  at 
temperatures  between  200  to  300°C  by  resistance  heating. 
Gaseous  reactants  were  5%  SiH4/N2,  5%  GeH4/N2,  and  pure  N20. 
The  input  window  was  purged  with  N2  or  He  gas  to  prevent 
opaque  deposition  from  forming,  but  the  output  window  was  not 
purged.  Typical  gas  flow  rates  used  for  the  deposition  are  Si 
source: 200  seem,  Ge  source:  0-50  seem,  N20:  1000-1500  seem, 
N2: 800-1500  seem,  and  He: 800-2000  seem.  In  a typical  run,  the 
partial  pressure  of  N20  ranges  between  4 to  8 Torr. 

The  flow  ratio  of  the  reactant  gases  was  determined  based 
upon  a consideration  of  the  absorption  cross-section  of  the 
laser  light  and  the  quantum  yield  of  the  photodissociation  of 
each  reactant  used.  Germane  gas  appears  to  have  a similar 
cross-section  to  N20,  and  has  a nearly  ten  times  larger  cross- 
section  compared  to  that  of  silane  gas.  N2  is  almost 
"transparent”  in  this  reaction.  Furthermore,  based  upon  these 
consideration,  it  is  believed  that  Si02  is  mainly  produced  by 
a direct  oxidation  reaction  between  SiH4  and  photodecomposed 
oxygen  from  N20.  On  the  other  hand,  Ge02  is  generated  by  a 
reaction  between  photodecomposed  oxygen  from  N20  and  Ge  atoms 
from  GeH4.  Therefore,  in  this  deposition  method,  it  is 
thought  that  the  Ge02-generating  process  is  limited  by  the 
GeH4  supply  and  that  the  Si02-generating  process  is  limited  by 
the  O supply. 
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With  this  method  single-mode  Ge02-Si02  waveguides  have 
been  prepared  and  a propagation  loss  as  low  as  1 dB/cm  at 
0 . 6328/im  for  low  Ge02  doped  film  was  achieved.  However,  the 
propagation  loss  increases  drastically  for  higher  Ge02  content 
films . 

3.2  Characterization 
3.2.1  Infrared  Spectroscopy 

Fourier-transformed  infrared (FTIR)  spectroscopy  (Nicolet 
20SXB  FT-IR  Spectrometer)  was  used  to  examine  the  synthesized 
alkoxides,  Ge02-Si02  gels  and  the  thin  film  samples. 

Transmission  spectra  of  germanium  alkoxides  were  obtained 
using  a liquid  cell (circle  cell,  Spectra  Tech  Inc.)  with  a 
ZnSe  crystal  which  gives  a transmission  range  from  4000  to  600 
cm'1.  The  spectra  were  obtained  in  a nitrogen  atmosphere  with 
32  accumulation  scans. 

Transmission  spectra  were  also  obtained  from  the  Ge02-Si02 
gels,  heat-treated  at  different  temperatures,  using  a diffuse 
reflectance  stage (Spectra  Tech  Inc.)  and  128  accumulation 
scans.  The  gel  was  ground  to  fine  powder,  packed  loosely  and 
leveled  to  the  top  of  the  stage  for  the  measurement. 

The  transmission  spectra  of  thin  films  produced  from  both 
LPCVD  and  sol-gel  methods  were  obtained  by  using  an  attenuated 
total  ref lectance (ATR)  stage (Spectra  Tech  Inc.)  with  a KRS-5 
crystal (a  mixed  crystal  of  ThI2-ThBr2  eutectics).  The  KRS-5 
crystal  was  sonicated  in  a methanol  bath  for  ten  minutes  and 


63 


dried  in  a N2  jet  before  use.  Each  thin  film  sample  was 
cleaned  with  ethanol  and  held  against  the  crystal  for  the 
measurement.  128  scans  were  used  to  obtain  each  spectrum. 

3.2.2  Thermal  Analysis 

Thermal  behavior  of  the  gels,  including  the  behavior  of 
organic  burn-out,  was  studied  using  differential  thermal 
analysis (DTA) , and  thermogravimetric  analysis (TGA) (Model  STA- 
409,  Netzsch  Inc.,  Exton,  PA.).  The  gels  were  dried  at  120°C 
under  vacuum  for  8 hours  and  ground  to  fine  powder.  Powdered 
samples  were  preferred  in  this  study  not  only  to  assure  a good 
thermal  contact  between  sample  and  crucible  to  increase  the 
accuracy  of  DTA  measurement,  but  also  to  ease  drying  and 
removal  of  physical  adsorbed  water  which  might  be  easily 
trapped  within  the  pores  in  a bulk  sample.  For  these 
measurements,  20  mg  of  the  sample  was  placed  in  an  Al203 
crucible  and  heated  up  to  900°C  in  02  atmosphere  at  a heating 
rate  of  5°C/min. 

3.2.3  ICP  Solution  Analysis 

Inductively  coupled  plasma (ICP)  emission  spectroscopy 
(Instrumentation  Lab.  Plasma  200)  was  used  to  study  the 
composition  of  the  gels  prepared  in  present  study.  After 
removing  excess  solvent  from  the  gel,  solution  samples  were 
prepared  by  dissolving  a predetermined  amount  of  wet  gel  in 
100  ml  of  deionized  water. 


A few  drops  of  hydrofluoric 
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acid(HF)  was  added  for  dissolving  gels  of  lower  germania 
content,  e.g.  10GS.  The  solution  was  then  diluted  for  the 
measurement.  Some  common  contaminants  of  the  silicate  system, 
such  as  Na,  Mg,  Al,  Fe,  and  Cu,  were  analyzed  to  check  the 
purity  of  the  gels.  All  the  samples  and  standard  solutions 
were  prepared  just  before  conducting  the  measurement. 

In  the  measurement,  the  solution  sample  was  fed  at  a 
constant  rate  into  the  system  by  a peristaltic  pump.  The 
solution  was  then  atomized  by  a nebulizer  and  sprayed  through 
an  Ar  plasma  area,  sustained  by  a rf  coil.  The  electrons  of 
the  elements  in  the  solution  were  excited  by  the  plasma, 
followed  by  various  atomic  transition  processes  which  emit 
characteristic  photons  corresponding  to  the  excited  state 
energy  structure  of  different  elements.  The  light  was 
detected  by  an  assembly  of  monochromater  and  photomultiplier 
tube (PMT) . 

3.2.4  X-ray  Photoelectron  Spectroscopy (XPS) 

XPS  (Perkin-Elmer  PHI  5100  ESCA  System)  was  used  to 
analyze  the  composition  of  the  thin  films  prepared  by  both 
sol-gel  and  LPCVD  processes.  A magnesium  anode (Mg  Ka=1253.6 
eV)  was  used  as  X-ray  source.  Photoelectrons  were  collected 
at  a take  off  angle  of  45°  and  spatially-dispersed  with  a 
hemispherical  electron  energy  analyzer  and  detected  with  a 
multichannel  detector. 

All  of  the  samples  were  examined  in  the  survey  scan. 
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which  is  the  low  resolution  high  sensitivity  mode,  with  a pass 
energy  of  89.45  eV,  a resolution  of  0.35  eV/step,  and  an 
accumulation  time  of  50  ms/ step,  to  determine  whether  there 
were  any  impurities  present  on  the  analyzed  surface  after 
achieving  a vacuum  of  10'9  Torr.  The  regions  of  interest  were 
studied  with  high  resolution  multiplex  scans,  with  a pass 
energy  of  8.95  eV,  a resolution  of  0.05  eV/step,  and  a 
accumulation  time  of  100  ms/step.  The  peak  areas  of  Ge3d  and 
Si2p  band  were  calculated  to  represent  the  atomic 
concentrations  of  Ge  and  Si  respectively  in  the  film. 

3.2.5  Electron  Microscopy 

Scanning  electron  microscopy (SEM)  (JOEL  JSM-6400  Scanning 
Microscope)  was  used  to  study  the  surface  morphology  of  the 
waveguides . 

Transmission  electron  microscopy  (JEOL  200CX  TEM)  was  used 
to  investigate  the  microstructure  of  the  films.  Plane  samples 
were  prepared  for  this  study  by  mounting  the  sample  on  a fused 
silica  slide  with  crystal  bond.  The  sample  was  then  polished 
from  the  side  of  the  substrate  to  a thickness  of  -150  /xm. 
Several  disk  samples  were  cut  from  the  polished  thin  sample. 
The  disk  was  then  dimpled (Dimpler,  VCR  Group,  San  Francisco, 
CA)  to  a thickness  of  around  10  jum.  An  ion-milling (Gatan 
Dual  Ion  Mill  Model  600)  step  was  followed  to  perforate  the 
sample  for  TEM  observation.  Selected  area  electron 
diffraction  (SAD)  was  used  to  study  the  possible  formation  of 
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small  crystallites  in  the  film. 

3.2.6  X-ray  Diffraction (XRD)  Measurement 

Phase  identification  was  performed  using  room  temperature 
X-ray  powder  diffraction  with  Cu  Ka  radiation  (Phillips  APD 
3720  X-ray  Diffractometer) . The  Cu  Ka  radiation  was  generated 
using  a fine  focus  tube  and  a graphite  monochromator  at  40  KV 
and  20  mA.  2d  peaks  between  5°  to  100°  were  scanned  at  a rate 
of  3°/min. 

3.2.7  UV-VIS  Spectroscopy 

Optical  characteristics  of  the  LPCVD  thin  films  in  the 
ultra-violet  range  were  investigated  by  using  a UV-VIS-NIR 
spectrophotometer  (Perkin-Elmer  Lambda  9) . 

3.2.8  Ellipsometric  Measurements 

An  ellipsometer (Gaertner  L117  Manual  Ellipsometer)  was 
used  to  measure  the  refractive  index  and  the  thickness  of  the 
deposited  films.  Reflection  ellipsometric  measurements 
involve  illuminating  the  surface  of  a sample  with  monoch- 
romatic light  of  known  wavelength  and  polarization  at  a fixed 
angle  of  incidence,  and  then  analyzing  the  change  of  polari- 
zation that  takes  place  on  reflection  due  to  the  difference  in 
attenuation  amplitude  and  phase  shift  experienced  by  the 
p (parallel)  and  s (perpendicular)  components. 
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Figure  3.5  shows  the  schematic  of  the  system  set-up.  The 
light  source  used  in  this  study  is  a He-Ne  laser  (632 . 8nm)  . The 
beam  is  circularly  polarized  at  the  laser  output  and  passes 
through  a polarizer.  In  passing  through  the  polarizer  the 
beam  polarization  is  converted  from  circular  to  linear.  This 
constant  intensity,  linearly  polarized  beam  is  then  converted 
to  one  of  elliptical  polarization  by  a quarter-wave  compen- 
sator (wave  plate  oriented  at  +4  5°)  . Upon  illuminating  the 
surface  of  the  sample,  the  reflected  light,  with  its  polari- 
zation altered  by  the  properties  of  the  sample,  passes  through 
an  analyzer  and  an  optical  interference  filter  and  is  sensed 
by  a photodetector.  Certain  azimuth  settings (P)  of  the  pola- 
rizer cause  the  light  reflected  from  the  sample  to  become 
completely  linearly  polarized.  When  the  polarizer  is  rotated 
to  one  of  these  settings,  the  analyzer  can  then  be  rotated  to 
a position (A)  where  almost  no  light  reaches  the  detector, 
thereby  indicating  the  lowest  reading  on  the  extinction  meter. 
This  is  the  condition  at  which  polarizer  and  analyzer  readings 
are  recorded. 

When  an  elliptically  polarized  light  is  reflected  from  a 
surface,  its  polarization  direction  is  changed  by  the  refra- 
ctive index  and  absorption  coefficient  of  the  reflecting 
surface.  If  the  surface  is  a thin  film,  as  shown  in  Figure 
3.6  of  a three  medium  system,  then  the  film  thickness  also 
enters  into  the  reflection  equation [Azz77 ] . The  Fresnel 
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Figure  3.5  Optical  system  functional  diagram  for 

reflection  Ellipsometry  measurement [Gae86 ] . 


the 
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Figure  3.6  Oblique  reflection  and  transmission  of  a plane 
wave  by  an  ambient (0) -film (1) -substrate (2)  system  with 
parallel-plane  boundaries [Aza79 ] . 
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equations  are  shown  as  follows: 


roi„  = 


WjCOS0o  -NoCOS0j 
W,COS0O  +WOCOS0, 


(3.2.1) 


W2cos0j  -W,cos02 
1?p  W2cos0,  +Nlcos(p2 


(3.2.2) 


roi,= 


Wocos0o  -W2cos0, 

NqCos<Pq  +N]COS0j 


(3.2.3) 


ri2,= 


_ WjCOS0j  -W2COS02 
N,COS0,  +W2COS02 


(3.2.4) 


where  r01p,  r01f,  r12p/  r12s  are  interface  Fresnel  reflection 

coefficients.  In  general  N;  is  a complex  number,  but  for  a 
non-absorbing  medium  it  is  real  number  nj.  N0  is  the 
refractive  index  of  the  incident  medium (usually  air) , Nj 
corresponds  to  the  unknown  film  and  N2  to  a known  substrate. 
0i  are  the  angles  in  each  medium  and  they  are  related  by 
Snells'  law. 


nosin0o  =nlsincpl  =n2sin02 


(3.2.5) 
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Figure  3.6  shows  that  multiple  reflections  will  occur  in  the 
thin  film,  thus  the  reflected  beam  is  a sum  of  infinite 
contributions  and  the  complex  amplitude  reflection  coeffi- 
cients (for  both  polarizations)  can  be  found  to  be[Azz77] 


R = 


i01  p 


+ r 


12pc 


, -tv 


1 +r. 


01/>ri2/>e 


-m 


(3.2.6) 


*s  = 


r01  s+ri2se'i2t> 
1 +r0l,ri2se~i2^ 


(3.2.7) 


where  /3  is  the  phase  change  that  the  multiply-reflected  wave 
inside  the  film  experiences  as  it  trverses  the  film  from  one 
boundary  to  the  other.  It  is  given  by 


P=2n( 


d 


)nlcos(pl 


(3.2.8) 


where  d is  the  film  thickness,  and  X the  wavelength  of  the 
incident  light. 

The  complex-amplitude  reflection  coefficients  may  be 
written  in  terms  of  their  absolute  values  and  angles [Azz77 ] : 


i A. 


I A. 
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From  measurements  of  the  incident  and  reflected  polarizations 
(i.e.  ellipsometry)  the  ratio 


P 


of  the  overall  system  is  determined.  Expressing  p in  terms  of 
the  ellipsometric  angles  $ and  A, 


P 


•A 

= tan$  e‘A  = — i 


(3.2.9) 


where  $ and  A are  ellipsometric  angles:  amplitude  ratio  change 


tan$  = 


(3.2.10) 


and  change  in  phase  difference 


A=Ap-A,  (3.2.11) 

In  the  ellipsometric  measurement,  two  sets  of  polarizer 
and  analyzer  readings,  i.e.  P,,  P2,  A,,  A2,  were  taken.  They 
are  related  by  the  following  relations 
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P2=Pl+90 


(3.2.12) 


7^  = 180  -A, 


(3.2.13) 


Thus  by  tracing  the  two  sets  of  data  systematic  errors  are 
reduced  and  better  statistics  are  obtained.  These  measured 
data  can  be  used  to  calculate  ellipsometric  angles  using  the 
following  relations: 


[180-  (Aj-A,)  ] 
2 


(3.2.14) 


A = 360  - (P,  +P2) 


(3.2.15) 


Parameter  Psi($)  varies  from  0 to  90,  while  Delta(A)  can 
assume  values  from  0 to  360. 

Equations  (3 . 2 . 7) - (3 . 2 . 11)  can  not  be  solved  analytically 
for  n,  and  d,  but  numerical  methods  have  to  be  applied.  In 
this  study,  a computer  program  was  developed  by  Ochoa  to 
calculate  the  n,  and  d.  By  using  the  program,  a quick 
determination  of  n,  can  be  made  by  graphically  comparing  the 
A and  $ measured  to  calculated  values. 

Once  N,  is  found,  d can  be  calculated  using  the  following 
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equation [Azz77] 


d = d0+  [1  {Ni-Nosin2(p0)-ll2]m\  m = 0 , ±1, 2 , 3 . . (3.2.16) 

Thus,  the  thickness  calculated  is  not  an  unique  value. 
However,  if  one  knows  the  approximate  thickness,  so  that  m may 
be  defined,  then  an  accurate  value  can  be  obtained.  The 
thickness  of  the  films  were  checked  by  using  other  techniques, 
such  as  prof ilometry , SEM,  and  interferometry,  whenever 
necessary. 

3.2.9  Raman  Spectroscopy 
3.2.9. 1 Raman  scattering 

In  a Raman  scattering  experiment,  visible  light  from  a 
polarized  monochromatic  laser  is  passed  through  the  sample. 
A small  fraction  of  light  is  scattered  by  atoms,  and  some  of 
this  fraction  of  light  interacts  with  the  sample  in  such  a way 
as  to  excite  a vibrational  mode.  When  this  occurs  the  energy 
of  the  scattered  light  is  reduced(e.g.  for  Stokes)  by  an 
amount  corresponding  to  the  energy  of  the  vibrational 
transition.  This  type  of  inelastic  scattering  is  known  as 
Raman  scattering.  The  energy  of  the  scattered  light  is 
analyzed  using  a spectrometer.  Raman  lines  appear  as  weak 
peaks  shifted  in  energy  from  the  Rayleigh  line.  The  positions 
of  these  Raman  peaks  about  the  incident  line  correspond  to  the 
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frequencies  of  Raman  active  vibrations  as  illustrated  as  in 
Figure  3 . 7 [ McM8  8 ] . 

3. 2. 9. 2 Measurement  on  thin  film  samples 

Raman  scattering  is  a useful  method  to  get  information 
about  vibrational  or  rotational  modes  in  molecular  systems. 
The  scattered  flux  which  enters  the  spectrometer  is 
proportional  to  the  intensity  of  the  exiting  beam  and  to  the 
number  of  molecules  which  are  included  in  the  illuminated 
sample  volume  which  is  defined  by  the  image  of  the  entrance 
slit  of  the  spectrometer. 

In  order  to  study  efficiently  Raman  scattering  in  a thin 
film,  the  scattering  from  the  film  must  be  maximized,  while 
scattering  from  the  substrate  must  be  minimized.  Several 
Raman  scattering  measurement  techniques,  such  as  the  waveguide 
method [ Lev74 ] , the  total  reflection  coupling  method [ Dup74 ] , 
and  the  Raman  microprobe [ She88 ] etc.,  have  been  used  to  study 
thin  film  samples.  Among  them,  the  waveguiding  method  is  the 
most  effective  and  has  been  successfully  used  in  the  studies 
of  various  waveguiding  films[Ste88,  Tal88,  Levy74,  Gal78]. 

In  this  study,  the  Raman  spectra  of  Ge02-Si02  films  were 
obtained  by  using  the  waveguiding  technique.  The  Raman 
spectra  were  obtained  by  using  the  experimental  set-up  as 
shown  in  Figure  3.8.  An  Ar+  ion  laser  beam (5 14  nm)  was  prism- 
coupled  into  the  thin  film  by  using  a high  index  glass  prism 
( n= 1.8054) . The  polarization  of  the  laser  beam  was  parallel 
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Figure  3.7  (a)  The  Raman  scattered  light,  (b)  The  Raman 
spectrum  consists  of  a Rayleigh  scattering  peak,  and  the  much 
weaker  Raman  shifted  lines,  (c)  The  energy  level  diagram  for 
the  Raman  scattering  [McM88]. 
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to  the  plane  of  the  film,  i.e.  TE  mode,  and  the  power  incident 
on  the  coupling  prism  was  typically  400mW.  The  scattered 
light  was  collected,  in  a nearly  90°  configuration  with 
respect  to  the  direction  of  propagation  of  guided  wave  in  the 
film,  into  the  slit  of  the  spectrometer  (Instruments  S.A. 
Ramanor  U1000  Monochomator) . The  spectrometer  has  a computer- 
controlled  1.0m  double  monochromator  equipped  with  two 
gratings ( 1800  grooves/mm)  and  a photomultiplier  tube.  For 
this  study,  the  spectral  resolution  of  the  system  is  5 cm'1, 
slit  size  is  500  /im  and  a integration  time  of  20  sec.  is  used. 

For  the  films  that  do  not  possess  waveguiding  quality, 
the  Raman  spectra  were  obtained  by  using  a microprobe 
configuration  as  shown  in  Figure  3.9.  A small  portion  of  the 
green  light  was  reflected  by  a beam  spliter(BS) , passing 
through  a microscope  objective (40x)  and  directed  to  the  film 
and  the  substrate.  The  scattered  light  was  then  collected  by 
the  same  microscope  objective.  The  objective  was  finely 
adjusted  to  focus  on  the  surface,  so  that  the  signal  from  the 
substrate  was  minimized  and  the  film  signal  was  maximized. 
The  power  of  the  laser  light  directed  to  the  sample  is  15  mW, 
the  spectral  resolution  used  was  2 cm'1,  integration  time  was 
5 sec.,  and  a silt  size  of  500  jum  was  used. 

3.2.10  Planar  Waveguide  Characterization 
3.2.10.1  Prism  coupling 

Prism  coupling  is  one  of  the  primary  diagnostic  tools  for 
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studying  the  propagation  characteristics  of  dielectric 
waveguides.  The  theory  of  prism  coupling  has  been  treated 
extensively  in  the  literature [Tie71 , Ulr70,  Tie70],  and  the 
basic  concepts  will  be  summarized  here.  In  its  simplest  form, 
the  coupler  makes  use  of  a high  refractive  index  prism  placed 
in  close  proximity  (air  gap<  X/4)  to  a planar  waveguide  as 
shown  in  Figure  3.10.  In  order  to  excite  all  possible 
waveguide  modes,  the  phase  matching  condition  between  prism 
and  film  has  to  be  met.  This  reguire  the  index  of  prism(np) 
to  be  larger  than  that  of  the  film(nf)  . 

When  an  optical  beam  passing  through  the  prism  is 
incident  upon  its  bottom  face  at  an  angle  exceeding  the 
critical  angle,  the  beam  will  be  totally  reflected  at  the  base 
of  the  prism.  Because  of  the  total  reflection,  the 
electromagnetic  field  in  the  prism  is  a standing  wave  that  may 
be  described  by  an  exponentially  decreasing  function  below  the 
base  of  the  prism.  The  part  of  the  field  that  extends  below 
the  prism  base  is  called  the  evanescent  field.  Since  it 
decreases  rapidly  away  from  the  prism,  it  does  not  represent 
free  radiation.  The  evanescent  fields  that  extend  below  the 
prism  base  can  penetrate  into  the  waveguide  by  a process 
called  optical  tunneling.  These  fields  are  capable  of 
transferring  power  between  the  incident  beam  and  a waveguide 
mode.  By  appropriate  choice  of  the  angle  of  incidence  and 
proper  coupler  design,  a significant  portion  of  the  power  in 
the  incident  beam  may  be  transferred  into  a single  chosen 
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Figure  3.10  Overlap  of  evanescent 
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waveguide  mode(called  excitation  of  a mode).  Further,  by 
reciprocity,  if  a second  identical  prism  is  placed  in  close 
proximity  to  the  waveguide  at  some  distance  away  from  the 
input  prism,  each  propagating  mode  will  be  coupled  out  of  the 
guide  at  an  angle  that  is  characteristic  of  the  particular 
mode. 

3.2.10.2  Mode  index  measurement 

By  utilizing  a prism  coupling  technique,  the  light  may  be 
coupled  into  the  waveguide,  as  shown  in  Figure  3.11.  Because 
the  transverse  resonance  condition[Tie71]  has  to  be  satisfied 
in  order  to  guide  the  light,  as  has  been  described  in  section 
1.3.1,  the  coupling  occurs  at  specific  discrete  angles.  The 
maximum  coupling  takes  place  when  the  phase  matching  condition 
between  the  light  within  the  prism  and  the  specific 
propagation  mode  within  the  waveguide  is  met.  This  will  be 
discussed  in  more  details  in  section  6.2.  The  mode  index  of 
the  propagating  mode  may  be  obtained  by  measuring  the  incident 
angle  of  the  laser  beam  into  the  coupling  prism,  which 
efficiently  excites  that  specific  mode  in  the  waveguide. 

To  measure  the  angle  of  incidence  to  the  coupling  prism, 
a precision  rotation  stage (Klinger) , as  shown  in  Figure  3.12 
is  used.  The  waveguide  (substrate  refractive  index=1.4579  at 
632. 8nm)  was  held  in  a holder  with  a high  index  glass 
prism (SF-6  n=  1.8052)  closely  attached  to  its  surface.  The 
waveguide  was  placed  so  that  the  corner  of  the  prism  was  at 
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Figure  3 . 11 


Ray  diagram  for  the  prism  coupler. 
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Figure  3 . 12 
P:  polarizer. 


Experimental  setup  for  mode  index  measurement, 
L:  lens,  and  D:  photodetector. 
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the  center  of  the  rotation  stage.  A He-Ne  laser (X=632 . 8nm, 
20mW,  Spectra-Physics)  was  focused  into  the  prism  with  a lens 
of  long  focal  length.  The  point  of  maximum  coupling, 
determined  as  the  input  incident  angle  is  slowly  increased  by 
rotating  the  stage.  The  light  passing  through  the  film  is 
collected  at  the  end  of  the  waveguide  by  a microscope 
objective (20x)  and  detected  by  a photodetector (Newport  818- 
SL) . The  incident  angle  that  gives  the  maximum  light  output 
is  recorded. 

3.2.10.3  Propagation  loss  measurement 

There  are  two  types  of  losses  which  contribute  to  the 
propagation  loss  in  optical  waveguides:  absorption  losses  and 
scattering  losses. 

The  absorption  losses,  at  various  specific  wavelengths, 
are  primarily  caused  by  the  presence  of  impurities;  especially 
the  transition  metal  and  hydroxyl  ions.  These  losses  may  be 
kept  small  by  controlling  the  composition  of  the  waveguide 
materials.  The  scattering  losses  are  of  two  types:  volume 
scattering  and  surface  scattering.  Volume  scattering  is 
caused  by  imperfections  such  as  pores,  inclusions (i. e.  density 
and/or  composition  fluctuation  regions)  etc.,  within  the 
volume  of  the  waveguide.  The  type  of  loss  depends  very 
strongly  on  the  relative  size  of  the  imperfection  in 
comparison  with  the  wavelength  of  light  in  the  material  as 
well  as  its  concentration.  Surface  scattering  loss  can  be 
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significant  particularly  for  the  fundamental  mode  in  an 
asymmetrical  waveguide  because  most  of  the  guided  energy  in 
this  case  is  confined  to  the  region  that  is  closer  to  the 
surface.  The  estimation  of  propagation  loss  is  important  in 
the  waveguide  design  and  fabrication,  because  from  an 
application  point  of  view,  the  waveguide  should  have  as  low  an 
insertion  loss  as  possible,  and  propagation  loss  is  one 
component  of  the  insertion  loss. 

Several  methods  of  measuring  the  attenuation  have  been 
reported [Tie71] . First,  laser  light  is  coupled  into  the  film, 
e.g.  with  a prism  coupler,  and  then  it  propagates  freely  in 
the  form  of  a streak  in  the  film.  Usually  , for  an  non- 
absorbing material,  as  long  as  the  attenuation  is  1 dB/cm  or 
more,  the  attenuation  in  the  film  can  be  determined  by 
measuring  the  intensity  of  the  scattered  light  as  a function 
of  distance  along  the  guided  beam.  This  method  has  been 
successfully  used  for  measuring  attenuation  as  low  as 
0. 3dB/cm[Nou85] . The  intensity  of  the  scattered  light  is  a 
product  of  the  local  intensity  of  the  scattering  efficiency. 
Because  the  latter  is  not  necessarily  homogeneous  along  the 
guide,  the  results  from  this  technique  become  ambiguous  for 
lower  losses. 

The  most  widely  used  method  for  measuring  the  attenuation 
of  a waveguide  is  to  measure  the  transmitted  power  as  a 
function  of  waveguide  length [Tie71 ] . In  this  method,  all  the 
energy  flowing  through  the  guide  must  be  collected  in  the 
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detector.  If  one  does  not  want  to  cut  the  waveguide  and 
polish  the  end  face,  one  has  to  use  an  output  coupling  prism 
to  extract  the  energy  out  of  the  guide  at  different  positions. 

The  plot  of  output  power  vs.  distance  has  a negative 
slope  as  some  power  is  either  absorbed  by  the  impurities  or 
lost  by  scattering.  The  magnitude  of  the  slope  is  the 
attenuation  of  the  waveguide.  To  have  a reliable  measurement, 
however,  it  is  necessary  to  keep  the  same  coupling  for  each  of 
the  prism  at  all  positions  of  the  output  prism.  This 
condition  is  rather  restrictive  as  the  input  prism  tends  to 
get  disturbed  while  moving  the  output  prism. 

Various  methods  have  been  used  to  assure  a high 
extraction  efficiency:  Weber  et  al.[Web73]  used  index  matching 
liquid  between  prism  and  guide.  But  when  several  modes  are 
present,  the  output  light  of  individual  mode  tends  to  overlap, 
which  makes  difficult  measurements  of  individual  mode  loss. 

When  no  index  matching  liquid  is  used,  the  prisms  must  be 
strongly  clamped  to  the  guide,  and  care  must  be  taken  not  to 
disturb  the  input  coupling,  not  to  spoil  the  guide  surface  at 
the  contact  point,  and  not  to  locally  change  intermodal 
coupling. 

In  this  study,  an  improved  three  prism  loss  measurement 
method [Won80 ] , in  which  the  results  are  independent  of  the 
prism  coupling  coefficients,  is  used  to  measure  the 
propagation  loss.  The  experimental  setup  is  shown  in  Figure 
3.13.  Prism  1,  2,  and  3 are  independently  clamped  to  the 
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Figure  3.13  Experimental  setup  for  three-prism  loss 
measurement;  PI,  P2  and  P3  is  glass  prism,  D is  photodetector. 
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waveguide.  Small  glass  prisms(base  length  5 mm,  n=1.85)  are 
used,  The  air  gap  between  the  film  and  the  prism  is  adjusted 
to  be  small  compared  to  the  wavelength  of  the  incident  light. 
The  sample  is  mounted  on  a rotational  stage.  The  whole  system 
is  adjusted  by  an  xyz  translation  stage  so  that  the  coupling 
point  remains  practically  stationary  on  the  prism  base  when 
the  stage  is  rotated  through  an  angular  range,  where  the 
coupling  is  possible.  The  linearly  polarized  laser 
beam (parallel  polarization,  He-Ne,  632. 8nm)  was  focussed  into 
the  input  prism  1,  so  that  the  beam  coincides  with  the 
coupling  spot  at  the  prism  base  to  obtain  an  efficient 
coupling.  During  the  measurement,  prisms  1 and  3 are  fixed, 
while  prism  2 is  moved  along  the  guide  in  the  z direction. 
The  output  power  from  prisms  2 and  3 are  detected  by 
photodetectors . 

3.2.10.4  Mode  profile  measurement 

Measurement  of  the  mode  profile  is  of  great  importance  in 
waveguide  characterization.  When  two  waveguides  are  connected 
end  to  end  for  coupling,  the  efficiency  of  power  transfer 
depends  upon  overlap  of  the  fields  in  these  waveguides.  Thus 
a variety  of  coupling  problems  arise,  such  as  fiber  to  fiber, 
laser  to  waveguide,  fiber  to  waveguide,  and  waveguide  to 
waveguide  coupling  etc.  It  is  therefore  necessary  to 
undertake  an  experimental  determination  of  the  mode  profile. 
Moreover,  in  the  cases  such  as  the  performance  of  the 
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directional  coupler,  and  bend-induced  losses,  all  depend  on 
the  mode  profile. 

Two  methods  are  commonly  used  to  carry  out  the  mode 
profile  measurement.  Both  use  an  imaging  lens  of  short  focal 
length  placed  near  the  waveguide (i.e.  near  field  method)  to 
magnify (250x-500x)  the  mode  profile.  For  the  detection 
system,  one  method  uses  a vidicon  camera  placed  at  the  image 
plane  and  an  oscilloscope  to  view  the  video  signal [Kor82 ] , 
while  the  other  way  employs  a suitable  semiconductor 
photodetector  with  a small  aperture  to  scan  along  the  image 
plane  manually [Kei79] . The  latter  is  slower,  while  using  the 
video-oscilloscope  combination  allows  the  entire  profile  to  be 
viewed  continuously  and  the  signal  to  noise  ratio  of  the 
signal  to  be  increased  by  averaging  the  signals  from  several 
successive  frames. 

In  this  study,  the  second  method  was  used  to  measure  the 
mode  profile  of  the  waveguides,  as  shown  in  Figure  3.14.  The 
multimode  planar  waveguide  is  polished  at  the  output  end,  down 
to  0.25um  diamond  paste,  to  avoid  scattering.  The  waveguide 
is  excited  using  a prism  coupler  at  the  input  end.  The  modes 
are  individually  excited  and  the  waveguide  output  is  collected 
by  a high  quality  planar  objective (20x)  lens  and  imaging  the 
waveguide  end  after  approximately  a distance  of  150cm  on  a 
photodetector (Newport  883-UV)  with  a digital  power  meter 
(Newport  Model  815)  . A lOum  slit  was  mounted  at  the  front  end 
of  the  detector.  The  intensity  profile  of  the  image  is 
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Experimental  setup  of  the  near  field  method  for 
measurement. 
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measured  by  scanning  the  slit/detector  combination  using  a 
precision  translator  along  the  image  plane.  The  magnification 
of  the  image  is  determined  by  imaging  an  electron  microscope 
aperture ( lOum  size)  at  the  detector  face  and  scanning  the 
detector.  For  a working  distance  of  150cm  the  magnification 
of  the  near  field  image  is  around  350. 

To  improve  the  sensitivity  of  the  method,  the  beam  is 
mechanically  chopped  at  a frequency  of  50  Hz,  and  the  detector 
output  is  measured  in  conjunction  with  a precision  lock-in 
amplifier  (EG&G  Princeton  Applied  Research  5208  Two  Phase 
Lock-in  Analyzer)  to  detect  the  intensity  profile  along  the 
image  plane. 

3.2.11  Third  Harmonic  Generation (THG)  Measurement 

Due  to  the  increasing  importance  of  understanding  the 
origin  of  the  nonlinear  optical  susceptibility  of  glasses  and 
its  dependence  on  the  constituents  present  in  the  material, 
the  nonlinear  optical  property  of  high  index  glass  films  of 
both  Ge02-Si02  and  Ge02-Ti02  systems  was  investigated.  A Third 
harmonic  generation (THG)  Maker  fringe  measurement  was  used  to 
measure  the  value  of  x(3)  of  films  prepared  by  the  sol-gel 
method.  The  measurement  was  conducted  and  analyzed  in  Dr. 
George  Stegeman's  laboratory  at  the  University  of  Arizona. 

THG  is  an  ideal  technique  for  measurements  of  relative 
small  values  of  xa)  in  films  deposited  on  transparent 
substrates.  The  samples  considered  here  consist  of  a thin 
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film  (~l/um)  deposited  on  a much  thicker  substrate.  The 
expected  THG  signal  from  the  thin  film  is  too  small  to  be 
measured  by  optical  techniques  such  as  Degenerate  Four  Wave 
Mixing  (DFWM) , or  Kerr  rotation  for  incident  directions 
outside  the  film.  The  THG  technique  also  has  the  advantage 
that  it  probes  purely  electronic  contribution  to  the 
nonlinearity  when  the  third  harmonic  frequencies  are  far  from 
resonance.  Therefore,  orientational  and  thermal  effects  as 
well  as  other  dynamic  nonlinearities  derived  from  excitations 
under  resonance  conditions  are  eliminated. 

3.2.11.1  Third  harmonic  generation 

When  an  intense  laser  pulse  is  focused  into  an  isotropic 
medium,  the  total  polarization  of  the  material  can  be  expanded 
into  a power  series  of  the  applied  optical  field.  Because  of 
the  existence  of  center  of  symmetry  in  a macroscopic  sense, 
the  first  nonlinear  process  to  occur  varies  as  the  cube  of  the 
electric  field  [Blo65] 

P(w»)  = Pl  + P nl  = £o  CX(1)E0“  + X(3)  (E0“)3]  (3.2.17) 

The  nonlinear  polarization  cubic  in  the  field  can  lead  to  many 
nonlinear  processes.  Typically  one  nonlinear  process 
dominates  over  the  others  because  phase  matching  occurs,  for 
example  in  DFWM  or  the  optical  Kerr  effect.  In  the  case  where 
the  experimental  geometry  does  not  produce  a nonlinear  phase 
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matched  effect,  Maker  fringes (i.e.  an  interference  fringes) 
can  be  obtained  by  varying  the  optical  path  in  the  nonlinear 
material,  either  by  thickness  or  refractive  index  changes. 
Figure  3.15  shows  the  schematic  of  the  formation  of  THG 
interference  fringes  in  a three-layer  system.  Kajzar  et 
al.[Kaj83]  studied  the  nonlinear  susceptibility  of 
polydiacetylene  films,  deposited  by  the  Langmuir-Blodgett (LB) 
technique  on  both  sides  of  the  substrate,  using  the  THG 
interference  fringes  method (similar  in  shape  to  the  Maker 
fringes  but  of  a different  origin) . These  authors  proposed 
that  the  interference  pattern  results  from  interference 
between  the  harmonic  light  emitted  by  the  film  on  the  forward 
face  with  that  emitted  by  the  film  on  the  backward  face.  The 
dispersion  of  the  refractive  index  of  the  silica  substrate 
introduces  a dephasing  factor  between  harmonic  waves  generated 
at  both  sides  of  the  substrate.  This  dephasing  factor  varies 
as  the  plate  is  rotated  and  the  interaction  length  in  the 
substrate  varied. 

According  to  Kajzar  and  Messier's  analysis [Kaj 86,  Kaj87], 
the  resultant  harmonic  field  for  linearly  polarized  incident 
light  is  equal  to  the  sum  of  contributions  from  consecutive 
nonlinear  media: 

Er3“  = T,E,3“  + T2Es3“  +T3E23u  = Z TjEj3"(rj,  t)  (3.2.18) 

To  calculate  the  radiated  harmonic  field  strengths,  E,3",  E,3", 
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and  E23",  associated  with  each  of  the  three  media,  multiple 
reflections  of  both  the  fundamental  and  harmonic  fields  are 
ignored.  This  can  be  justified  for  angles  of  incidence  away 
from  the  normal.  Each  harmonic  field  at  the  output  of  a 
nonlinear  plane  parallel  slab  of  thickness  d is  given 
by[Kaj86] : 


£3“  = 4tt( 


) ATe‘  *H[e‘  - 1 ] 


(3.2.19) 


where  P^l/4  xa)  (E0“)3  is  the  nonlinear  polarization  inside  the 
slab,  T(T,,  Ts,  T2)  is  a total  transmission  factor,  A(Alf  A,,  A2) 
is  a coupling  factor  determined  by  boundary  conditions,  E0"  is 
the  fundamental  laser  beam  electric  field,  Ae  is  the 
dielectric  constant  dispersion ( As  = (n“) 2- (n3") 2,  where  n"  and 
n3“  are  the  refractive  indices  at  the  fundamental  and  harmonic 
frequencies  respectively)  , and  At p(  A\pf  and  A\p,)  , is  the  phase 
mismatch  between  the  fundamental (F)  and  harmonic (H)  waves 
inside  the  slab: 


A\J/  = \p  -\p  = 67rd(n"cos0"-n3“cos03“)  (3.2.20) 


where  0",  and  03“  are  the  propagation  angles  of  the  fundamental 
and  harmonic  waves,  respectively,  inside  the  slab,  and  xw(- 
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3oj,w,gj,w)  is  the  third-order  susceptibility  tensor.  All  these 
quantities  will  be  complex  in  an  absorbing  medium  for  both  the 
fundamental  and  harmonic  waves. 

Further,  the  resulting  harmonic  light  intensity  can  be 
expressed  as  follows[Kaj85] : 


r3U=  ce  (3 03) 
87 T 


| £ T,  £/"<!•,,  t)P 


(3.2.21) 


therefore, 


-3<j  _ 


64nJ 


,0) 


) +T2(e‘^-1)  + p T3  e ‘ ( ^ * 0)  ( e ‘ ^ - 1 ) ] | (X^) 3 


(3.2.22) 


which  leads  to  a complex  interference  spectrum  as  a function 
of  the  angle  of  incidence  when  the  substrate  is  rotated  along 
an  axis  perpendicular  to  the  incident  light.  In  equation 


(3.2.22)  : 


(X(3)/A  e)f 
(X^/Ae), 


’A  ~ iAf  + ’Aw  + lA// 
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T12,3  are  transmission  factors  arising  from  the  boundary 
conditions,  0 is  the  phase  of  film  susceptibility  provided 
that  of  substrate  is  equal  to  zero  (in  the  opposite  case,  it 
is  the  difference) , and  I0“  is  the  incident  light  intensity. 

Because  of  the  small  film  thickness ( 1<<1C,  where  1c=tt/|AJc| 
=Xu/6  (n3“-n") ) , the  phase  of  x/3)  relative  to  that  of  silica  does 
not  significantly  change  the  harmonic  light  intensities  but 
only  the  positions  of  the  interference  maxima  and  minima. 

Equation  (3.2.22)  was  used  to  fit  the  experimentally 
recorded  interference  fringes.  From  the  shift  and  the  offset 
of  the  interference  fringes  between  the  clean  silica  substrate 
and  the  sample(i.e.  films  with  substrate),  both  the  magnitude 
and  the  phase  of  xp)  could  be  obtained [Tor9 la] . 

3.2.11.2  THG  interferometry  measurement 

The  THG  experimental  setup  is  shown  in  Figure  3.16.  The 
fundamental  beam  of  light  obtained  from  a Q-switched  Nd : YAG 
laser  operating  at  1.064/xm  is  frequency  doubled  to  0.532/un  by 
a phase-matched  KDP  single  crystal.  The  frequency  of  the 
outgoing  laser  beam  is  subsequently  shifted  in  a high 
pressure (~56bars)  hydrogen  cell  by  4155cm’1  or  its  multiple 
through  a stimulated  Raman  scattering  process.  The  dispersive 
prism  permits  a separation  between  other  Raman  frequencies  and 
the  fundamental  one.  A system  of  filters  allows  elimination 
of  the  parasitic  light.  In  order  to  obtain  the  Maker  fringes, 
the  vertically  polarized  near  infrared  radiation  generated  by 
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Figure  3.16  Experimental  setup  for  THG  interferometry 
measurement. 
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the  H2  shifter(i.e.  953nm,  1064nm,  and  1579nm)  was  focused 
with  a 20cm  long  focal  length  lens  at  the  center  of  a 200 
mTorr  vacuum  chamber.  The  sample  was  rotated  with  an  IBM  PC- 
controlled  Oriel  stepper-motor  rotation  stage.  The  pump  beam 
was  filtered  with  K65  Schott  filters  and  a 1/4  m Jarrell-Ash 
spectrometer.  The  signal  was  detected  with  a Hamamatsu  R-928 
photomultiplier  tube,  preamplified,  and  sent  to  a SRS  Boxcar 
integrator.  The  data  was  taken  by  typically  averaging  30 
pulses  per  data  point  and  scanning  every  half  degree  between  - 
30  and  30  degrees  measured  relative  to  the  normal  direction  to 
the  film.  The  DC  output  signal  of  the  Boxcar  was  digitized 
and  recorded  by  the  IBM  PC  as  a function  of  angle. 

Two  sets  of  Maker  fringes  were  typically  recorded  for 
each  sample  at  every  wavelength.  After  fringes  of  the  clean 
Si02  substrate  were  recorded,  the  sol-gel  film(e.g.  30Ge02- 
70Si02  , 50Ge02-50Si02,  3 0Ge02-7  0Ti02,  and  50Ge02-50Ti02)  was 
deposited  on  both  sides  of  the  substrate,  and  the  fringes  of 
the  sample  were  taken.  The  measurement  was  performed  in 
vacuum  to  eliminate  the  influence  of  the  harmonic  field  that 
can  also  be  generated  from  surrounding  air  during  the 
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CHAPTER  IV 


Ge02-Si02  FILMS  BY  SOL-GEL  METHOD 
4.1  Introduction 

The  sol-gel  process  offers  a new  approach  to  the 
synthesis  of  oxide  coatings.  In  such  a process,  molecular 
precursors  are  transformed  into  an  oxide  network  via  inorganic 
polymerization  reactions.  The  molecular  precursor  is  usually 
a metal  alkoxide,  M(OR)4,  and  polymerization  occurs  via 
hydrolysis  and  condensation.  The  overall  reaction  can  be 
described  as  follows: 

hydrolysis:  M(OR)4  + x H20  ==>  M(OH)x  + x ROH 
condensation:  M(OH)x  ==>  MO^  + (x/2)  H20 

When  working  with  a multicomponent  sol,  the  reaction 
rates  of  various  components  have  to  be  comparable  to  prevent 
preferential  reaction  of  one  component  which  could  lead  to 
premature  precipitation  of  one  component  from  the  homogeneous 
sol.  In  order  to  obtain  a homogeneous  Ge02-Si02  sol,  the 
germania  sol  was  prepared  using  germanium  alkoxide  which  was 
modified  by  higher  alkyl  groups  using  an  alcohol-interchange 
reaction  as  described  in  section  3.1.1.  It  was  then  mixed 
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with  the  silica  sol  which  was  prehydrolyzed  with  HC1/H20, 
according  to  the  procedures  given  in  Figure  3.2. 

This  chapter  will  focus  on  fabrication  of  the  Ge02-Si02 
films,  which  may  be  used  as  planar  waveguides,  using  the  sol- 
gel  method.  Results  related  to  the  preparation  of  films  as 
well  as  to  various  film  characterizations  will  be  presented. 

4.2  Results  and  Discussions 

4.2.1  Synthesis  of  Alkoxide  Precursors 

(1)  . Ammonia  method:  preparation  of  germanium  ethoxide [Bra56] . 

Anhydrous  ammonia  gas  was  slowly  bubbled  into  a solution 
of  germanium-tetrachloride,  ethanol  and  hexane  as  described  in 
section  3.1.1.  There  was  immediate  formation  of  ammonium 
chloride,  in  both  vapor  phase  and  in  solution,  and  an 
appreciable  amount  of  heat  was  liberated.  The  supply  of  gas 
was  discontinued  and  when  the  reaction  mixture  cooled  to  room 
temperature.  After  filtration  and  distillation,  a colorless 
liquid  of  germanium  ethoxide  was  obtained.  The  identification 
of  the  compound  was  performed  by  using  FTIR  and  a ZnSe  liquid 
cell,  the  spectra  of  germanium  ethoxide  and  pure  ethanol  are 
shown  in  Figure  4.1  and  Figure  4.2,  respectively.  Most  of  the 
absorption  peaks  in  ethanol  match  those  in  germanium  ethoxide 
except  the  broad  band  that  occurs  at  -3  3 60  cm'1  corresponding 
to  a vibration  mode  related  to  -OH  groups,  which  does  not 
exist  in  germanium  ethoxide,  Ge(OC2H5)4. 

(2) .  The  alcohol-interchange  method:  preparation  of  germanium 
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Figure  4 . 2 
ethanol . 
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compounds  with  higher  alkyl  groups. 

The  alcohol  interchange  method  has  been  used  by  Bradley 
et  al.[Bra53]  to  prepare  various  germanium  esters.  This 
technique  was  utilized  to  modify  the  starting  germanium 
ethoxide  to  replace  the  smaller  groups  by  higher  alkyl  groups. 
The  experimental  procedure  was  described  in  detail  in  section 
3.1.1.  The  replacement  of  the  ethyl  groups  with  higher  alkyl 
groups  is  expected  to  increase  the  stability  of  the  germanium 
alkoxide  against  hydrolysis/condensation  reaction  due  to 
steric  effects.  The  success  of  this  modification  was  analyzed 
by  using  FTIR  with  a ZnSe  liquid  cell.  Figure  4.3  and  Figure 
4 . 4 show  the  transmittance  spectra  of  germanium  ethoxide 
modified  by  2-methoxyethanol  and  pure  2-methoxyethanol, 
respectively.  When  compared  with  the  spectrum  of  germanium 
ethoxide,  the  absorption  bands  corresponding  to  ethyl  groups 
were  replaced  by  the  features  related  to  2-methoxy-ethyl 
groups,  indicating  the  high  efficiency  of  the  alcohol 
interchange  reaction.  In  the  present  study,  propanol,  n- 
butanol,  2-methoxy-ethanol , and  2 -ethoxy-ethanol  have  been 
used  as  modifying  alcohols. 

4.2.2  Sol  Preparation 

The  Ge02-Si02  sols  were  prepared  by  mixing  a prehydrolyzed 
TEOS,  and  germania  sols  together.  The  germania  sols  have  been 
prepared  using  various  modified  germanium  alkoxides  as  shown 
in  Table  4.1.  Because  of  the  reduction  of  the  hydrolysis/ 
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Figure  4 . 4 
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Table  4.1  The  Effect  of  Alcohol  Interchange  to  the  Quality  of 
Germania  Sols. 


STARTING 

COMPOUND 

MODIFYING 

ALCOHOL 

SOLVENT 

CATALYST 

WATER 

REMARK 

60(002^)4 

No 

Ethanol 

No 

No 

Cloudy 

GefOCgHs^ 

Propanol 

Ethanol 

No 

No 

Cloar 

Unstablo 

60(002*5)4 

60(002*5)4 

Butanol 

2-methoxy- 

ethanol 

Ethanol 

Ethanol 

No 

HCI 

No 

Yos 

Cloar 

Stable 

Clear 

Stable 

Ge(OC2H5)4 

2-ethoxy- 

ethanol 

Ethanol 

HCI 

Yas 

Clear 

Stable 
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condensation  reaction  rates  by  using  alkoxides  with  larger 
alkyl  groups,  i.e.  steric  hindrance  effect,  clear  and  stable 
germania  sols  can  be  obtained  by  using  germanium  compounds 
with  higher  alkyl  groups.  The  achievement  of  a clear  sol  that 
remains  stable  for  a long  enough  time  is  important  for  a 
multiple  dip-coating  process.  This  ensures  the  homogeneous 
composition  of  the  sol,  which  in  turn  allows  depositing  films 
of  uniform  composition  without  degrading  the  film  quality.  In 
the  current  study,  germanium  2-methoxy-ethoxide  was  used  as 
the  germanium  precursor  to  prepare  germania  sols  for  making 
Ge02-Si02  waveguides  of  different  compositions.  The 
preparation  of  the  multicomponent  sol  was  conducted  according 
to  the  procedure  given  in  Figure  3 . 2 as  described  in  section 
3.1.1. 

4.2.3  Film  Formation 

One  of  the  major  limitations  of  the  sol-gel  coating 
process  for  practical  application  is  the  small  thickness  of 
each  coating,  which  makes  multiple  coatings  necessary  for 
building  up  the  required  thickness  for  specific  applications. 
Therefore,  doing  multiple  coatings  without  degrading  the  film 
quality  is  the  major  challenge  of  the  sol— gel  coating  process. 
In  the  dip-coating  process,  the  withdrawal  speed  of  substrate 
and  viscosity  of  sol  are  usually  low.  The  deposited  film 
thickness (h)  can  be  expressed,  according  to  the  relationship 
derived  by  Landau  and  Levich[Lan42 ] as  follows: 
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0-94(7^) 2/3  (4.2.1) 

7 iv(Pgr) 1/2 


where  ij=sol  viscosity,  U0=withdrawal  speed,  7LV=liquid-vapor 
surface  tension,  p=density  of  the  sol,  and  g=acceleration  of 
gravity.  The  thickness  of  the  coating,  therefore,  can  be 
increased  by  increasing  the  sol  viscosity  and/or  pulling 
speed.  In  the  present  study,  the  maximum  pulling  speed  of 
the  dip-coating  system  is  usually  used  (i.e.  3.8  cm/min.),  and 
several  organic  additives,  such  as  polyvinylpyrrolidone (PVP  k- 
30,  GAF) , cp2(Daxad  cp2 , W.R.  Grace  Co.)  and  hydroxypropyl- 
cellulose (HPC,  Aldrich) , were  used  in  an  attempt  to  adjust  the 
viscosity  of  the  sol.  Among  them,  only  HPC  was  found  to  be 
compatible  with  the  Ge02-Si02  sols,  thus  it  can  be  added  to  the 
sol  without  causing  precipitation  or  degradation  of  the  sol. 
The  amount  of  HPC  used  is  usually  less  than  2 wt%.  A thicker 
film  with  improved  uniformity  was  obtained  by  using  HPC  as  an 
independent  viscosity  adjusting  agent.  Good  quality  films 
made  with  20-72  dippings  were  obtained  by  this  method.  The 
coatings  were  fired  after  every  10  to  12  dippings  depending  on 
the  sol  composition. 

Table  4.2  shows  the  results  of  composition  analysis  by 
using  ICP  solution  analysis  of  the  gels  that  obtained  from  the 
sols  used  for  depositing  thin  films.  The  compositions 
measured  are  close  to  the  target  compositions  except  for 
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Table  4.2  ICP  Solution  Analysis  of  Various  Si02-Ge02  Gels  and 
Trace  Impurities  in  These  Gels. 


SAMPLE 

COMPOSITION(mole%) 

SiO  2 

GeO  2 

10GS 

90.5 

9.5 

30GS 

71.7 

28.3 

50GS 

57.4 

42.6 

75GS 

33.0 

67.0 

(b) 


SAMPLE 

TRACE  IMPURITIES 

(ppm/ml) 

Na 

Mg 

Al 

Fe 

Cu 

10GS 

<0.7 

<0.01 

<1.8 

<0.2 

<0.1 

30GS 

M 

M 

H 

N 

M 

50GS 

M 

M 

H 

N 

M 

75GS 

II 

M 

II 

II 

II 
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higher  germania  contents,  e.g.  50GS  and  75GS.  The  major 
reason  for  the  deviation  is  believed  to  be  due  to  the  weighing 
error  during  sol  preparation,  because  very  small  amounts  of 
alkoxides  were  used  for  mixing  the  sol.  Another  possible 
source  of  error  is  the  existence  of  a small  amount  of  residual 
solvent  in  the  germanium  precursor,  which  was  not  removed 
completely  during  vacuum  distillation. 

Several  possible  common  contaminants (Na,  Mg,  Al,  Fe,  and 
Cu)  that  exist  in  silicate  glass  were  analyzed.  All  of  the 
elements  analyzed  were  below  the  detection  limit  of  the  ICP 
system  used.  The  value  reported  in  the  Table  4.2(b)  is  100 
times  of  the  detection  limit,  therefore,  these  impurities  are 
in  the  ppb  level,  if  present. 

The  composition  of  the  films  was  also  analyzed  by  XPS 
which  confirms  the  results  obtained  by  ICP  analysis.  Fig  4.5 
shows  the  XPS  spectrum  of  the  3 0GS  Ge02-Si02  thin  film.  No 
extra  peak  except  Si,  Ge,  0,  and  trace  amount  of  C,  which 
commonly  exist  on  material  surfaces,  was  observed. 

4.2.4  Thermal  Analysis 

Thermal  behavior  of  the  gels,  including  the  behavior 
during  organic  burn-out,  was  studied  using  TGA/DTA  analysis. 
The  result  of  TGA  analysis  on  30GS  is  shown  in  Figure  4.6(a). 
The  TGA  curve  clearly  shows  two  distinct  stages  of  weight  loss 
from  the  sample  while  heating  from  room  temperature  to  1175° 
C.  About  8%  weight  loss  occurs  upon  heating  to  150°  C.  The 
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XPS  composition  analysis  of  30GS  film 
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Figure  4 . 6 


(a)  TGA  and  (b)  DTA  results  for  3 0GS  gel 


prepared  with  1.5  moles  water  per  mole  of  alkoxide. 
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DTA  curve(Figure  4.6(b))  shows  a broad  and  small  endotherm  in 
the  same  temperature  range,  suggesting  that  the  weight  loss  in 
this  temperature  range  is  primarily  due  to  the  elimination  of 
residual  solvent (ethanol)  , physically  adsorbed  water  and 
possibly  chemically  bound  water (i.e.  dehydroxylation) . At 
higher  temperatures,  a strong  exotherm  is  observed  which  peaks 
around  335°  C.  The  TGA  curve  shows  a corresponding  weight 
loss,  i.e.  18.5%,  in  the  same  temperature  range  and  it  is 
believed  that  this  peak  is  primarily  associated  with 
combustion  of  the  organic  binder (i.e.  HPC)  as  well  as 
carbonaceous  material  that  may  be  present  due  to  incomplete 
hydrolysis  of  the  alkoxide.  This  interpretation  is  consistent 
with  the  observation  that  powder (or  films)  removed  from  the 
furnace  in  this  temperature  range,  i.e.  200  to  600°  C,  has  a 
slight  brown  color  characteristic  of  material  with  carbon 
chars.  This  organic  burn-out  process  should  be  completed  when 
the  sample  is  heated  up  to  600°  C.  The  TGA  curve  shows 
basically  a flat  line  in  the  temperature  range  above  600°  C, 
and  the  sample  removed  from  furnace  after  being  heated  above 
600°  C has  a clear  appearance.  Similar  thermal  behavior  was 
observed  for  10GS  gel  as  shown  in  Figure  4.7. 

4.2.5  Consolidation  of  the  Film 

In  order  to  achieve  a coherent  film  without  cracking 
during  organic  removal  and  densif ication  processes,  the  films 
were  heat  treated  following  a program  as  shown  in  Figure  4.8. 
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Figure  4.7  (a)  TGA  and  (b)  DTA  results  for  10GS  gel 

prepared  with  2 moles  of  water  per  mole  of  alkoxide. 
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Figure  4.8  Typical  temperature  profile  for  heat-treating 

30GS  sol-gel  film. 
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The  upper  densif ication  temperature  ranges  from  700°C  to  1110°C 
depending  on  the  composition  of  the  film.  Figure  4.9  shows 
the  DTA  curve  of  50GS  gel,  heat  treated  using  the  same 
temperature  profile  as  shown  in  Figure  4.8.  Following  the 
small  endotherm  below  100°C,  the  combustion  of  organic 
additive  occurs  in  approximately  three  distinct  temperature 
ranges,  between  around  180°C  and  600°C,  with  a total  weight 
loss  of  about  45%. 

Since  the  porosity  and  the  refractive  index  of  the  porous 
film  can  be  related[Yol80]  , the  consolidation  of  the  films 
with  various  germania  contents  was  monitored  by  measuring  the 
change  of  refractive  index,  after  heat  treatment  in  oxygen  at 
various  temperatures  for  10  minutes,  using  a heating  rate  of 
4°/min.  The  consolidation  curves  for  all  Ge02-Si02  films 
studied,  as  shown  in  Figure  4.10,  show  a decrease  in  the 
index  of  refraction  in  the  temperature  range  between  200°C  and 
600°C , following  an  initial  increase  by  heating  to  200°C.  This 
is  consistent  with  the  thermal  analysis  data.  The  decrease  in 
index  in  this  temperature  range  is  due  to  the  loss  of  organic 
additives,  before  the  beginning  of  sintering.  The  index 
increases  slowly  after  600°C  and  peaks  at  different 
temperatures  for  different  film  compositions.  The  film  with 
higher  germania  content  appears  to  density  at  a lower 
temperature.  The  reason  for  the  subsequent  small  drop  in 
index,  in  films  of  higher  Ge02  content,  at  high  temperatures 
(>800°C)  is  not  well  understood.  However,  according  to  FTIR 
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Figure  4 . 9 


DTA  result  for  50GS  gel 
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studies,  as  shown  in  Figure  4.11,  it  may  be  related  to  the 
loss  of  a large  amount  of  hydroxyl  groups (i.e.  dehydroxylation 
process)  which  occurs  between  800°  and  900°C  in  the  30GS  gel. 

If  the  film  were  heat  treated  at  excessively  high 
temperatures  or  held  too  long  at  high  temperatures, 
crystallization  would  take  place.  Figure  4.12(a)  shows  the  X- 
ray  diffraction  spectrum  of  the  50GS  waveguiding  film  (treated 
between  650°C~850°C  for  6 hrs.)  , which  appears  to  be  amorphous. 
The  Ge02  hexagonal  phase  crystallizes,  as  shown  in  Figure 
4.12(b),  following  a heat-treatment  in  the  same  temperature 
range  for  19  hours.  The  observation  of  the  formation  of  the 
hexagonal  phase  agrees  with  Mukerjee's  work  [Muk84,  Muk86]  on 
Ge02-Si02  monolithic  gel  glasses,  in  which  he  observed  the 
appearance  of  hexagonal  Ge02  at  high  Ge02  concentrations. 

Figure  4 . 13  shows  the  variation  in  refractive  index  of 
the  film  as  a function  of  germania  content.  This  large  change 
in  refractive  index  allows  the  preparation  of  waveguides  with 
various  numerical  apertures,  which  are  compatible  with  single 
mode  optical  fibers.  The  experimentally  measured  values  are 
in  good  agreement  with  the  values  calculated  by  using  the 
Lorentz-Lorenz  equation [Huan78 ] , as  shown  by  the  line  in 
Figure  4.13,  to  be  described  in  the  following  sections. 

In  general,  the  density  of  glass  depends  upon  the  type  of 
ion  present,  the  way  these  ions  are  arranged  and  the  thermal 
history  of  the  glass [Hug43 ] . In  the  case  of  silica  glass,  the 
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Figure  4.11  FTIR  transmittance  spectra  showing 
variation  of  hydroxyl  content  in  the  3 0GS  gel  with 
treatment  temperature. 


the 

heat 
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Figure  4.12  XRD  spectra  of  50GS  film  heat-treated  between 
650°~850°C  for  (a)  6 hr.  and  (b)  19  hr.,  compared  to  that  of 
Ge02  hexagonal  crystal. 
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Figure  4.13  Variation  of  refractive  index  of  Ge02-Si02  films 
with  Ge02  content.  + : Ellipsometry  data,  solid  line: calculated 
values  using  the  Lorentz-Lorenz  model. 
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oxygen  and  silicon  atoms  are  linked  together  by  ionic  forces, 
with  the  silicon  ion  electrostatically  balanced  by  four 
bridging  oxygen  ions  forming  a tetrahedral  structural 
unit [Bor 69].  Germania  glass  has  the  same  atomic 
arrangement [Rie68 ] . To  a first  approximation,  it  is  reasonable 
to  assume  that  for  the  entire  composition  range  of  binary 
germania-silica  glasses,  the  Si04  and  Ge04  units  within  the 
glass  retain  their  tetrahedral  structure  and  form  a solid 
solution,  and  the  effective  volume  of  the  ions  contained  in 
the  glasses  remains  invariant.  Consequently  the  specific 
volume  of  germania-silica  glass  is  additive  and  thus  the 
density  can  be  written  as[Huan78]: 


, M.  M, 

IV  N 2'  N ' 


(4.2.2) 


where  M,— molecular  mass  of  Si02,  M2=molecular  mass  of  Ge02, 
N=Avogadro ' s number , and  Nj  and  N2  are  the  numbers  of  molecules 
per  cubic  centimeter  of  Si02  and  Ge02,  respectively.  The 
concept  of  structural  compactness  leads  to  the  following 
expression: 


, N.M . NJi, 

(-L2)  +(^Li)  =i 


(4.2.3) 
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where  and  d2  are  the  densities  of  Si02  and  Ge02, 
respectively.  For  a glass  composition  that  is  expressed  in 
terms  of  mole  fraction  of  germania,  x. 


W2 

N}+N2 


(4.2.4) 


From  eguation  (4.2.2)  (4.2.3)  and  (4.2.4)  the  density  of  the 

Ge02-Si02  glass  can  be  written  as 


[#2*+  (l-x)W,] 
[W,(  l-x)/d,]  +xM2/d2 


(4.2.5) 


The  refractive  index  n is  of  prime  importance,  when  the 
glass  is  used  as  a dielectric  waveguide.  The  dielectric 
constant  of  glass  contains  contributions  from  four 
polarization  processes:  (l)  space  charge (or  interfacial 
polarization)  (2)  orientation  polarization  (3)  atomic 
polarization  and  (4)  electronic  polarization.  The  Clausius- 
Mosotti  equation,  which  describes  the  relationship  between  the 
dielectric  constant,  e,  of  a material  and  electronic  ae, 
atomic  a,,  orientational  ad,  and  interfacial  a;,  polarizability 
can  be  expressed  as 
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(±L)  (N'Ct'  + N'd'+N/Xt+Nfr) 


(4.2.6) 


where  Ne,  N,,  N,,  and  Nj  are  the  number  per  unit  volume  of 
polarizable  electronic,  atomic,  orientational  and  interfacial 
species,  respectively. 

In  the  optical  frequency  range,  the  dielectric  constant 
arises  entirely  from  the  electronic  polarizability,  because 
the  atomic,  orientational  and  interfacial  polarizations  are 
unable  to  respond  to  the  frequency  of  the  external  field  due 
to  the  inertia  of  the  atoms  and  molecules [Pye72 ] . 
Furthermore,  utilizing  Maxwell's  relation  e=n2,  equation 
(4.2.6)  becomes 


(n^-D  ,(«»)„  (4.2.7) 
(n2  + 2)  3 ' ' * 

Equation  (4.2.7)  is  the  Lorentz-Lorenz  equation  and  is  valid 
only  for  values  of  the  n and  molecular  electronic 
polarizabilities  at  the  same  optical  frequency.  By  applying 
the  additive  model  for  binary  Ge02-Si02  glasses,  equation 
(4.2.7)  can  be  expressed [Con67 ] as 


(n2- 1) 

(n2  + 2 ) 


4 IT 

= ) («.f  Nl+ae2N2) 


(4.2.8) 
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where  ael  and  ote2  are,  respectively,  the  molecular  electronic 
polarizability  of  v-si02  and  v-Ge02;  N,  and  N2  are  the  numbers 
of  molecules  per  cubic  centimeter  of  Si02  and  Ge02, 
respectively.  The  refractive  index  of  the  binary  glass  system 
can  be  expressed  as  [Huan78] 


2c  + 1 (4.2.9) 

1 -c 


where 


c-iE[ 

3 L 


N(l-x)aeI 
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) (^i)  (i-x)  + (!!**) 
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where  N is  Avogadro's  number  (N=6 . 02xl023/mole)  . The  refractive 
indices  for  binary  Ge02— Si02  glass  films  were  calculated  from 
equation  (4.2.9),  using  d,  = 2.202  g/cm3  , d2  = 3.604  g/cm3, 
^ei=^  * 952x10 34  cm3,  ae2=3 . 962xl0'24  cm3,  Mi=60.09  g/mole,  and 
M2=104 . 59  g/mole [Huan78 ] . 

The  formation  of  hexagonal  Ge02  structure  from  the  binary 
Ge02-Si02  glass,  especially  for  the  ones  of  higher  Ge02 
contents,  suggests  the  possibility  of  coexistence  of  both 
four-coordinated  and  six-coordinated  Ge  species  in  the 
material.  The  possible  existence  of  six-coordinated  Ge 
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species  indicates  an  increased  packing  efficiency  of  0 about 
Ge  to  a coordination  of  6 which  in  turn  would  be  expected  to 
raise  the  index.  This  effect  may  be  the  reason  for  the  data 
of  higher  Ge02-content  samples  being  higher  than  the  above 
Equation. 

4.2.6  Microstructure  of  the  Film 

The  surface  morphology  of  the  waveguiding  film(e.g.  50GS) 
was  examined  under  SEM.  Figure  4.14  is  the  SEM  micrograph  of 
the  50GS  film  surface.  Two  types  of  'defects'  were  observed: 
(1)  There  are  many  humps,  as  shown  in  Figure  4.14(a)  on  the 
film,  which  may  be  caused  by  the  fact  that  the  film  was 
deposited  in  an  open  environment  where  dust  or  other  extrinsic 
particulates  may  be  introduced  into  the  sol  or  the  film  during 
the  multiple  coating  process.  Most  of  these  humps  are  around 
0.1  /xm  size  range.  (2)  A few  pinholes,  as  shown  in  Figure 
4.14(b),  were  also  observed  in  the  film,  although  the 
concentration  is  very  low.  These  defects  of  the  film  can 
cause  scattering  of  light  when  the  film  is  used  for 
waveguiding. 

The  microstructure  of  the  sol-gel  prepared  thin  film  were 
investigated  using  TEM.  The  film  consists  of  very  fine  parti- 
cles in  the  size  range  below  lOnm  as  shown  in  Figure  4.15.  in 
some  areas,  very  small  crystallites  were  observed  which  show 
up  as  dark  spots  in  the  bright  field  image,  Figure  4.16(a). 
These  small  crystalline  clusters (<10nm)  are  better  shown  in  a 
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Figure  4.14  SEM  micrographs  of  50GS  waveguiding 
showing  "defects"  in  film  formation,  (a)  surface  humps, 
(b)  pinhole. 
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Figure  4.15  TEM  microgragh  of  50GS  waveguiding 
showing  the  film  consists  of  very  fine  structures  at  the 
of  -lOnm. 


film, 

scale 
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Figure  4.16 
image,  (b) 
pattern. 
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dark  field  image  of  the  (200)  ring  in  SAD  pattern,  Figure 
4.16(b),  in  which  case  these  crystallites  appear  as  bright 
spots  due  to  intense  diffraction  which  occurs  at  the  small 
crystallites.  Selected  area  electron  diffraction (SAD)  was 
used  to  analyze  these  crystalline  phases.  The  SAD  pattern  was 
given  in  Figure  4.16(c),  and  shows  several  diffuse  rings, 
indicating  the  formation  of  a crystalline  phase  in  the  50GS 
film,  although  it  was  not  detectable  by  X-ray  diffraction 
studies  as  shown  Figure  4.12(a)  . The  analysis  of  the  SAD  ring 
pattern  is  summarized  in  Table  4.3.  These  crystallites  are 
identified  to  be  the  Ge02  hexagonal  phase [JCP88 ] . 

4.2.7  Vibrational  Spectroscopy 

Infrared  and  Raman  spectroscopy  were  used  to  investigate 
the  structure  of  Ge02-Si02  gel  glasses  as  well  as  thin  films. 
Figure  4 . 17  shows  the  FTIR  transmittance  spectra  of  Ge02-Si02 
gels  with  various  Ge02  contents,  namely  10GS,  30GS,  and  50GS, 
after  being  heat  treated  to  900°C  in  02.  The  peaks  around 
580cm1,  which  were  identified[Bor69]  as  the  Ge-O-Ge  stretching 
modes  increase  with  increasing  Ge02  content.  The  peaks  around 
670  cm'1  and  1020cm1,  which  were  identified  by  Borrelli  as  Si- 
O-Ge  bond  stretching  modes,  also  increase  with  increasing  Ge02 
content,  indicating  the  formation  of  a homogeneous  solid 
solution  of  Ge02-Si02  glass  with  Si-O-Si,  Si-O-Ge,  and  Ge-O-Ge 
linkages  that  interconnect  into  a continuous  3-dimensional 
network,  instead  of  segregation  of  the  Si02  and  Ge02  phases. 
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Table  4.3  The  Analysis  of  SAD  Pattern  of  50GS  waveguide 


ring 

diameter (r) 
(cm) 

d spacing 
(A) 

Ge02-hex. 
(h  k 1) 

1 

0.975 

2.159 

2 0 0 

2 

1.125 

1.870 

112 

3 

1.590 

1.324 

10  4 

4 

1.860 

1.133 

2 2 2 

135 


o 


Figure  4.17  FTIR  transmittance  spectra  of  Ge02-Si02  sol-gel 
glass  after  being  fired  to  900°C. 
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Another  interesting  feature  of  these  spectra  is  the  reduction 
of  the  intensity  of  the  800cm'1  band  due  to  Si-O-Si  bond 
stretching  which  appears  more  pronounced  than  the  reduction  of 
both  the  1180cm'1  and  1090cm'1  bands,  which  are  assigned  to  the 
Si-O-Si  LO  anti -symmetric  stretching  mode  and  TO  anti- 
symmetric stretching  mode  respectively.  The  implications  of 
this  finding  are  not  fully  understood.  In  the  present  study, 
the  Ge02-Si02  gel  with  -75%  Ge02  shows  some  crystallization  of 
the  Ge02  hexagonal  phase  after  it  is  heat  treated  to  900°C  in 
02 . For  pure  Ge02  gel,  the  formation  of  Ge02  crystalline 
clusters  was  found  in  as-prepared  samples.  The  pure  germania 
gel  retains  a large  amount  of  hydroxyl  groups  even  after  being 
heat  treated  to  1000°C  in  02. 

FTIR  transmittance  spectra  of  Ge02-Si02  waveguides  of 
various  Ge02  content  were  obtained  using  a ATR  stage.  The 
results  were  shown  in  Figure  4.18.  Comparing  to  a pure  Si02 
film,  the  germania-doped  film  has  absorption  bands  at  around 
670cm'1  and  1020cm1  corresponding  to  different  Si-O-Ge 
stretching  modes,  suggesting  the  formation  of  a homogeneous 
glass  network  of  interconnected  Si04  and  Ge04  tetrahedra 
instead  of  segregation  of  the  Si02  and  Ge02  phases  in  the 
coatings.  A limitation  to  the  measurements  using  the  ATR 
stage  results  from  the  film  being  so  thin  that  interference 
from  the  Si02  substrate  is  possible.  Moreover,  the  films, 
which  were  used  as  waveguides,  have  different  thicknesses  for 
different  Ge02  content.  Therefore,  only  qualitative 


TRANSMITTANCE 


137 


WAVENUMBER 


Figure  4.18 
thin  films. 


FTIR  transmittance  spectra  of  Ge02-Si02  sol-gel 
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information  may  be  obtained  from  these  measurements. 

Raman  spectroscopy  was  also  used  to  investigate  the 
distribution  of  tetrahedra  in  multicomponent  network 
structures  in  the  Ge02-Si02  glass  films.  The  results  are  shown 
in  Figure  4.19.  When  compared  to  the  spectrum  of  pure  Si02 
film,  several  differences  are  observed  in  Ge02-doped  silicate 
films : 

(1) .  The  intensity  of  sharp  bands  of  the  Si02  glass  spectrum 

at  490  cm1 (D1  band)  and  606  cm1(D2  band)  decreases  with 
increasing  Ge02  content. 

(2)  . The  strong,  broad  and  asymmetric  Si02  glass  envelope  with 

the  maximum  at  -4  3 0cm'1  appears  to  narrow  down  in  band 
width  and  develop  a weak  shoulder  at  -580  cm'1;  the  band 
maximum  shifts  to  -418  cm'1. 

(3)  . The  -800  cm1  envelope,  which  is  a doublet  composed  of 

bands  at  797cm'1  (TO)  and  -830cm'1  (LO)  due  to  the  Si-O-Si 
stretching  mode,  decreases  in  intensity  with  increasing 
Ge02  addition. 

(4) .  A small  weak  band  appears  at  -675  cm'1  due  to  the  Si-O-Ge 

stretching  mode,  which  increases  with  increasing  Ge02 
content. 

The  identification  of  Raman  spectra  in  the  current  study 
was  based  on  the  work  of  Sharma  et  al.[Sha84]  and  Shibata  et 
al.  (Shi81)  . Sharma  et  al.  studied  Ge02-Si02  bulk  glasses  with 
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Figure  4.19  Raman  spectra  of  Ge02-Si02  thin  films  prepared 
by  sol-gel  process. 
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different  Ge02  contents,  prepared  by  a melting  technique, 
while  Shibata  et  al.  made  their  samples  using  a VAD  technique. 

The  disappearance  of  D1  and  D2  silica  defect  bands  in 
Ge02-Si02  Raman  spectra  has  been  observed  by  many  researchers 
[Bat7  2 , Sch7  3 , Spr81,  Chm89],  The  490cm'1  band  in  the  spectrum 
of  Si02  glass  has  been  assigned  to  the  presence  of  four- 
membered  rings  of  Si04  tetrahedra  stabilized  by  the  defect 
centers  responsible  for  the  606  cm'1  band[Sha81].  The  absence 
of  the  490  cm1  band  in  the  spectrum  of  Ge02-Si02  glass  films 
indicates  that  the  addition  of  Ge02  to  Si02  makes  the  four- 
membered  rings  of  Si02  tetrahedra  unstable.  The  instability 
of  the  four-membered  ring  structures  may  be  due  to  elimination 
of  defect  centers  in  the  Si02  glass  by  Ge  atoms. 

The  -430cm1  band  attributed  to  the  Si-O-Si  symmetric 
stretching  mode  shifts  toward  lower  frequency  with  the 
increase  of  Ge02  content  in  the  glass.  This  band  also  shows 
a decrease  in  its  halfwidth  with  increasing  Ge02  concentration 
in  the  glass,  indicating  that  the  structural  fluctuation  in 
the  T-O-T  bond  angles  in  the  network,  where  T=Si  or  Ge, 
decreases  in  Ge02-rich  glasses [Wal75 ] . 

As  it  has  been  suggested [ Sha84  ] , the  presence  of  a single 
*'.(T-°-T)  stretch  band  in  the  spectra  of  binary  Ge02-Si02  glass 
indicates  that  the  modes  corresponding  to  j>s (Si-O-Si)  , ys(Si-0- 
Ge)  and  ^(Ge-O-Ge)  are  strongly  coupled.  This  experimentally 
supports  the  theoretical  prediction  of  Bell  and  Dean[Bel72] 
that  the  vibrational  modes  in  the  low  frequency  region  of  the 
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vibrational  spectra  of  the  network  of  tetrahedral  glasses  are 
delocalized. 

The  formation  of  Si-O-Ge  linkages  in  Ge02-Si02  glass 
network  in  the  films  studied  was  supported  by  the  observation 
of  the  formation  of  -675  cm'1  band. 

In  the  higher  frequency  region,  i.e.  850cm1-1200cm'1, 
because  of  the  small  amount  of  material  present  in  the  thin 
films  and  the  superposition  of  many  vibration  modes  in  this 
frequency  region,  a conclusive  analysis  is  difficult. 


CHAPTER  V 


Ge02-Si02  FILMS  BY  LPCVD  METHOD 
5.1  Introduction 

The  rapid  acceptance  of  LPCVD  in  the  microelectronic 
industry  was  due  to  its  extremely  favorable  economics  compared 
to  the  conventional  CVD  method.  The  major  reason  for  this 
cost  breakthrough  is  the  very  high  packing  density  available 
with  stand-up  closely-packed  substrate  loading. 

The  use  of  a low  pressure (0 . 1-1  torr)  process  results  in 
a reactor  condition  such  that  mass  transfer  limitations  at 
very  close  substrate  spacing (0. I»h0. 3")  are  unimportant 
compared  to  the  rate  of  the  chemical  reaction  at  the  substrate 
surface.  With  the  surface  reaction  rate  controlling  the 
deposition,  the  uniform  temperature  inherent  in  a resistance- 
heated  furnace  leads  to  excellent  uniformity  of  film  thickness 
and  composition. 

The  current  study  is  the  first  attempt  at  the  preparation 
of  planar  waveguide  of  Ge02-doped  silicate  glass  film,  using 
the  LPCVD  process.  In  this  chapter,  the  processing  conditions 
as  well  as  the  results  of  various  characterization  of  the 
prepared  films  will  be  presented. 
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5.2  Results  and  Discussions 
5.2.1  Film  Formation 

The  LPCVD  system  and  the  deposition  conditions  used  in 
the  present  study  have  been  described  in  details  in  section 
3.1.2.  Conducting  the  CVD  process  at  a low  pressure  has 
several  advantages  when  compared  to  atmospheric  CVD 
process [Tani71,  Ros77]: 

(1)  The  mean  free  path  of  a reactant  molecule  is  longer  and 
the  diffusivity  of  gaseous  species  is  greater. 

(2)  The  homogeneous  reaction ( i . e.  soot  formation)  may  be 
depressed  and  heterogeneous  reaction (or  surface  reaction) 
becomes  dominant. 

These  features  lead  to  the  deposition  of  uniform  thin  films  on 
the  substrates. 

The  details  of  the  surface  reaction  mechanisms  and 
kinetics  in  LPCVD  process  are  not  fully  understood.  in 

general,  the  gas  flow  within  the  reactor  was  adjusted  so  that 
the  chemical  reaction  is  depressed  in  the  vapor  phase (i.e. 
depressed  homogeneous  reaction)  and  a laminar  flow  region (or 
boundary  layer)  with  a thickness  of  l is  achieved  next  to  the 
hot  substrate,  as  shown  in  Figure  5.1.  Within  the  laminar 
flow  region,  the  molecules  transport  by  diffusion  due  to  the 
concentration  gradient  across  the  boundary  layer  and  the  flux 
of  certain  gas  species  toward  the  substrate 


can  be  expressed 
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bulk  gas  phase 
(transport  by  convection) 


i 

^ Laminar  flow  region 

Q. 

(or  boundary  layer) 

(transport  by  diffusion) 

hot  substrate  < 

Figure  5.1  Different  gas 
substrate  in  CVD  process. 


flow  regions  adjacent  to  a hot 
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as  follows: 


dc. 

J‘-~D 


(5.2.1) 


where  c is  the  concentration  of  molecular  species  i,  and  D;  is 
the  diffusion  coefficient  of  i.  Assuming  a linear 
concentration  gradient  across  the  boundary  layer  (5.1)  yields 

J;=  -Dj  (ACj/£ ) (5.2.2) 

where  £=thickness  of  the  boundary  layer.  For  a steady  state 
AC;  is  constant.  J;  is  dependent  upon  i and  D;.  Due  to  the 
increase  in  the  diffusivity  of  the  gaseous  species  at  low 
pressure  condition,  the  diffusion  flux  is  much  higher  compared 
to  that  of  atmospheric  CVD. 

During  the  deposition,  many  reactions  occur  which  can 
generally  be  illustrated  as  shown  in  Figure  5.2.  For  a 
temperature  range  between  350°C  to  500°C  the  overall  equations 
of  the  oxidation  reactions  are[Str68]: 


SiH4(g)  + 02W  ===>  Si02(s)  + 2H2(g) 

GeH4(g)  + 02(g)  ===>  Ge02(s)  + 2H2(g) 

According  to  the  works  of  various  researchers [Gra85,  Hen72, 
Yu72]  on  the  deposition  of  Si02  by  CVD  method. 


the 
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Figure  5.2  Various  reactions  occur  during  LPCVD 

deposition. 
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hypothetical  mechanisms  for  the  surface  reaction  of  silica 
film  on  pure  silica  substrate  could  be  illustrated  as  shown  in 
Fi.cfur'e  5.3.  Since  the  difference  in  silicon  and  oxygen 
electronegativities  is  large (1.8  for  Si  against  3.5  for 
0[Pau60] ) , the  Si- O bonds  have  a partially  ionic  character. 
At  room  temperature  the  Si  bonds  at  the  surface (dangling 
bonds)  are  saturated  with  hydroxyl  groups.  However,  at 
temperatures  above  300°C  the  surface  undergoes  a dehydration 
process [ Sny66 ] . Therefore,  the  silicon  dioxide  surface 
consists  of  positively  charged  Si  atoms  and  negatively  charged 
0 atoms.  Grabiec  and  PrzyLuski  hypothesized  that  silane 
reactions  (i.e.  adsorption (a)  followed  by  dissociation (b) , 
Oxygen  reactions  (i.e.  adsorption (e)  followed  by 
dissociation^ ) ) , and  dehydration (c  and  g)  occur 
independently.  For  the  hydride  LPCVD  process  to  prepare  Ge02- 
doped  Si02  films,  the  GeH4  molecule  behaves  similarly  to  the 
silane  molecule  and  forms  a binary  Ge02-Si02  glass  network. 

The  various  characteristics  of  the  LPCVD  films  are  summarized 
in  Table  5.1. 

5.2.2  Microstructure  of  the  Film 

Surface  morphology  of  the  LPCVD  films  was  examined  using 
SEM.  Figure  5.4  shows  the  SEM  micrographs  of  the  films.  All 
the  films  appear  to  consist  of  very  fine  particulates.  The 
LP1  and  LP3  films  have  similar  morphology  containing 
particulates  of  around  30nm  size 


range,  while  the  LP2  film 
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Figure  5.3  The  hypothetical  scheme  of  surface  reactions  in  SiO, 

LPCVD  process [Gra85 ] . 
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Table 

5.1  Characteristics 
Process. 

of 

Thin  Filins 

Prepared 

by  LPCVD 

Sample 

Appearance 

[Ge]/[Ge+Si] 

Index 

Thickness(um)  Abs.  band  Defect  structure 

Ipl 

light  brown 

4.24 

1.465 

0.660 

232nm 

185nm 

=Ge-Ge= 

Ip2 

slightly 

cloudy 

0.32 

1.451 

2.855 

238nm 

= Si-Si  E 

Ip3 

clear 

0.32 

1.468 

5.740 

238nm 

E Si-SiE 
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Figure  5. 
different 


SEM  micrographs  of  LPCVD  films  deposited  usinq 
processing  conditions. 
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contains  larger  features  around  lOOnxn  size  range.  The  fact 
that  the  LP2  film  contains  larger  particulates  is  believed  to 
cause  its  slightly  hazy  appearance  due  to  the  increased 
scattering  of  visible  light. 

The  microstructure  of  the  LP3  film  was  investigated  using 
TEM,  and  the  results  are  shown  in  Figure  5.5.  The  film 
consists  of  much  coarser  particulates  when  compared  to  the 
film  prepared  by  sol-gel  method,  see  Figure  4.15.  Some  large 
inclusions  in  the  film  have  also  been  observed,  as  shown  in 
Figure  5.5(b).  This  particle  size  variation  within  the  film 
could  enhance  the  scattering  loss  for  the  light  propagating  in 
the  film  when  it  is  used  for  waveguiding. 

5.2.3  XPS  Chemical  Analysis 

The  composition  of  LPCVD  films  deposited  using  various 
conditions  was  analyzed  by  using  XPS.  Figure  5.6  shows  the 
XPS  spectrum  of  LP1  film.  Only  Si,  o,  small  amount  of  Ge  and 
surface  carbon  were  detected  in  the  sample.  The  atomic 
concentration  of  Ge  and  Si  was  determined  by  calculating  the 
peak  area  of  the  Si2p  and  Ge3d  peak  respectively,  as  shown  in 
Figure  5.7(a)  and  (b)  . The  results  of  composition  analysis  of 
the  various  films  are  shown  in  Table  5.1.  The  LP1  film 
contains  4.24  at%  Ge  while  LP2  and  LP3  contain  only  0.32  at% 
of  Ge  in  the  film.  The  reason  for  the  low  Ge02  content  in 
these  films  is  believed  to  be  due  mainly  to  the  use  of  very 
dilute  germane  gas  (i.e.  0.2  wt%  GeH4/Ar)  and  pure  silane  gas 
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Figure  5.5  TEM  micrographs  of  LP3  film.  (a)  The 
contains  fine  particles  typically  around  lOOnm  scale. 
Large  particle  may  be  observed  in  some  areas. 


film 

(b) 


N(E>, 
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Figure  5.6 


XPS  composition  analysis  of  LP1  film 
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Figure  5.7  (a)  Si2p  and  (b)  Ge3d  band  in  XPS  spectrum  are  used 

to  calculate  the  atomic  concentration  of  the  LP3  film. 
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as  reactant  gases.  This  result  agrees  with  the  observation 
reported  by  Huffman  et  al.  [Huf86]  that  due  to  the  differences 
in  deposition  rates  in  LPCVD  process,  higher  Ge02  content 
silicate  films  are  difficult  to  obtain. 

5.2.4  UV/VIS  Spectroscopy 

Optical  losses  in  glass  waveguides  arise  from  absorption, 
Rayleigh  scattering,  and  waveguide  imperfections.  Absorption 
may  be  caused  by  molecular  transitions  involving  vibrational 
and  electronic  excitation  of  the  intrinsic  glass  material  and 
impurities  or  color  centers  that  are  incorporated  in  the 
processing.  In  general,  electronic  transitions  are  observed  in 
the  UV  region,  while  the  molecular  vibrations  dominate  in  the 
IR  region.  Absorption  in  the  UV-Vis  region (200~700nm)  by 
Ge02-doped  silicate  films,  prepared  by  LPCVD  method,  was 
investigated  in  the  current  study.  The  spectra  of  the  films 
are  compared  to  those  of  Suprasil  (Heraeous) , after  being 
normalized  to  the  same  thickness.  The  results  are  given  in 
Figure  5.8.  All  the  samples  have  an  absorption  band  around 
240nm,  namely  238nm  for  LP2  and  LP3  films,  and  232nm  for  LPl 
film.  The  origin  of  the  absorption  band  at  ~240nm  is  very 
complicated  despite  the  fact  that  it  has  been  observed  for  a 
long  time[Mit56,  Arn73,  Gar58a].  Tohmon  et  al.[Toh89]  studied 
pure  Si02  samples  using  UV  and  vacuum  UV  optical  absorption, 
photoluminescence (PL) , ESR,  and  MO  calculation.  These  authors 
found  that  there  are  two  types  of  absorption  bands 


near  245nm, 
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Figure  5.8  UV-Vis  spectra  of  various  LPCVD  films  and  fused 

silica  substrate (Suprasil) . The  spectra  are  normalized  to  the  same 
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each  with  unique  peak  wavelength,  half width,  PL  spectra,  and 
reaction  to  heat  treatment.  One  of  these  is  attributed  to  the 
oxygen  vacancy  defect (=si-si=) . The  existence  of  =Si-Sis  can 
cause  absorption  at  245nm  and  163nm  due  to  different 
electronic  transitions.  Furthermore,  Garino-Carina  first 
studied  extensively  the  band  near  240nm  observed  in  Ge02  and 
Ge02-Si02  glasses  and  attributed  this  band  to  the  presence  of 
GeO [Ga258 ] . This  conclusion  was  supported  by  Cohen  and 
Smith [ Coh58 ] . Although  it  is  agreed  that  the  absorption  at 
~240nm  in  Ge02  is  related  to  an  oxygen-deficient  defect,  the 
exact  defect  structure  is  still  an  issue  under  debate [Yue82 , 
Rad86 , Rau87 , Gri90] . However,  this  band  is  more  commonly 
refered  to  as  an  oxygen  vacancy  defect (sGe-Ges) , and  this 
oxygen  deficient  center  in  Ge02  or  Ge02-Si02  glasses  will  cause 
the  absorption  at  242nm  and  325nm[Row74,  Lev81],  The 

intensity  of  the  242nm  band  has  been  reported  to  be  nearly 
1000  times  more  intense  than  the  325nm  band[Yue82],  The 
absorption  at  ~l85nm  is  attributed  to  Ge02[ Yue82 ] . Therefore, 
the  absorption  band  for  LP1  is  assigned  to  oxygen  deficient 
defects  in  Ge02,  because  the  tail  of  the  absorption  at  ~185nm 
appears  with  the  232nm  band.  The  absorption  band(238nm)  for 
both  LP2  and  LP3  film  is  attributed  to  oxygen  deficient  defect 
in  Si02  because  the  tail  of  the  absorption  band  appears  to 

extend  further  down  into  vacuum  UV  range (most  likely  related 
to  163nm  band) . 
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5.2.5  Vibrational  Spectroscopy 

The  refractive  index  and  the  film  thickness  of  LPCVD 
samples  were  measured  using  ellipsometry  and  profilometry  and 
the  results  are  given  in  Table  5.1. 

Infrared  and  Raman  Spectroscopy  were  used  to  investigate 
the  structure  of  Ge02-doped  silicate  thin  films.  FTIR 
transmittance  spectra  of  various  films  were  obtained  using  an 
ATR  stage,  as  shown  in  Figure  5.9.  Because  the  Lpl  film  is  so 
thin  ( 0 . 66|xm)  its  spectrum  doesn't  show  features  corresponding 
to  Ge02.  The  spectra  of  both  LP2  and  LP3  have  a broad  and 
small  absorption  band  between  515cm1  to  585cm1,  corresponding 
to  Ge-O-Ge  symmetric  stretching  mode[Bor69].  The  formation  of 
Si-O-Ge  linkage  in  these  Ge02— Si02  films  can  be  supported  by 
the  observation  of  the  peak  around  1000cm1 (i. e.  960-1020cm1) 
which  has  been  assigned  to  Si-O-Ge  antisymmetric 
stretching [Bor69] . 

Raman  spectroscopy  was  also  used  to  study  the 
distribution  of  tetrahedra  in  Lp3  waveguiding  film.  The 
result  was  shown  in  Figure  5.10.  Compared  to  Suprasil  pure 
silica  and  Si02  film(made  by  sol-gel  method) , LP3  film 
consists  of  structures  that  are  very  similar  to  those  of  Si02 
after  neutron  irradiation,  as  shown  in  Figure  5.11[Sto70], 
Basically  there  are  some  important  features  in  LPCVD  film: 

(1) . The  major  low-frequency  band  lies  near  470cm1  and  is 
considerably  narrower  than  in  bulk  Si02. 
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Figure  5.9 


FTIR  transmittance  spectra  of  LPCVD  thin  films 


Normalized  Intensity 
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Figure  5.10  Comparison  of  Raman  spectra  of  the  LP3  film 
sol-gel  Si02  film  and  the  fused  silica  (Suprasil) . 


the 
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Figure  5.11  Room  temperature  Stokes  Raman  spectra  for  (A) 
unirradiated  and  (B)  neutron- irradiated  Optosil  fused  silica 
[ Sto7  0 ] . 
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(2) .  The  defect  bands  at  490cm'1(Dl  band)  and  606cm*1  (D2  band) 

appear  to  be  the  dominant  structure. 

(3) .  A peak  appears  at  970cm1  corresponding  to  the  Si-OH 

stretching  mode. 

The  decrease  of  the  halfwidth  of  the  main  peak  indicates 
a narrower  distribution  of  the  T-O-T  bond  angles  in  the 
network,  where  T=Si  or  Ge[Wal75,  Huf85]. 

Although  the  structural  model  for  D1  and  D2  defect  bands 
are  still  controversial  in  the  literature [Rev83 , Sha81,  Phi82, 
Sil84 , Gal85 ] , it  is  more  commonly  believed  that  D1  is  related 
to  the  symmetric  stretching  mode  of  4-membered  ring  of 
tetrahedra [ Gal8 5 ] and  D2  is  related  to  the  breathing  mode  of 
a planar  3-fold  of  bonds [Gal85 ] . The  possible  reason  for  the 
dominance  of  these  defect  structures  in  the  LPCVD  film  is  the 
low  deposition  temperature  used  in  the  process (e.g.  450°C  far 
lower  than  800-1150  °C  used  than  in  the  sol-gel  films)  . 
Therefore  the  material  doesn't  have  enough  energy  to  relax  its 
defects  to  a structure  that  has  a lower  energy.  According  to 
Huffman  and  McMillan's  work  on  deposition  of  Si02  by  LPCVD 
method [Huff 85 ] , the  structure  of  as-deposited  Si02  may  relax 
to  near  the  bulk  glass  value  and  distribution  on  annealing  at 
900°C  for  2 to  3 days. 

The  existence  of  SiOH  in  the  LP3  film  is  believed  to  be 
caused  by  the  reaction  of  hydrogen (H2)  , which  could  be  formed 
after  the  oxidation  of  hydride  reactants,  with  Si02-Ge02  glass 
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film.  The  formation  of  hydroxyl  groups (SiOH  and/or  GeOH)  in 
optical  fibers  due  to  the  reaction  with  H2  have  been  reported 
extensively  in  the  literature [Ues83 , Tan84,  Ple84,  Tani84, 
Pit84  ] . 


CHAPTER  VI 


PLANAR  WAVEGUIDE  CHARACTERIZATION 
6 . 1 Introduction 

Important  waveguide  characteristics  such  as  mode  index 
measurement,  propagation  loss  as  well  as  mode  profile 
measurement  are  investigated  by  using  a prism  coupler  as  a 
light  coupling  technigue  as  was  described  in  section  3.2.10. 

The  planar  waveguides  used  in  the  present  study  are 
prepared  by  using  both  sol-gel  and  LPCVD  processes.  The  film 
properties,  such  as  thickness,  refractive  index,  refractive 
index  difference (An)  between  film  and  substrate,  and 
processing  condition  are  reported  in  Table  6.1.  The 
refractive  index  difference (An)  is  calculated  by  using  the 
following  equation [Nat82 ] 


An  = 


ni  n2 


2 nl 


(6.1.1) 


where  n,=  index  of  film(i.e.  core  layer),  n2=  index  of  v-Si02 
substrate ( i . e . 1.457  for  632. 8nm).  The  value  An  can  be  used 
to  calculate  an  effective  numerical  aperture (NA)  for  the 
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Table  6.1  Planar  Waveguides  Fabricated  Using  Both  Sol-Gel  and 
LPCVD  Method. 


Sample 

Process 

Temp. 

(°C) 

Number  of 
dippings 

Thickness 

(um) 

Index 

(n,) 

An 

(%) 

10GS 

sol-gel 

1110 

72 

2.615 

1.485 

1.87 

30GS 

sol-gel 

900 

56 

1.570 

1.502 

3.01 

50GS 

sol-gel 

700 

22 

0.667 

1.550 

5.82 

LP3 

LPCVD 

450 

— 

5.74 

1.468 

0.75 

waveguide  with  the  following  equation 
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NA  = 


(6.1.2) 


6.2  Mode  Index  Measurement 

If  we  represent  the  incident  wave  on  the  prism  by  a wave 
vector  K0/  its  propagation  in  the  prism  can  be  characterized 
by  a magnitude  k0rip  which  can  be  decomposed  into  a horizontal 
component  k0npCosc  and  a vertical  component  kQnpSinc,  as  shown 
in  Figure  6.1.  The  boundary  conditions  of  the  electromagnetic 
fields  at  the  prism  base  require  that  the  field  below  and 
above  the  prism  base  have  the  same  horizontal  wave  component. 
Therefore,  maximum  tunneling  interaction  takes  place  when  the 
horizontal  component  of  the  wave  vector (along  the  z direction) 
is  equal  to  the  propagation  constant  of  a given  mode  among  the 
many  possible  modes  of  the  waveguide ( i. e.  phase  matching 
condition  is  satisfied) . Under  this  condition  that  particular 
mode  gets  excited [Ulr73 ] . 

Thus, 


0 =kanpcosc 


(6.2.1) 


However , 
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np  > nf  > ns 


Figure  6.1 


Ray  diagram  for  the  prism  coupler. 
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J 3=kaN 


(6.2.2) 


so  that 


N =npcosc 


(6.2.3) 


where  N is  the  mode  index  (or  effective  index)  of  a given  mode, 
rip-  refractive  index  of  the  coupling  prism.  Now  considering 
A ABC  within  the  coupling  prism,  as  shown  in  Figure  6.1. 


so 


and 


/a+/c+/r  + 90°  = 180° 


/c  = 90°-/a-/ 


ripCos  c = rip  sin(/  a+/_  r) 


(6.2.4) 


(6.2.5) 


(6.2.6) 


From  Snells'  law  at  point  B 


and 


sin  Im  = np  sin  r 


Z r = sin1  (sinlm/np) 


(6.2.7) 


(6.2.8) 


Thus, 
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N = rip  sin  [a  + sin1  (sinlm/np)  ] (6.2.9) 

The  mode  indices (NJ  of  both  30GS  and  LP3  waveguides  are 
calculated  from  the  experimentally  measured  angle  of 
incidence (IJ  for  various  modes  and  the  results  are  given  in 
Table  6.2  under  the  column  Im  and  Nm/  respectively.  By 
determining  these  strong  coupling  angles,  we  can  find 
experimentally  the  characteristic  propagation  constants  of 
various  propagation  modes  in  a given  film,  relative  to  the 
propagation  constant  of  free  space,  k0=  oj/c,  according  to  the 
following  relation[Tie71] 

6=  Mm  = K nf  sin#  (6.2.10) 

where  the  0 is  the  ray  propagating  angle  within  the  film  as 
shown  in  Figure  6.6.  This  measurement  will  also  lead  to  the 
determination  of  the  phase  velocity,  vm,  of  the  various 
characteristic  modes  of  propagation  in  the  guide [Ulr73]  by 
using  the  following  equation 

vm  = C/Nm  = u/6  (6.2.11) 

For  a planar  waveguide,  in  order  to  guide  the  light  in 
the  structure  the  transverse  resonance  condition  has  to  be 
met[Tie71,  Kog79]  as  has  been  described  in  section  1.2.1. 
That  is  the  total  phase  shifts  introduced 


as  the  ray  travels 
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Table  6.2  Characteristics  of  the  Sio2-Geo2  Waveguides. 


SAMPLE 

MODE 

1 

m 

N _ 
m 

N calc. 

A7E 

EF 

30GS 

TEO 

19.833° 

1.494 

1.493 

7.947 

0.970 

TE1 

17.300° 

1.469 

1.467 

7.947 

0.824 

TEO 

16.833° 

1.465 

1.467 

34.898 

0.993 

LP3 

TE1 

16.417° 

1.461 

1.465 

34.898 

0.969 

TE2 

16.000° 

1.456 

1.461 

34.898 

0.905 
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a full  zig-zag  cycle  need  to  be  an  integral  multiple  of  2n , as 
shown  below 


where  kQ=  2n/\0  = propagation  constant  in  free  space, 

X0=  free  space  wavelength, 

k0nf  cos 6 = propagation  constant  in  waveguide, 
m=  integer (or  mode  number) 

<p=  phase  shift  occurs  at  film/cladding  interface, 

<p=  phase  shift  occurs  at  film/ substrate  interface. 

For  TE  modes,  the  phase  shift  0C  and  <pt  can  be 
expressed [Kog7 9]  in  terms  of  film  parameters  as 


For  a set  of  given  film  parameters  of  nf,  ns,  nc,  and  d, 
and  a given  wavelength (Xj  there  may  have  different  values  of 
9,  corresponding  to  different  mode  numbers (m)  for  which  the 
eguation  6-14  is  satisfied.  For  each  value  of  9 there  is  a 


2dk0nf  cos 6 = 2mn  + 2 0,  + 20, 


(6.2.12) 


(6.2.13) 


(6.2.14) 
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corresponding  value  of  incident  angle (Im)  at  which  a given 
mode  is  excited  in  the  waveguide [ Kog74 ] . 

As  mentioned  in  Equation  (6.2.10),  the  propagation 
constant  for  a guided  mode  (3  is  given  as 


and  its  value  lies  between  the  plane  wave  propagation 
constants  in  the  substrate  and  the  film[Kog74] 


When  /3= k0ns/  the  mode  is  said  to  be  cutoff.  Mode  results  for 
slab  waveguides  are  generally  expressed  in  terms  of  normalized 
dimensionless  parameters [Kog74 ] V and  b.  The  V parameter  is 
called  the  normalized  frequency ( or  normalized  film  thickness) 
and  is  defined  as 


k0n,  < jS  < k0nf 


(6.2.15) 


X 


and  b is  called  the  normalized  mode  index (or  normalized 
propagation  constant)  and  is  defined  as 
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b = 


Wm  ~ns 

n}-nt ? 


(6.2.16) 


The  index  b is  0 at  cutoff(i.e.  Nm  = n.)  and  approaches  1 when 
far  away  from  it(i.e.  Nm=  nf)  . 

A measure  of  the  asymmetry  of  the  waveguide  structure  is 
also  def ined[Kog74]  for  both  TE  and  TM  modes  as  follows 


aTE  ~ 


2 2 
^ s nc 

~2  ~2 

nf  ~ns 


(6.2.17) 


aTM 


( nf}4  ( ns~nc ) 

n<  (nj-n]) 


(6.2.18) 


The  value  of  a ranges  from  0(i.e.  perfect  symmetry,  n=nc)  to 
infinity(i.e.  strong  asymmetry  nanf)  . Combining  Equation 
(6. 2. 13) -(6.2.18)  and  substituting  into  Equation  (6.2.12),  the 
normalized  dispersion  relation  may  be  expressed  as[Kog74] 


\Jv  (1-b)  =m7r  + tan'1 


1-b 


+ tan 


-l 


a +jb 

X=jb 


(6.2.19) 


For  a given  value  of  m the  normalized  mode  index  b depends 
only  on  V and  a.  The  calculated  normalized  film  thickness  as 
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a function  of  normalized  mode  index ( or  b-V  characteristics) 
for  the  first  four  TE  modes  of  3 0GS  sol  gel  waveguide  (aTC= 
7.947)  is  given  in  Figure  6.2.  According  to  the  film 
parameters  of  30GS  film  (n^l.502,  n =1.454,  n=l,  d=1.566/un, 
and  Xo=0.6328Mm) , V parameter  is  calculated  to  be  5.845.  A 
dotted  line  is  drawn  across  the  plot  at  constant  value  of 
5.845,  the  two  intercepts  with  m=0  and  m=l  curves  give  the 
expected  b values  for  these  two  modes  for  30GS  waveguide  which 
in  turn  may  be  used  to  obtain  Ncal  utilizing  Equation  (6.2.16). 
The  calculated  mode  indices  (Ncal  ) according  to  the  step-index 
waveguide  model,  i.e.  Equation  (6.2.18),  are  listed  in  Table 
6.2.  The  calculated  mode  indices (N^)  are  very  close  to  the 
experimentally  measured  values (Nm). 

Similarly,  the  normalized  film  thickness (V)  versus 
normalized  mode  index (b)  curves  for  the  first  four  possible  TE 
modes  in  LP3  waveguide (aTE=34 . 898)  are  calculated  and  plotted 
in  Figure  6.3.  According  to  the  LP3  film  parameters  (n^l.468, 
n =1.457,  n=l,  d=5 . 74  nm,  X =0 . 6328/ixn)  , the  V parameter  is 
calculated  as  10.233.  By  drawing  a line  across  the  plot  at 
constant  value  of  V=10.233,  the  line  intercepts  first  three 
curves (i.e.  m=0,l,and  2),  indicating  the  three  possible  TE 
modes  that  can  propagate  in  LP3  waveguide.  The  calculated 
mode  indices  (N^  ) for  LP3  sample  can  be  obtained  from  the 
corresponding  b values  on  the  plot  and  using  Equation 
(6.2.16).  The  results  of  NCJl  are  given  in  Table  6.2  and  are 
very  comparable  with  those  obtained  experimentally.  This 
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6‘ /ttx  Normalized  mode  index  (b)  versus  normalized  film 
thickness (V)  characteristics  for  the  30GS  film,  modeled  as  a 
step-index  waveguide. 


Normalized  Film  Thickness(V) 
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Figure  6.3  b-V  characteristics  of  the 
a step  index  waveguide. 


LP3  film  modeled 


as 
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indicates  that  the  waveguides (both  sol-gel  and  LPCVD  samples) 
used  in  present  study  may  possess  approximately  a step- 
refractive  index  profile. 

6.3  Propagation  Loss  Measurement 

In  this  study,  a three-prism  loss  measurement 
method [Won80 ] , in  which  the  results  are  independent  of  the 
prism  coupling  coefficients,  is  used  to  measure  the 
propagation  loss  in  the  waveguide.  During  the  measurement, 
prisms  l and  3 are  fixed  as  shown  in  Figure  6.4,  while  prism 
2 is  moved  along  the  guide  in  the  z direction.  The  output 
power  from  prism  2 and  3 are  detected  by  photodetectors. 

Considering  P2  and  P3  as  the  output  powers  of  prism  2 and 

3,  respectively,  r2  and  r3  are  the  coupling  coefficients  at  the 
output  prisms,  then 


P2=r2P(z2) 


(6.3.1) 


P3=r3P(z3) 


(6.3.2) 


where  P(z2)=  power  propagating  at  point  z2. 

Let  a be  the  attenuation  coefficient  in  the  waveguide,  then 
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/N  X 


input  prism 


O — > z 


Figure  6.4 
measurement 


Three-prism  method 


for  propagation  loss 
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(6.3.3) 

P(z3)  =[P(z2)  - P2  ] e 

(6.3.4) 

P3-r3[P(z2)  -r2P(z2)  ]e'ad 


For  the  case  that  r2— 0,  i.e.  prism  2 is  removed  and  the  power 
output  at  prism  3 is  P30: 


P30~r3[P(z2)  ]e~ad 


(6.3.5) 


M>3  = P3-P3o=r3r2P(z2)e-“i 


(6.3.6) 


P3P3o 


A P 


_r2P(z2)r3P(z2)e~arf 

r^PCZzJe-^ 

= P(z2) 


(6.3.7) 


From  the  last  equation  , it  is  seen  that  the  expression  for 
F5 ( z2 ) / the  power  propagating  at  z2,  has  no  r dependence  in  it. 
Therefore,  P(z2)  can  be  estimated  without  the  knowledge  of 
coupling  coefficients,  r2  and  r3.  Thus  by  moving  the  prism  2, 
i.e.  by  changing  z2,  one  can  find  the  power  along  the 
different  points  in  the  waveguide.  The  loss  rate  was  obtained 
from  the  plot  of  relative  power (dB)  vs.  distance,  and  the 
relative  power  was  calculated  by  using  the  following  equation: 
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Prel  (dB)  = 10  log (Z)  (6.3.8) 

Po 

where  P0  is  the  propagating  power  level  at  the  reference 
point,  and  P is  the  various  propagating  powers  at  different 
distances  along  the  waveguide. 

Figu^-"®  6.5  shows  the  plot  of  relative  power  level  versus 
propagation  distance  for  both  50GS  and  LP3  waveguide.  The 
slope  of  each  curve  gives  the  propagation  loss  rate.  The  50GS 
sol-gel  glass  waveguide  has  a loss  rate  of  3.31  dB/cm,  and  the 
LP3  LPCVD  glass  waveguide  has  a loss  rate  of  2.59  dB/cm. 
These  loss  rates  are  comparable  to  those  have  been  reported  in 
the  literature  for  various  thin  films  used  as  waveguides  in 
integrated  optics [Ostr78 ] , as  shown  in  Table  6.3. 

The  result  from  the  sol-gel  film  is  especially 
encouraging  and  open  a new  application  for  this  multiple  dip- 
coating approach.  The  film  exhibits  good  optical  quality  and 
low  loss,  although  each  layer  was  prepared  in  an  open 
environment  without  special  substrate  preparation,  solution 
filtration,  and  clean  room  facilities. 

There  are  two  major  types  of  losses  which  may  contribute 
to  the  propagation  loss  in  optical  waveguides:  absorption 
losses  and  scattering  losses.  The  materials  prepared  in  the 
present  study  are  of  very  high  purity  and  the  major  impurity 
ions  are  the  hydroxyl  ions,  which  don't  absorb  at  632. 8nm. 
Moreover,  the  waveguide  length  used  in  the  present  study  is 
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Figure  6 . 5 
and  the  LP3 


Propagation  loss  measurement  for  both  the  50GS 
films  at  632. 8nm. 


Table  6.3  Thin  Filins  for  Integrated  Optics (Ost78)  . 
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Material 

Form 

Fabrication 

technique 

Index 
(6328  A) 

Loss  at 
6328  A 
(dB/cm) 

7059  glass 

Amorphous 

Sputtering 

1.53-1.59 

1-2 

1-5 

1 

1 

20 

Barium  silicate 

Amorphous 

Sputtering 

1.57-1.67 

Ta,Os 

Amorphous 

Sputtering 

2.214 

Nb.Oj 

Amorphous 

Sputtering 

2.276 

ZnO 

Polycrystalline 

Sputtering 

1.97 

ZnS 

Polycrystalline 

evaporation 

2.342 

5 

LiTaO, 

Polycrystalline 

Sputtering 

1.95-2.08 

2 

Glass 

Aluminosilicate 

Amorphous 

AgNOj  exchange 

1.52-1.65 

0.1 

glass 

Amorphous 

KNOj  exchange 

1.52-1.56 

0.1 

Silica 

Amorphous 

Li  implantation 

1.47 

0.2 

Characteristics  of  Organic  Films 


Material 

Fabrication 

technique 

Index 
6328  A 

Loss  at 
6328  A 
(dB/cm) 

Polyurethane 

Spin  or  dip  coating 

1 .54- 1 .58 

1 

0.3 

Epoxy  (Araldite) 

Spin  or  dip  coating 

1 .581 

KPR 

Spin  coating 

1 615 

7 

Vinyltrimethylsilane 

Polvphenylsiloxane 

Plasma  polymerization 
S pi  neon  ting 

1.537 
1 ‘'b 

0.0! 
i i 

Methyl/glycidyl 

methacrylate 

Dip  coating 

1.515 

1 -4 
0.2 

183 


much  shorter  compared  to  that  usually  used  in  optical  fiber 
studies.  Therefore,  the  major  cause  for  the  loss  in  these 
samples  is  attributed  to  scattering  processes.  Several  types 

of  scattering  processes  may  contribute  to  the  propagation  loss 
in  these  cases: 

(1).  Intrinsic  Rayleigh  scattering (y^)  which  is  caused 
by  the  refractive  index  fluctuations (e.g.  fluctuations  in 
density,  composition,  or  structure)  that  are  small  with 
respect  to  wavelength  of  light,  and  can  be  expressed  in  terms 
of  fundamental  material  parameters [Tic90] 


( B2Z2jTg ) 

XX4 


(6.3.9) 


where  Z,  is  the  anion  valence,  0 is  the  isothermal 
compressibility,  Tg  the  glass  transition  temperature,  E0  the 

average  energy  gap,  and  B is  a dimensionless  structure  factor 
given  by 


B=Nad 3 


(6.3.10) 


where  N.  is  the  volume  density,  and  d is  the  first-neighbor 
bond  length. 

(2).  Volume  scattering  which  is  caused  by 
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inhomogeneities  such  as  pinholes,  micropores,  and  interfaces, 
within  the  volume  of  the  waveguide,  that  introduce  variations 
refractive  index  along  the  path  of  the  propagating  light. 
The  variation  of  refractive  index  causes  the  deflection  of 
light,  and  the  stray  light  may  leak  out  of  the  waveguide 
because  the  change  in  propagation  direction  is  no  longer 
confined  to  the  waveguide  structure [Kin76 ] . 

(3)  . Surface  scattering  that  is  caused  by  irregularities 
in  core/cladding  boundary.  When  the  boundary  is  not  perfectly 
flat,  the  light  traveling  in  the  core  is  reflected 
irregularly.  Some  of  the  rays  come  off  at  different  angles, 
giving  rise  to  mode  conversion,  while  others  radiate  away  into 
the  cladding  because  their  angles  of  incidence  exceed  the 
critical  angle  for  total  ref lection[Iga90] . 

According  to  the  electron  microscopy  study,  which  was 
presented  in  last  two  chapters,  the  films  prepared  by  both 
sol-gel  and  LPCVD  method  are  composed  of  small  particulates 
instead  of  containing  a homogeneous  and  continuous  phase. 
Therefore  interfaces  between  particles,  and  surface 
irregularities  have  been  observed.  Moreover,  pinholes  and 
some  large  inclusions  that  were  introduced  into  the  film 

during  film  deposition  may  act  as  scattering  centers  in  the 
waveguide. 

6.4  Mode  Profile  Measurement 

To  calculate  the  field  profiles.  Maxwell's  equations  have 
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to  be  used  under  the  appropriate  boundary  conditions. 
Considering  a three-layer  step-index  planar  waveguide  with 
light  propagating  in  yz  plane  along  z direction  with  a 
propagation  constant  j8  (or  with  a phase  velocity  vp=«/j8)  , as 
is  shown  in  Figure  6.6.  By  assuming  the  light  wave  in  the 
film  to  be  infinitely  wide  in  the  y direction  so  that  d / dy  =0, 
the  transverse  electric  field  component  for  the  TE  mode  Ey  can 
be  shown  to  satisfy  the  transverse  Helmholz  eguation[Tie71] 

-^3  + (k2on2-02)E  =0  (6.4.1) 

L y 


The  boundary  conditions  imposed  by  the  waveguide  structure 
demand  that  Ey(and  thereby  automatically  Hx)  and  dEy/dx(and 
thereby  Hz)  be  continuous  across  the  film  boundaries  at  x=0 
and  x=d [Kog79 ] . For  guided  TE  modes,  the  solution  of 
Equ. (6.4.1)  will  be  an  exponentially  decaying  function  for 
both  x<0  and  x>d,  and  an  oscillatory  function  within  d>x>0, 
which  can  be  expressed  as  follows [Kog79 ] 

EyC(x)=  Ef  cos0c  e0"  for  x<0  (cladding)  (6.4.2) 

Eyf(x)=  Ef  cos  (kx-0c)  for  o<x<d  (film)  (6.4.3) 


Ey,(x)=  Ef  cos  (kd-0c)  e'oa(x-d) 


for  x>d  (substrate) (6.4.4) 
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y O — > z 


V 

X 


cladding  n 


x=0 


x=d 


Figure  6 . 6 


Ray  path  in  a three  layer  step  index  waveguide. 
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where  ac=  absorption  coefficient  of  the  cladding  layer,  and  a= 
absorption  coefficient  of  the  substrate.  The  various 
constants  may  be  expressed  in  terms  of  film  parameters  (nf,  ns, 
nc,  d)  [Ram89] 


=kJn]-N2 


ac=kjN2-n 


as  = k^N2-nl 


(6.4.5) 


(6.4.6) 


(6.4.7) 


Recalling  the  definition  of  various  normalized  parameters (i . e. 
b and  V)  and  asymmetric  parameter (a^) 


b = 


N2-nl 


0<b<l 


(6.4.8) 


V‘k „d  yjnj-rf 


(6.4.9) 


a 


TE  ~ 


n f -n 


2 

c 

1 

s 


0<a<oo 


(6.4.10) 
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and  by  defining  x'=x/d  and  assuming  Ef=l,  the  Equ. (6.4.2) 
-(6.4.4)  can  be  expressed  as  normalized  field  distribution 
equations  in  terms  of  these  parameters (i . e.  V,  a,  b,  and 
x' ) [ Ram8  9 ] 


Eyc(x)  =cos<pcevi^x'  for  x'<0  ( cladding ) 


(6.4.11) 


Eyf{x)  =cos(0mx/ -0c)  for  0<x'<l  {film) 


(6.4.12) 


Eys(x)  = (-l)Mcos<^evJB'(1"l,> 


(6.4.13) 

for  x'>\  ( substrate ) 


Various  phase  shift  terms  (i.e.  <pm,  <p%,  and  0C)  in  Equ. 

(6.4.11)-  (6.4.13)  may  also  be  expressed  in  terms  of  film 

parameters [ Ram8  9 ] 


<pm=mn  +<ps  + (pc 


(6.4.14) 


0,=tan->[  (1+g)  ( 

> 


)] 


1-jb 


(6.4.15) 
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0e  = tan-1[(_L^)  ( 

1 -aq 

where  q=0  for  TE  modes,  and  q=(nf2/n,2)  -1  for  TM  modes. 
Therefore  the  normalized  field  distribution  equations 
(6. 4. 11) -(6. 4. 13)  can  be  expressed  completely  in  terms  of  film 
parameters [Ram8 9 ] 


a+b  ^ 1 (6.4.16) 

T^b  1 1 


E^x) 


1-b 
1 +a 


0 V'J  15  +bj~  x ' 


for  X' <0  ( cladding ) 


(6.4.17) 


E^x)  = [ 


cos(Wd-ib) 


x') 


a+b 
1 +a 


for  (Kx'ci 


sin(vV(l-b)  x')  ] 

(film)  (6.4.18) 


Eys(x)  - (-l)m\/l  -b  (l'*)  for  x'>  1 ( substrate ) 

(6.4.19) 

By  inputing  various  film  parameters (nf,  n„,  nc,  d,  and  X0)  , 
which  may  be  obtained  experimentally,  the  field  distribution 
across  the  transverse  direction (i.e.  x direction) , both  within 
and  beyond  the  waveguide  structure  can  be  calculated  using 
Equ. (6.4.17)-  (6.4.19).  The  results  of  calculated  field 


190 


distribution  for  TEO  and  TE1  modes  of  3 0GS  waveguide  are  given 
in  Figure  6.7. 

The  instantaneous  power  flux,  both  magnitude  and 
direction,  of  electromagnetic  energy  is  usually  described  by 
the  Poynting  vector  (pj  [Lee86] . For  a step-index  slab 
waveguide,  the  Poynting  vector  for  TE  mode  can  be  expressed  as 

P=JL\e\2  (6.4.20) 

which  is  basically  taking  the  square  of  the  field  distribution 
functions  and  then  multiplying  with  a constant. 

The  total  power  flow  in  a dielectric  waveguide  can  be 
obtained  by  integrating  from  -«  to  +oo  in  x direction, and  it 
can  be  expressed  as  follows[Lee86] 


p< = 2^, 1 r, 1 v 1 2ax  * jo"  i i ■ e*i  * jr  i e,.\ i i < * • 4 • 2 1 > 

The  mode  profile  of  the  planar  waveguide  is  measured 
using  a near-field  method,  as  described  in  section  3.2.10. 
This  method  consists  of  measuring  the  intensity  distribution 
along  the  transverse  direction  (i.e.  x direction)  of  the  near- 
field image  of  a propagating  mode,  which  was  projected  to  the 
image  plane  from  the  end  face  of  the  waveguide  by  using  a 
microscope  objective.  The  results  for  the  experimentally 
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Figure  6.7  Calculated  field  profiles  for  m=0,  and  1 TE 
inodes  of  the  30GS  film  modeled  as  a step-index  waveguide. 
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measured  TEO  and  TE1  mode  are  shown  in  Figure  6.8  and  6.9, 
respectively.  The  fluctuation  of  the  intensity  in  Figure  6.8 
near  the  edge  of  the  waveguide  is  attributed  to  the  scattering 
of  light  by  the  defect  at  the  end  face  of  the  waveguide 
introduced  during  polishing.  The  intensity  profile  of  TE1 
mode,  see  Figure  6.9,  contains  two  peaks  as  expected  from  the 
results  of  the  calculated  mode  profile,  as  shown  in  Figure 
6.10.  However,  the  relative  intensities  of  these  two  peaks 
are  not  comparable  as  predicted  from  Figure  6.10,  which  is 
based  on  a perfect  step-index  profile  assumption.  This 
indicates  a slight  variation  of  the  index  profile  across  the 
actual  waveguide.  The  compositional  uniformity  across  the 
30GS  film  was  analyzed  by  a Scanning  Auger  Electron 
Spectroscopy  (AES)  , with  Ar  ion  gun  sputtering  for  depth 
profiling.  The  intensity  of  Si  , Ge,  and  O were  monitored 
with  sputtering  time.  All  the  elements  studied  appear  to  be 
guite  constant  across  the  film  even  after  five  hours  of 
sputtering , as  shown  in  Figure  6.11.  However,  this  analysis 
was  complicated  by  the  charging  problem  that  is  commonly 
encountered  for  dielectric  materials,  making  extensive 
sputtering  studies  a difficult  task.  A possible  reason  for 
the  deviation  from  a perfect  step-index  profile  may  be  due  to 
film  microstructure,  which  may  vary  between  the  layers  close 
to  and  far  from  the  substrate.  This  could  be  examined  by 
cross-sectional  TEM.  However,  the  problem  was  not  considered 
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Figure  6.8 


Experimentally  measured  TEO  mode  prof  ile  (350x)  . 
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Figure  6 . 9 


Experimentally  measured  TE1  mode  prof  ile(350x)  . 
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Figure 


Calculated  mode  profiles  for  the  30GS  film  as 
a perfect  step-index  slab  waveguide. 
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Figure  6.11  Scanning  Auger 
analysis  of  the  30GS  film. 
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significant. 

Similarly,  the  field  profile  for  TEO,  TE1 , and  TE2  modes 
in  LP3  film  have  been  calculated  and  are  shown  in  Figure  6.12. 
Each  mode  profile  is  experimentally  measured  and  the  results 
are  given  in  Figure  6.13,  6.14,  and  6.15  for  TEO,  TE1,  and  TE2 
mode  respectively.  These  mode  profiles  appear  as  predicted 
by  calculated  results  of  mode  profiles,  i.e.  one  peak  for 
fundamental  mode,  two  peaks  for  TE1  mode,  and  three  peaks  for 
TE2  mode  as  is  shown  in  Figure  6.16.  The  relative  intensity 
of  the  two  peaks  in  TE1  mode  profile  are  not  comparable  as 
predicted  by  the  step-index  model.  Similar  deviation  is  also 
observed  in  the  measurement  for  TE2  mode,  as  shown  in  Figure 
6.15. 

The  power  confinement  efficiency (EF)  of  each  mode  is 
obtained  by  calculating  the  ratio  between  the  area  below  the 
intensity  profile  within  the  waveguide  to  that  of  the  total 
intensity  profile  (i.e.  including  cladding,  waveguide,  and 
substrate)  which  may  be  expressed  as 


(1 

EF  = J 

0Vj2 

(6.4.22) 

<J 

The  results  are  listed  in  Table  6.2  under  the  column  EF. 

The  fact  that  EF  decreases  with  increasing  mode  number  is 
because  more  power  is  extending  out  beyond  the  boundaries  of 
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Figure  6.12  Calculated  field  profiles  of  TEO,  TE1  and  TE2 
inodes  for  the  LP3  film  as  a step  index  waveguide. 
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Figure  6 . 13 


Experimentally  measured  TEO  mode  prof  ile (350x) 
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Figure  6.14  Experimentally  measured  TE1  mode  prof ile (350x)  . 
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Figure  6.15  Experimentally  measured  TE2  mode  prof ile (350x)  . 
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Figure  6.16  Calculated  mode  profiles  for  TEO,  TE1,  and  TE2 
modes  of  the  LP3  film  modeled  as  a perfect  step— index  slab 
waveguide. 
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the  guiding  film,  which  is  clearly  shown  in  the  calculated 
mode  profile  plots  in  Figure  6.10  and  Figure  6.16  for  30GS  and 
LP3  waveguide,  respectively. 


CHAPTER  VII 


NONLINEAR  OPTICAL  MEASUREMENT 
7 . 1 Introduction 

The  nonlinear  properties  of  glasses  are  of  paramount 
significance  in  areas  where  light  intensities  are  in  the  order 
of  l GW/cm2.  Two  of  these  areas  are  high  power  laser- 
fusion  (Fel73 ) and  all-optical  switching  for  information 
processing ( Fri87 ) . In  the  former  the  nonlinear  properties 
have  to  be  reduced  (n2  must  be  small)  in  order  to  limit  the 
damage  to  optical  glass  components  by  self -induced  focusing. 
For  the  latter  case  there  is  an  increasing  need  for  trans- 
parent materials  in  the  1-2  jum  region  with  highly  intensity 
dependent  nonlinear  refractive  index  coefficients  which  can 
withstand  high  optical  peak  powers. 

The  nonlinear  optical  phenomena  in  glasses  may  be  divided 
into  two  categories.  The  first  category  includes  phenomena  in 
which  the  glass  itself  is  responsible  for  the  nonlinearity 
( intrinsic  effect) , such  as  nonlinear  refractive  index, 
multi-Photon  absorption,  and  stimulated  Raman  and  Brillouin 
scattering.  The  other  category  includes  glasses  containing 
rare  earth  ions,  semiconductor  microcrystals,  or  organic  dyes, 
where  the  nonlinearity  is  associated  with  the  dopant  and  the 
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glass  acts  as  a host  (extrinsic  effect)  and  its  contribution  is 
secondary.  For  intrinsic  effects  the  composition  and 
structure  of  the  glass  are  of  primary  importance  in 
controlling  the  nonlinearity. 

In  the  present  study,  we  investigate  the  optical 
nonlinearity  of  glass  films  with  various  compositions  prepared 
by  the  sol-gel  method.  The  samples  used  are  listed  in  Table 
7.1.  The  Ge02-Ti02  films  were  prepared  by  using  a procedure 
similar  to  that  for  fabricating  the  Si02-Ge02  system  as 
described  in  section  3.1.1,  except  that  titanium 
isopropoxide  (Ti  (OC3H7)4,  Aldrich)  and  germanium  2-methoxy- 
ethoxide (Ge (OCH2CH2OCH3) 4)  were  used  as  precursors. 

Both  Ge02  and  Ti02  are  important  oxide  components  used  in 
various  optical  materials  to  increase  the  refractive  index. 
By  varying  the  amount  of  Ge02  or  Ti02  added  to  Si02  we  may  span 
the  linear  optical  properties,  refractive  index  and  Abbe 
number,  from  Crown  glasses ( pd>50)  to  Flint  glasses ( j>d<50)  . 


6.2  Semi-Empirical  Model:  BGO  Equation 

The  most  commonly  used  model  for  predicting  the  nonlinear 
third  order  optical  susceptibility  of  transparent  optical 
glass  was  developed  by  Boling,  Glass  and  Owyoung(Bol78) . This 
model  assumes  a single  component  which  determines  linear  and 
nonlinear  indices.  For  oxide  glasses  this  component  is 
oxygen.  A semi-empirical  equation (i.e.  BGO  equation)  relating 
linear  refractive  index  and  dispersion  to  nonlinear  index  was 
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Table 

7.1  Sol-Gel  Thin 

Films  for  THG 

Measurement. 

Sample 

Temp.(°C) 

IndexOh) 

Thickness  (/xm) 

B30GS(22x) 

700,  02 

1.44 

1.327 

C50GS(20x) 

700,  02 

1.53 

0.76 

D30GT(22x) 

200,  02 

1.86 

0.25 

E50GT(22x) 

200,  02 

1.84 

0.45 
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derived  as  follows 


n2(10~,3esu)  = 


*(7iD-l)  (nD  + 2) 


[1.517  + (lip  + 2)  (nD  + 1)_L] 

6 nD 


(7.2.1) 


where  K=68,  j/D=Abbe  number  = (nD-l)  / (nF-nc)  , Xp=486nm,  X„=587nm, 
and  Xc=656nm.  The  equation  has  proven  accurate  in  predicting 
and  describing  the  NLO  coefficient  of  many  materials (Ada87) . 

Figure  7.1(Sim91)  shows  the  dependence  of  the  nonlinear 
index  on  linear  index  and  dispersion  calculated  from  the  BGO 
equation.  The  figure  clearly  shows  a strong  dependence  on 
dispersion,  indicative  of  the  effect  of  a nearby 
resonance (e.g.  possibly  near  1/2  Eg  or  Eg)  . 

In  order  to  calculate  n2  using  the  BGO  equation,  both  the 
linear  refractive  index (n,)  and  the  Abbe  number (j/D)  must  be 
known.  An  additive  model  developed  by  Huggins  and  Sun(Hug43) 
is  applied  to  compute  the  material  density,  linear  refractive 
index  as  well  as  Abbe  number  of  the  glass  films  we  investigate 
in  the  present  study. 

In  the  Huggins  method,  the  oxides  composing  a given  glass 
are  expressed  in  parts  by  weight (fM)  instead  of  weight 
percent (wt% ) , the  former  is  defined  as 
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Figure  7.1  Nonlinear  refractive  index  (n2)  versus  linear 
refractive  index (ni)  for  various  Abbe  numbers,  calculated 
using  BGO  equation [ Sim91 ] . 
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fM  = 


Pm 

100 


*E^  = 1 


(7.2.2) 


A part  by  weight  fM  is  therefore  equal  to  one  hundredth  of  the 
amount  of  oxide  expressed  in  weight  percent,  pM. 

The  density (d)  of  the  oxide  glass  of  an  arbitrary 
composition  may  be  calculated  from  the  following  equation 


4=E  (vm*m)  (7.2.3) 

where  vM  is  specific  volume  of  each  oxide  component.  In  the 
present  work  it  has  been  assumed  that  for  the  entire 
composition  range  of  binary  oxide  glass  (e.g.  Si02-Ge02,  Ge02- 
Ti02,  or  Si02-Ti02)  the  structural  units  remain  unchanged  and 
form  a solid  solution  as  has  been  described  in  section  4.2.5. 
Therefore  the  specific  volume  of  each  oxide  component  is 
assumed  to  be  constant  and  is  calculated  from  the  reported 
density  data  in  the  literature,  namely,  dsiO2=2.202  g/cm3, 
dOcO2=3.604  g/cm3[Huan78  ] , and  dTiQ2=3 . 8 g/cm3  [Vol88].  The  vM 
values  which  are  calculated  using  the  relation  vM=l/d  are 
given  in  Table  7.2. 

The  linear  refractive  index  n^  for  an  arbitrary 
wavelength,  X,  can  be  computed  from  the  following 
equation (Hug4 3) 
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Table  7.2  Various  Constants  Used  to  Calculate  the  Refractive 
Index  and  the  Abbe  Number  for  Si02-Ge02,  Ge02-Ti02, 
and  Si02-Ti02  Sol-Gel  Glass  Systems. 


Oxide 

vM 

eM*  10-8 

§m*  10-8 

Si02 

0.4541 

29.81 

145.69 

Ge02 

0.2775 

13.90 

85.99 

Ti02 

0.2632 

11.49 

39.6 

* After  Huggins  and  Sun(Hug43) 
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nX  = l+tf(ErAo/J 


(7.2.4) 


where  rMX  is  defined  for  most  of  the  oxides (except  PbO,  CaO, 
and  BaO)  as 


The  values  of  eMxl08  and  gMxl08  are  reproduced  from  Huggins 
and  Sun's  paper (Hug43)  and  listed  in  Table  7.2. 

The  Abbe  number (or  relative  dispersion)  of  the  glass  can 
be  obtained  by  using  the  following  equation 


where  nD=refractive  index  for  sodium  D line  at  589nm,  nF  and 
nc  are  index  for  Hydrogen  F and  C lines  at  486nm  and  656  nm 
respectively,  which  can  be  computed  using  equation  (7.2.4). 

Figure  7.2  shows  the  calculated  linear  refractive  index 
at  632 . 8nm (HeNe  laser  wavelength  used  in  the  ellipsometry 
measurements)  by  using  the  additive  model  for  three  binary 
oxide  glass  systems(e.g.  Si02-Ge02,  Ge02-Ti02,  and  Si02-Ti02)  of 
various  compositions.  The  curves  are  nearly  linear  for  these 


v.=e»( — i-_  -48xl0-6X!] 

9“'V 


(7.2.5) 


(7.2.6) 
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Ge02(wt%) 


Figure  7.2  Calculated  refractive  indices  for  various 
binary  oxide  glass  systems  using  an  additive  model;  □:  data 
for  Ge02-Ti02  THG  samples,  *:  data  for  Si02-Ge02  THG  samples. 
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binary  systems  with  an  index  of  1.457  for  pure  glassy  Si02, 
1.60  for  pure  Ge02  glass,  and  2.176  for  pure  Ti02  in  a glassy 
state.  Titania  has  the  smallest  specific  volume  among  the 
oxides  calculated,  which  results  in  a rapid  increase  of  the 
density  and  therefore  an  increase  in  the  linear  refractive 
index  with  the  addition  of  Ti02  in  the  material. 

The  experimentally  measured  refractive  index  of  the  sol- 
gel  binary  glass  film  used  in  the  THG  measurements  are  shown 
in  Figure  7.2  as  asterisk  symbols  for  Ge02-Si02  samples  and  as 
hollow  squares  for  30GT  and  50GT  films.  It  appears  that  only 
50GS  film  is  close  to  the  theoretical  value,  while  the  other 
samples  exhibit  indices  lower  than  the  predicted  values, 
indicating  that  these  samples  (30GS,  30GT  and  50GT)  may  not  be 
fully  densified.  This  is  not  unexpected  considering  that  the 
heat  treatment  temperatures  for  these  films  are  not  high, 
especially  for  Ge02-Ti02  samples  (only  heated  at  200°C)  . 

The  computed  results  of  Abbe  numbers  of  these  binary 
oxides  are  shown  in  Figure  7.3.  Titania  has  the  largest 
dispersion  which  leads  to  a small  Abbe  number.  The  addition 
of  Ti02  into  the  material  causes  a rapid  decrease  of  Abbe 
number.  The  computed  Abbe  number  for  pure  glassy  oxides  are 
^sio2=  67.86,  ^0e02=41.88,  and  ^Tio2=18.78. 

6.3  Third  Harmonic  Generation (THG)  Measurement 

As  mentioned  in  section  3.2.11,  a Third-Harmonic- 
Generation  (THG)  interference  fringe  method  was  used  to 
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Ge02(wt%) 


Figure  7.3  Calculated  Abbe  numbers  for  various 

oxide  glass  systems  using  an  additive  model. 
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measure  the  nonlinearity  of  the  films.  The  measurement  and 
the  analysis  of  the  THG  fringes  were  conducted  by  Dr.  W.E. 
Torruellas  in  Dr.  G.  Stegeman's  group  at  the  University  of 
Arizona.  To  optimally  use  THG  method  it  is  desirable  to  have 
the  product  of  the  THG  coherence  length  (lc)  of  the  substrate 
and  the  substrate  nonlinearity  be  comparable  to  the  product  of 
the  film  thickness  and  its  nonlinear  susceptibility (x®) • The 
coherence  length  in  this  case  means  the  distance  that  the 
fundamental  wave  and  the  third  harmonic  wave  may  travel 
without  coming  out  of  phase,  as  the  beginning  of  destructive 
interference  diminishes  the  amplitude  of  the  TH  wave.  Since 
the  THG  coherence  length  of  fused  silica  at  1.579/xm  is 
approximately  15/im[Tor91a] , the  third  harmonic  optical  fields 
for  the  film  and  the  substrate  are  comparable  for  -1/xm  films. 
Both  the  magnitude  and  the  phase  of  x(3)  could  be  measured  from 
the  shift  and  the  offset  of  the  interference  fringes [Tor91a] , 
an  example  of  which  for  sample  B30GS  is  shown  in  Figure  7.4. 

Figure  7.5  shows  the  THG  results  measured  at  1.579jum  for 
various  sol-gel  film  compositions.  The  data  for  Si02-Ti02  are 
taken  from  Torruellas  et  al.  's  recent  work[Tor91b] . As  can  be 
seen,  the  results  of  the  Si02-Ge02  and  Si02-Ti02  films  follows 
Boling's  model  appropriately.  All  the  n2  values  measured  are 
within  the  region  between  the  curves  corresponding  to  vD=20 
and  vd=60,  which  is  the  typical  range  of  Abbe  numbers  for 
these  binary  oxide  glasses. 

The  results  of  the  calculated  nonlinear  refractive 


T H G s i g n a 


216 


Figure  7.4  THG  interference  fringes  of  30GS  film  at  1.579 
/im.  : the  fringe  of  bare  Si02  substrate,  ■:  the  experi- 
mental data  points,  and  the  solid  line  is  the  best  fit  to  the 
interference  fringe [Tor91 ] . 
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Figure  7.5  n2  is  plotted  as  a function  of  n,.  Solid  line: 
calculated  values,  markers: the  experimental  data,  n,  for  SiO,- 
Ti02  are  calculated  using  the  additive  model. 


218 


index  (n2B)  using  the  BGO  equation  and  that  from  THG 
measurements  for  1.579/im  (n2THO)  are  given  in  Table  7.3.  The 
linear  refractive  index  and  Abbe  number  of  each  film  are 
computed  using  the  additive  model  described  previously.  The 
THG  data  for  Si02(i.e.  0.85xl0'13  esu)  is  taken  from  Weber  et 
al.'s  results  measured  at  1 . 064jLtm(Web78)  . The  results  of  the 
THG/ Si02  ratio  are  taken  from  Torruellas's  work [Tor 9 la ] . The 
comparison  between  the  calculated  n2  values  (n2B)  and  the 
measured  nonlinear  indices  (n2THO)  are  plotted  as  shown  in 
Figure  7.6.  Although  the  measured  values  seems  to  follow  the 
predicted  values  quite  closely,  in  most  of  the  cases  the 
measured  values  are  slightly  higher  than  the  calculated 
values,  except  30GS  and  two  Ge02-Ti02  samples.  The  fact  that 
both  Ge02-Ti02  films  possess  much  lower  values  compared  to  the 
predicted  ones  is  attributed  to  the  possibility  that  these 
films  are  still  quite  porous,  since  they  were  only  heat 
treated  at  200  °C  to  prevent  crystallization  of  the  film.  The 
measured  values  of  n2  for  the  Si02-Ti02  system  appears  to  be 
slightly  higher  than  the  predicted  values  by  the  Boling  model. 
This  can  be  due  to  the  fact  that  the  model  assumes  oxygen  ions 
as  the  only  dominant  constituent  for  both  linear  and  nonlinear 
indices.  When  the  material  contains  heavy  ions,  especially 
ions  with  unfilled  d— shells,  the  d— electrons  will  make  a 
significant  contribution  to  the  linear  and  nonlinear 
polarizabilities [Bol78 ] . Therefore,  the  addition  of  titania 
to  the  glass  may  cause  some  deviation  from  the  prediction  of 
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Table  7.3  Comparison  Between  the  Calculated  Nonlinear  Refractive 
Index  (n2B)  Using  BGO  Equation  and  from  THG  Measurement 

( n2THG ) • 


Si02 

(wt%) 

Ge02 

Ti02 

nd 

u* 

n2B(esu) 

*THG/Si02 

100 

0 

0 

1.459 

67.86 

0.87x1 0'13 

1 

0.85x10'13 

57.3 

42.7 

0 

1.523 

56.54 

1 .4x1  O'13 

2 

1.7x1  O'13 

36.5 

63.5 

0 

1.537 

49.15 

1 .8x1  O’13 

4 

3.4x1 013 

0 

56.7 

43.3 

1.868 

23.79 

13x1  O'13 

7 

6.0x10'13 

0 

36.0 

64.0 

1.982 

20.5 

20x1  O'13 

15 

1 2.7x1  O'3 

93.5 

0 

6.5 

1.474 

61.49 

l.lxlO13 

5 

4.2x1  O'13 

63.7 

0 

36.3 

1.625 

32.01 

4.4x1  O'13 

9 

7.6x1 013 

42.9 

0 

57.1 

1.748 

25.75 

8.5x1’3 

12 

1 0.2x1 013 

24.4 

0 

75.6 

1.918 

21.74 

1 6x1  O'13 

21 

1 7.8x1  O'13 

* Data  reported  by  W.E.  Torruellas.  The  Si02-Ti02  samples  were  prepared  at  the 
University  of  Arizona(Z.  Osborne  and  B.J.  Zelinski) 

@ Data  calculated  by  J.H.  Simmons. (Sim90) 
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Figure  7.6  Comparison  between  the  predicted  n2  values 
according  to  BGO  model  and  that  measured  by  THG  measurement. 
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the  model. 

Figure  7.7[Tor91a]  shows  the  relative  magnitude  of  x0)(3(i>) 
compared  to  pure  fused  silica  at  three  different  wavelengths 
in  the  near  infrared  for  50GT  and  50ST  films.  In  the  case  of 
the  Si02-Ti02  sample  (whose  bandedge  is  at  ~350nm)  the  value  for 
XQ)  remains  nearly  constant  until  the  frequency  approaches  a 
resonance  at  3(t)09Si.  Figure  7.8  shows  calculated  dispersion 
curves  for  the  THG  samples  using  the  additive  model.  As  can 
be  seen  in  the  Figure  7.8,  all  the  Ti02-containing  glasses 
start  to  approach  resonance  at  around  350nm  while  Si02-Ge02 
glasses  resonate  at  frequencies  below  200nm (unless  some  defect 
bands  exist  in  the  material  e.g.  245nm,  325nm  etc.).  This  can 
also  be  seen  in  the  UV-VIS  absorption  spectra  in  Figure  7.9 
for  some  of  the  films.  The  resonance  enhancement  on  the 
measured  xa)  has  been  observed  in  various  measurements  on 
organic  f ilms [Kaj83 ] . 

All  the  other  samples  containing  Ge02  in  this  study 
possess  similar  dispersion  as  film  50GT  as  are  shown  in  Figure 
7.10.  Due  to  the  fact  that  only  a limited  number  of 
wavelengths  were  measured  on  these  films,  the  reason (s)  for 
the  drop  of  x( ) 0.953/xm  for  all  the  Ge02— containing  samples 
is  not  known.  It  has  been  observed  [Tor91a]  that  the  Ge02 
based  samples  have  a stronger  tendency  to  damage  at  0.957/im 
than  the  Si02-Ti02  ones.  The  pump  energy  used  for  those 
samples  thus  was  near  the  detection  limit  of  the  THG 
measurement  at  that  wavelength.  Therefore  the  error  of  the 
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Figure  7.7  The  relative  magnitude  of  x(3)(3cj)  compared  to 
pure  fused-Si02  is  plotted  as  a function  of  the  fundamental 
wavelength  for  50GT  and  50ST  samples,  at  0.957jum,  1.064um  and 
1. 597/im[Tor91]  . 
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Figure  7.8  Calculated  dispersion  curves  of  various  samples 
used  in  THG  measurement  using  an  additive  model. 
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LINEAR  ABSORBANCE  normalized  to  lOOOnm 


Figure  7.9  The  linear  absorbance  for  various  THG  samples 
measured  by  UV-Vis  spectroscopy,  spectra  are  normalized  to  the 
same  thickness [Tor 91c] . 
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Figure  7.10  The  relative  magnitude  of  x<3)(3w) 
pure  fused  silica  is  plotted  as  a function  of  the 
wavelength  for  various  THG  samples [Tor91c] . 
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measurement  at  that  wavelength  may  be  large.  If  the  drop  of 
XP)  at  0.957  /xm  for  Ge02-containing  samples  is  real,  we 
speculate  that  it  is  possible  that  the  THG  process  is 
exhibiting  a similar  nonlinear  index  oscillation  near  the 
resonance  frequencies  as  seen  in  a two-photon-absorption 
process (see  Figure  7 . ll [Sim91 ] ) , and  we  are  measuring  the  dip 
in  that  dispersion  curve.  However,  more  complete  dispersion 
studies  need  to  be  conducted  in  order  to  draw  a final 
conclusion. 


Non-Linear  Index  (nano-esu) 
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Figure  7.11  Plot  of  the  two-photon-absorption  induced 
nonlinearity  for  materials  with  different  bandgap  and  impurity 
absorption  energies  [Sim91]. 


CHAPTER  VIII 


SUMMARY  AND  CONCLUSIONS 


The  present  study  has  for  the  first  time  successfully 
demonstrated  the  feasibility  of  employing  both  sol-gel 
multiple  dip-coating  and  low  pressure  chemical  vapor 
deposition (LPCVD) , in  the  production  of  high  quality  Ge02-Si02 
thin  films  with  various  germania  contents.  Planar  view 
transmission  electron  microscopy  has  shown  that  the 
microstructure  of  the  sol-gel  film  consists  of  very  fine 
particles  in  the  size  range  of  lOnm,  while  the  LPCVD  films  are 
composed  of  small  particulates  of  around  30nm  to  lOOnm  size 
range  depending  upon  the  deposition  condition. 

In  the  sol-gel  technique,  several  germanium  alkoxides 
with  higher  alkyl  groups  have  been  synthesized  from  GeCl4 
using  a combination  of  an  ammonia  method  and  an  alcohol- 
interchange  technique.  Germanium  2-methoxy-ethoxide  was 
routinely  used  in  the  current  study  to  prepare  stable  binary 
Ge02-Si02  sols  for  film  deposition,  with  Ge02  content  between 
10mol%  to  50  mol% . A small  amount  of  hydroxy— propyl- 
cellulose  (HPC)  was  added  to  the  sols  as  an  independent 
viscosity  adjusting  agent  to  facilitate  the  multiple  dipping 
process.  This  resulted  in  films  of  improved  thickness 
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uniformity. 

According  to  the  thermal  analysis  of  the  organic  additive 
in  Ge02-Si02  gels,  the  organics  may  be  removed  during  heat- 
treatment  in  a temperature  range  of  200  to  600°C. 
Densif ication  behavior  of  the  films  was  followed  by  monitoring 
the  change  in  refractive  index  as  a function  of  the  heat 
treatment  temperature.  All  of  the  sol-gel  samples  exhibited 
an  index  decrease  in  the  temperature  range  between  200  to 
600°C  due  to  the  removal  of  the  organics.  This  is  then 
followed  by  an  increase  in  the  refractive  index  indicating 
densif ication  of  the  film.  The  densif ication  temperature 
appears  to  increase  with  decreasing  Ge02  content.  The 
refractive  index  difference  between  the  sol-gel  films  and  the 
Si02  substrate  may  be  tailored  by  controlling  the  film 
composition  and  the  extent  of  densif ication.  The  refractive 
indices  of  the  Ge02-Si02  films  span  from  1.458  to  1.61, 
corresponding  to  pure  Si02  and  pure  Ge02  glass,  respectively. 
Its  composition  dependence  follows  that  predicted  by  the 
Lorentz-Lorenz  model. 

Both  single  and  multimode  waveguides  have  been  fabricated 
by  the  sol-gel  process.  Using  a three-prism  loss  measurement 
technique,  the  propagation  loss  of  the  50Ge02-50Si02  single 
mode  waveguide  was  found  to  be  3.31  dB/cm.  These  losses  are 
attributed  to  various  scattering  processes,  caused  by  surface 
roughness  and  pinholes  in  the  film. 

Vibrational  spectroscopy  revealed  the  existence  of  Si-O- 
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Ge  linkages  in  Ge02-Si02  glass  network.  The  addition  of  Ge02 
in  Si02  caused  a decrease  in  the  size  of  both  the  D1  and  D2 
defect  bands  in  the  Si02  Raman  spectra.  A preferential 
decrease  of  the  800  cm'1  absorption  band  in  the  FTIR  spectra, 
due  to  the  Si-O-Si  stretching  mode,  was  also  observed. 

In  the  LPCVD  process,  dilute  GeH4/Ar  and  pure  silane (SiH4) 
gasses  were  used  as  reactants.  The  Ge02-Si02  film  was 
deposited  by  an  oxidation  reaction  with  02  in  the  temperature 
range  of  400-450°C  . Appropriate  film  deposition  conditions 
were  found  to  produce  waveguiding  films  using  reasonable 
deposition  rates.  The  Ge02  content  of  the  films  prepared 
using  this  technique  were  limited  due  to  the  use  of  a very 
dilute  GeH4  reactant  and  an  inherently  low  deposition  rate  of 
Ge02  in  this  process. 

Analysis  of  the  UV-Vis  spectra  revealed  that  these  films 
possess  a high  concentration  of  oxygen  deficient  centers 
compared  to  fused  Si02.  This  has  been  extensively  reported 
for  oxide  films  deposited  under  vacuum. 

Vibrational  spectroscopic  studies  confirmed  the  existence 
of  Si-O-Ge  linkages  in  the  binary  glass  network.  The 
structure  of  the  LPCVD  films  appears  to  be  dominated  by  D1  and 
D2  defect  structures,  when  comparing  the  Raman  spectra  of  the 
LP3  film  to  those  of  Si02  after  neutron  irradiation.  This  is 
sttributed  mainly  to  the  low  deposition  temperature  used  in 
the  process. 

The  propagation  loss  of  the  LP3  waveguide  is  2.59  dB/cm. 
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The  scattering  processes  caused  by  surface  roughness  and  some 
large  particulate  inclusions  in  the  film  are  believed  to  be 
responsible  for  this  loss. 

The  mode  indices  of  the  waveguides  were  measured  using  a 
prism  coupling  technique.  In  comparing  the  experimentally 
measured  mode  index  to  that  calculated  based  upon  a step-index 
profile  assumption,  good  agreement  was  found.  Thus  the 
refractive  index  profile  of  the  films  prepared  by  both  the 
sol-gel  and  the  LPCVD  process  may  be  approximated  by  a step- 
index  profile. 

The  theoretical  electromagnetic  field  distribution 
profile  for  a step-index  planar  waveguide  has  been  calculated. 
The  propagating  mode  profile  in  the  Ge02-Si02  waveguide  is 
experimentally  measured,  using  a near  field  technique,  and 
compared  to  the  computed  mode  profiles.  The  analysis  revealed 
that  the  measured  mode  profiles  substantiate  the  step-index 
planar  waveguide  model.  However,  slight  discrepancies  in  peak 
heights  between  the  calculated  and  experimentally  determined 
profiles  have  been  observed.  The  deviations  in  the  mode 
intensity  may  be  caused  by  the  surface  scattering  at  the  end 
face  of  the  waveguide  due  to  defects  introduced  during 
polishing  and/or  possible  small  index  fluctuation  across  the 
deposited  film. 

The  nonlinear  refractive  index  of  the  sol-gel  films 
(Ge02-Si02,  and  Ge02-Ti02)  was  measured  using  a THG 
interferometry  fringe  technique  in  cooperation  with  Dr.  W.E. 
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Torruellas  and  Dr.  G.  Stegeman  at  the  University  of  Arizona. 
Analysis  of  both  THG(n2THO)  and  linear  index (n,)  measurement 
results  shows  that  the  relation  between  n^o  and  nt  follows 
that  predicted  by  the  empirical  BGO  equation. 

The  BGO  model  has  been  used  to  predict  the  nonlinear 
refractive  indices  (n2B)  of  various  binary  oxide  glasses  (i.e. 
Ge02-Si02/  Ge02-Ti02,  and  Si02-Ti02)  . In  order  to  compute  n2B/ 
an  additive  model  was  used  to  calculate  linear  refractive 
indices  and  Abbe  numbers  of  the  materials,  assuming  a constant 
specific  volume  for  each  individual  oxide  component  (i.e. 
Si02,  Ge02,  and  Ti02)  . The  deviation  of  the  n2THO  values  of 
Ti02-containing  glasses  from  the  n2B  values  is  attributed  to 
the  high  polarizability  of  Ti  cations. 

Based  upon  the  limited  measurement  of  the  x®  dispersion 
in  these  binary  films,  the  Si02-Ti02  films  show  a sharp 
increase  of  x®  in  the  low  wavelength  side,  which  is  attributed 
to  a resonant  enhancement  of  the  third  harmonic  wave.  This 
agrees  with  the  calculated  n,  dispersion  curves  using  the 
additive  model.  It  is  also  confirmed  by  the  UV-Vis 
measurement  which  showed  a large  absorption  band  starting  at 
around  350nm. 

The  unusual  decrease  of  x®  for  both  Ge02-Si02  and  Ge02-Ti02 
systems  on  the  short  wavelength  region  is  not  understood. 
However,  the  THG  process  is  speculated  to  exhibit  a similar  n2 
oscillation  near  resonance  frequencies  similar  to  that  seen  in 
the  Two-Photon-Absorption  process.  Thus,  the  approach  to 
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resonance  can  cause  the  measured  dip  in  the  n2  dispersion 
curve.  However,  a more  systematic  dispersion  study  is 
required  before  a final  conclusion  can  be  drawn. 

In  conclusion,  it  is  possible  to  form  suitable  optical 
waveguides  in  the  Ge02-Si02  system  by  both  sol-gel  multi- 
coating and  LPCVD  approaches.  The  sol-gel  process  offers  a 
broad  range  of  compositions  and  therefore  numerical  apertures 
while  the  LPCVD  is  limited  to  low  Ge02  concentrations.  We 
have  demonstrated  that  the  sol-gel  multiple  coating  process 
with  up  to  72  layers  can  be  used  to  form  films  of  reasonable 
uniformity  when  using  an  organic  additive  to  control  solution 
viscosity  during  dipping.  The  losses  measured  in  waveguides 
made  by  the  sol-gel  process  are  mostly  due  to  scattering  and 
could  be  reduced  by  (1)  processing  under  clean  room  conditions 
and  (2)  using  a highly  automated  dip-coating  process  to 
maintain  uniform  thickness  over  the  film.  The  losses  in  the 
LPCVD  waveguide  appear  to  result  from  scattering  losses  also, 
although  hydroxyl  content  may  also  be  a problem.  Both  may  be 
reduced  by  high  temperature  anneals. 
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